matlab好心人帮帮忙,我弄了好几个星期了,function dydt=newday(t,y)global k r10 rou3 rou1 miu3 miu1 n w m s gama p0 pv pa f sig1 sig2;dydt=[y(2);((-y(2)^2)*0.5*rou3*y(1)*(3*(y(1)^3)+4*k)*(1/((y(1)^3)+k)^s)-1.5*rou1*y(1)^2-4*miu3*y(2)*(y(1)
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 17:15:08
![matlab好心人帮帮忙,我弄了好几个星期了,function dydt=newday(t,y)global k r10 rou3 rou1 miu3 miu1 n w m s gama p0 pv pa f sig1 sig2;dydt=[y(2);((-y(2)^2)*0.5*rou3*y(1)*(3*(y(1)^3)+4*k)*(1/((y(1)^3)+k)^s)-1.5*rou1*y(1)^2-4*miu3*y(2)*(y(1)](/uploads/image/z/3163756-4-6.jpg?t=matlab%E5%A5%BD%E5%BF%83%E4%BA%BA%E5%B8%AE%E5%B8%AE%E5%BF%99%2C%E6%88%91%E5%BC%84%E4%BA%86%E5%A5%BD%E5%87%A0%E4%B8%AA%E6%98%9F%E6%9C%9F%E4%BA%86%2Cfunction+dydt%3Dnewday%28t%2Cy%29global+k+r10+rou3+rou1+miu3+miu1+n+w+m+s+gama+p0+pv+pa+f+sig1+sig2%3Bdydt%3D%5By%282%29%3B%28%28-y%282%29%5E2%29%2A0.5%2Arou3%2Ay%281%29%2A%283%2A%28y%281%29%5E3%29%2B4%2Ak%29%2A%281%2F%28%28y%281%29%5E3%29%2Bk%29%5Es%29-1.5%2Arou1%2Ay%281%29%5E2-4%2Amiu3%2Ay%282%29%2A%28y%281%29)
matlab好心人帮帮忙,我弄了好几个星期了,function dydt=newday(t,y)global k r10 rou3 rou1 miu3 miu1 n w m s gama p0 pv pa f sig1 sig2;dydt=[y(2);((-y(2)^2)*0.5*rou3*y(1)*(3*(y(1)^3)+4*k)*(1/((y(1)^3)+k)^s)-1.5*rou1*y(1)^2-4*miu3*y(2)*(y(1)
matlab好心人帮帮忙,我弄了好几个星期了,
function dydt=newday(t,y)
global k r10 rou3 rou1 miu3 miu1 n w m s gama p0 pv pa f sig1 sig2;
dydt=[y(2);((-y(2)^2)*0.5*rou3*y(1)*(3*(y(1)^3)+4*k)*(1/((y(1)^3)+k)^s)-1.5*rou1*y(1)^2-4*miu3*y(2)*(y(1)^2)*(1/(k+y(1)^3))-4*miu1*y(2)*(1/y(1))-p0+pv-pa*sin(2*pi*f*t)+pg0*(r10^w)*(1/y(1)^w)-gama*pg0*(r10^w)*(1/y(1)^w)+gama*pg0*(r10^n)*(1/y(1)^n)-2*sig1*(1/y(1))-2*sig2*(1/y(2)))*((k+y(1)^3)^m)*(1/(rou3*y(1)^2-rou1*y(1)*(k+y(1)^3)^m))];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
global k=4.0703e-19;r10=3e-6;rou3=-100;rou1=1100;miu3=-0.279;miu1=0.28;n=6.21;w=3.21;m=1/3;s=4/3;gama=1.07;p0=1.013e5;pv=2.33e3;pa=1.7e5;f=2.62e4;sig1=0.04;sig2=0.073;
[t,y]=ode45(@newday,[0 40],[]);
plot(t,y(:,1),'r--')
matlab好心人帮帮忙,我弄了好几个星期了,function dydt=newday(t,y)global k r10 rou3 rou1 miu3 miu1 n w m s gama p0 pv pa f sig1 sig2;dydt=[y(2);((-y(2)^2)*0.5*rou3*y(1)*(3*(y(1)^3)+4*k)*(1/((y(1)^3)+k)^s)-1.5*rou1*y(1)^2-4*miu3*y(2)*(y(1)
你这样写错的很离谱,而且我看不出来你到底想求什么
你说一下你要求解的微分方程和初值吧,不要用程序说,最好用公式编辑器