龙格库塔解方程组
答:龙格库塔法是求解微分方程或微分方程组的。如果你的问题指的是齐次线性微分方程组,那是能够求解的。如果问的是线性代数里面的线性齐次方程组,恕我直言,你得好好看书了,不是一个领域内的问题。
答:求解二阶微分方程,初始条件还需要给出y1'(0)和y2'(0)。这里暂时按照0处理。function zd530003514 a=0.1;b=0.1;Y0 = [b-1; 0; b; 0];解方程 [t,Y]= ode45(@ode,[0 10],Y0);y1=Y(:,1);y2=Y(:,3);绘图 subplot 211 plot(t,y1);subplot 212 plot(t,y2);微分方程定义...
答:分析了题主的问题,1、缺初始条件(如,x、x'、y、y'、xb、xb'、yb、yb');2、给出的单位有点乱(不统一);3、k(y-yb)与ks((x1-xb)^2+(y1-yb)^2)*(x-xb))的单位不一致 所以,在用龙格库塔法求解时,遇到无法接受的结果。请题主认真核查一遍。把正确的问题贴出来,以便帮助你。
答:四阶龙格-库塔是求解常微分方程(常微分方程组)精度最高的一种数值方法。四阶龙格-库塔迭代公式为 根据四阶龙格-库塔迭代公式,可用matlab软件编写得到如下计算代码:调用方法:clc;close all;y0=1;h=0.1;a=0;b=1;[x,y] = runge_kutta(@func,y0,h,a,b);disp(' x y')[x',y']...
答:'用四阶龙格库塔法求解一阶常微分方程组 Function longkuta(a() As Integer, b() As Integer) As Single Dim X() As Single, x1() As Single, z() As Single, e() As Double Dim i As Integer, j As Integer ReDim Preserve a(0 To n)ReDim Preserve b(0 To n)ReDim X(1 To...
答:……所谓龙格库塔法,通俗地说,就是把一个n阶的常微分方程,整理成n个形如 f'(t)=g(t,f(t)) (注意此时右侧不含 f(t) 的导数)的一阶常微分方程组再加以求解的方法。你的方程整理成龙格库塔所需要的形式就是:x'[t] = y[t]y'[t] = (-c y[t] - k x[t] - F[t])/m ...
答:解方程很简单,调用ode45即可 关键是求运动时间,假设z方向坐标变为0即为落至地面运动结束,可求得运动时间的近似值及落地点坐标 代码如下:clear all clc c=0.07;u0=25;a=pi/6;f=@(t,x)([x(2);-c*sqrt(x(2)^2+(x(4)-7-.35*x(5))^2+x(6)^2)*x(2);x(4);-c*sqrt(x...
答:printf("用二阶龙格-库塔方法\n");RungeKutta(1,0,1,5,x,y,2,function);for(int i=0;i<6;i++)printf("x[%d]=%f,y[%d]=%f\n",i,x[i],i,y[i]);printf("用三阶龙格-库塔方法\n");RungeKutta(1,0,1,5,x,y,3,function);for(i=0;i<6;i++)printf("x[%d]=%f,y...
答:c1=t*(-y(i-1));c2=t*(-y(i-1)-c1/2);c3=t*(-y(i-1)-c2/2);c4=t*(-y(i-1)-c3);注意加减号问题。RK的精度还是比较高的。
答:程序开始,定义求解区间与步长 clear;clc;h=0.1;% 步长 a=0;% 求解区间下限 b=1;% 求解区间上限 x=a:h:b;% 变量取值 y_analy=sqrt(1+2*x);% 解析解 plot(x,y_analy,'b');hold on; %绘制精确解曲线 n=(b-a)/h;y0=1;Runge-Kutta 2 二阶龙格-库塔法 y_rk2(1)=y0;for ...
网友评论:
浦券18056475694:
matlab用四阶龙格库塔法解微分方程组; -
56878宰史
: function dy=test(t,y) dy=[-12*cos(y(2))-120*cos(208*2*pi/360-3*y(2));(12*sin(y(2))+120*sin(208*2*pi/360-3*y(2)))/y(1);];[t,y]=ode45('test',[0.01,1],[1,1])plot(t,y(:,1),t,y(:,2)); x=3000*sin(70*2*pi/360)-y(1).*sin(y(2)); z=3000*cos(70*2*pi/360)-...
浦券18056475694:
请问用四阶龙格库塔法解二阶微分方程的思想是什么?最近我遇到了一个难题!就是求一个二阶微分方程,形式如:y''=f(y),初始条件是y(0)=0,y'(0)=0,y是t的... -
56878宰史
:[答案] 令y1=y';y2=y''=(y1)';得到了二阶微分方程的一阶形式,然后可以按照普通的一节微分方程组的四阶龙格库塔法求解了;这样说的有点抽象,你可以把方程复制出来 我有空帮你编写一下
浦券18056475694:
急!!!求matlab 用四阶龙格 - 库塔法求解常微分方程 -
56878宰史
: 建立.m文件 --------------------------------------------- function theta=danbai(t,X) x=X(1); dx=X(2); ddx=-sin(x); theta=[dx;ddx]; ---------------------------------------------- 命令窗口输入 >> [t,Y]=ode45(@danbai,[0 6],[pi/3 -1/2]); >> plot(t,Y(:,1),'ro-',t,Y(:,2),'bv-'); >> legend...
浦券18056475694:
Matlab用四阶龙格库塔法求解不可微分方程组的初值u'=x+u+v u(0)=0 0 -
56878宰史
:[答案] fun=@(x,uvw)([x+uvw(1)+uvw(2);-x+uvw(1);uvw(1)+uvw(3)]);[x,uvw]=ode45(fun,[0,1],[0 1 1]);plot(x,uvw)legend('u','v','w')%下面是相图figureplot3(uvw(:,1),uvw(:,2),uvw(:,3))
浦券18056475694:
用四阶龙格库塔法求解矩阵微分方程
56878宰史
: global R M U syms t R=[ 0.0247,0,0,0,0,0; 0,0.0247,0,0,0,0; 0,0,0.0247,0,0,0; 0,0,0,0.0193,-0.0193,0; 0,0,0,0,0.0193,-0.0193; 0,0,0,1,1,1 ]; M=[ 0,0,0,-15727/10000*sin(5/12*pi+80*pi*t)*pi,15727/10000*cos(1/4*pi+80*pi*t)*pi,15727/10000*sin(1/12*...
浦券18056475694:
龙格库塔法能求解线性齐次方程组吗 -
56878宰史
: 龙格库塔法是求解微分方程或微分方程组的.如果你的问题指的是齐次线性微分方程组,那是能够求解的.如果问的是线性代数里面的线性齐次方程组,恕我直言,你得好好看书了,不是一个领域内的问题.
浦券18056475694:
用龙格库塔法求解微分方程 -
56878宰史
: 功能:用四阶 Runge-Kutta 法求解常微分方程 --------------------------------------------- function R=Rungkuta4(f, a, b, n, ya)% f:微分方程右端函数句柄% a,b:自变量取值区间的两个端点% n:区间等分的个数% ya:函数初值 y(a)% R=[x',y']:自变量...
浦券18056475694:
急急急!求助matlab用龙格 - 库塔方法求解方程组 -
56878宰史
: function df=ode45_fun(t,xyzuvw) %%注意小写的v和大写的V %常数(请修正) R_0=1; rho_0=1; beta=1; G=6.67e-11; M=1.5e24;x=xyzuvw(1); y=xyzuvw(2); z=xyzuvw(3); u=xyzuvw(4); v=xyzuvw(5); w=xyzuvw(6);R=sqrt(x*x+y*y+z*z); V=sqrt(u*u+v*v+w...
浦券18056475694:
Matlab用四阶龙格库塔法求解不可微分方程组的初值y1'=120 - 2*y1+2*y2 (0 -
56878宰史
:[答案] %% function dx=aaa(t,x)%定义matlab函数来描述方程 dx=[120-2*x(1)+2*x(2);2*x(1)-5*x(2)] %% x0=[0;0] [t,y]=ode45('aaa',[0,1],x0); plot(t,y)
浦券18056475694:
跪求用龙格库塔求解微分方程的matlab程序 -
56878宰史
: >> f=inline('1/6-y/30','t','y');>> [t,y]=ode45(f,[0,5],[0]);>> plot(t,y) 另外,由高数的知识,可以得到解析解为:y=5-5*e^(-t/30) 把上述的解析值也放到上图中.>> hold on>> plot(t,5-5*exp(-t/30),'r*') 可以看到二者相当吻合.