00问答网
所有问题
当前搜索:
matlab求一阶微分方程
如何用
matlab
解线性
微分方程
答:
把X表达成连续变量,则 f(x)=(
1
/2)δ(x)+(1/2)δ(x-1) --- 这里δ是脉冲函数.f(y)=u(y)-u(y-1) --- 这里u是阶跃函数.X,Y 独立,Z=X+Y,所以f(z)是f(x)和f(y)的卷积.f(z)=f(x)*f(y) --- * 是卷积. x处要带入z. y处也要带入z.f(z)=(1/2)(δ(z)+...
matlab
求解
微分方程
数值解 求代码
答:
结果:代码:clear allclc f=@(x,y)([y(2); y(
1
)/3536.4*(y(2)^2+14.142*y(2)+70.71)]);[x,Y]=ode45(f,[0 20],[3 0]);plot(x,Y(:,1),x,Y(:,2)),legend('y','dy/dt')xlabel('x')
用
Matlab
求解
微分方程
(组)
答:
这是阿波罗卫星轨迹例子 第二个
方程
应为y' = -2*x' + y ...,以为符号 代码如下 function solveode function fun = OdeFunc(t,y)mu =
1
/82.45;mu_star = 1 - mu;r2 = sqrt((y(1) - mu_star)^2+y(3)^2);r1 = sqrt((y(1) + mu)^2+y(3)^2);fun = [y(2);2*...
如何用
matlab
求解
微分方程
组的数值解
答:
在
Matlab
下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 function y=zhidao_feiying(t,x)a=1;b=2;c=3;d=4;比方说 f=a*x+y-z;g=b*sin(x*y)-c*cos(z);h=d*y-a*x;注意x用x(1)代,y用x(2)代,z用x(3)代 f=a*x(1)+x(2)-x(3);g...
matlab
求解
微分方程
的通解问题
答:
先分析,先把矩阵算一下,化简整理下 x‘-4t^2x=0 其中x是t的函数 编程如下;dsolve('dy-4*t^2*y=0')因为
matlab
中默认的dy表示y对t的
一阶
导,d2y表示y对t的二阶导,没有dx表示x对t的一阶导 所以会这样编程 运行结果:c2*exp((4*t^3)/3)+c1 ...
怎么用
matlab
解这个
微分方程
答:
求解常
微分方程
(组),一般可以用dsolve函数求其解析值,但有的微分方程只能通过ode函数求出其数值解。题主给出的微分方程,如用dsolve函数求解,则会出现错误,所以该方程可以用ode45函数得到其数值解。求解过程如下:a0=0;tspan=[0,5];fun=@(t,a)2/(-3+2*(
1
-a)^(-1/3)+4*(1-a)^(-...
用
matlab求微分方程
解析解,题目如下:
答:
给你参考一下我当年做的一个类似的程序 clear dt=0.004;t=0:dt:5;a=-pi:0.2:pi;l=100;r=8;d=10; %导弹尺度参数:长度、半径、火焰尺度 L=l;D=l/3;c=pi/4;R=L*tan(c);%飞机外形尺度参数:长度、尾翼宽、翼展(半)角度、翼(半)展 size(t);n=ans(
1
,2);for J=1:n ...
MATLAB
微分方程
组求解
答:
除了有参数,还要有x,y,x',y'的初始值 看你的
方程
,类似有阻力的抛体运动 这里随便给了个初值,x(0)=0 y(0)=10 x'(0)=10 y(0)'=0 可以理解为物体从10米高的高度,以10米每秒的速度做平抛运动 m=
1
;k=0.2;g=10;%参数 ts=0:0.01:2;%时间区间 初始值 x(0)=0 y(0)=10...
急!!!帮忙用
matlab
解一个
微分方程
(80分)
答:
clc;clear p=dsolve('Dp=lumda*(p1-p)-K*p-a*p')C1=solve(p-'p','C1')pt0='p1*((lumda/(K+lumda))-(lumda/(K+lumda))*p0/p1-(K/(K+lumda))*p0/p1)*exp(-(K+lumda)*t0)'C1=subs(C1,'t','t0')C1=subs(C1,'p',pt0)p=subs(p,'C1',C1)p=simplify(p)结果...
编程高手进
matlab
求解
微分方程
组
答:
函数中必须显式给出 dy(
1
)=...dy(2)=...dy(3)=...———函数没问题了,剩下的问题是 [t,y]=ode45(@weifensan,[0,3000],[10000;10000;0]);plot3(y(:,1),y(:,2),y(:,3))分号和括号要在英文半角状态下输入,而你用的是中文 ———hold on ...
棣栭〉
<涓婁竴椤
3
4
5
6
8
7
9
10
11
12
涓嬩竴椤
灏鹃〉
其他人还搜