作业帮 > 综合 > 作业

非线性微分方程的求解方程组-g*sinθ-k*v^2=dv/dt___(1)g*cosθ=v*dθ/dt________

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/07/17 16:46:55
非线性微分方程的求解
方程组
-g*sinθ-k*v^2=dv/dt___(1)
g*cosθ=v*dθ/dt________(2)
初值条件
v=v0
θ=θ0
g和k均为常数
解析法或数值法均可,数值解法最好给出matlab或C语言程序
知道dsolve函数就好求常微分方程或方程组了:)
>> s=dsolve('Dv=-k*v-g*sin(a),Da=g*cos(a)/v','v(0)=v0,a(0)=a0');
>> a=s.a
a =(-g*cos(a0)+a0*v0*k)/v0/k+g*cos(a0)/v0/k*exp(k*t)
>> v=s.v
v =1/cos(a0)*v0/exp(k*t)*cos((-g*cos(a0)+a0*v0*k)/v0/k+g*cos(a0)/v0/k*exp(k*t))