我一开始的想法是分别利用龙格库塔法解每一个微分方程,但变量很多,算法会比较复杂。后来明白可以把多变量看作是一个变量,利用matlab的feval函数进行代入变量的函数运算。 matlab中feval函数的作用:feval(f,x,y);将x,y代入函数f中。 四阶龙格-库塔法: 需要解的四个微分方程组为: 算法代码: 主函数main.m main....
求matlab编程源代码用四阶龙格库塔法解如下微分方程组dx1tdtx2tdx2tdtx3tdx3tdt60x1t25x2t10x3t1初始条件x10x20x300求解区间05h01functionrrk4fabyan为步长ya为初值hban 四阶龙格库塔法解如下微分方程组的MATLAB编程 求MATLAB编程源代码用四阶龙格库塔法解如下微分方程组 dx1(t)/dt=x2(t) dx2(t)/dt=x3(t...
matlab用四阶龙格库塔法解微分方程组;dy(1)=-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));dy(2)=(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);以上述计算结果再算一下表达式;x=3000*sin(70*2*pi/360)-y(1)*sin(y(2));...
*sin(y(2)); z=3000*cos(70*2*pi/360)-12*t-y(1).*cos(y(2)); plot(t,x,t,z)最后生成的图像效果不太理想。你可能需要重新调整参数或改进绘图方式。不过，所用的四阶龙格-库塔法是正确的。建议你可以在MATLAB中文论坛上提问，那里有很多高手可以提供帮助。如果你遇到的是绘图方面的问题...
matlab用四阶龙格库塔法解微分方程组;dy(1)=-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));dy(2)=(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);以上述计算结果再算一下表达式;x=3000*sin(70*2*pi/360)-y(1)*sin(y(2));z=3000*cos(70*2*pi/360)-12*t-y(1)*co...
matlab数值分析 常微分方程初边值问题数值解 标准龙格库塔四阶四段公式 欧拉法 1、标准龙格库塔四阶四段公式 function y=rk4(fun,a,b,y0,n) h=(b-a)/n; y(1)=y0; for k=1:n x=a+(k-1)*h; k1=h*feval(fun,x,y(k)); k2=h*feval(fun,x+h/2,y(k)+k1/2); ...
functiondy=test(t,y) dy=[-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));(12*sin(...
matlab用四阶龙格库塔法解微分方程组;dy(1)=-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));dy(2)=(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);以上述计算结果再算一下表达式;x=3000*sin(70*2*pi/360)-y(1)*sin(y(2));z=3000*cos(70*2*pi/360)-12*t-y(1)*co...
利用matlab求一阶变系数微分方程组的数值解。 已经有7人回复 急!求大神用Matlab四阶龙格库塔解个方程! 已经有5人回复 请问什么是分段龙格库塔的方法?需要用matlab软件采用这种方法求二阶微分方程的解 已经有3人回复 mathcad怎样用四阶龙格库塔法求解二阶微分方程? 已经有8人回复 求解一个初中的二元四次方程组,都...
利用matlab求一阶变系数微分方程组的数值解。 已经有7人回复 急!求大神用Matlab四阶龙格库塔解个方程! 已经有5人回复 请问什么是分段龙格库塔的方法?需要用matlab软件采用这种方法求二阶微分方程的解 已经有3人回复 mathcad怎样用四阶龙格库塔法求解二阶微分方程? 已经有8人回复 求解一个初中的二元四次方程组,都...