matlab好心人帮帮忙,我弄了好几个星期了,
来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/09/29 17:28:14
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--')
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--')
你这样写错的很离谱,而且我看不出来你到底想求什么
你说一下你要求解的微分方程和初值吧,不要用程序说,最好用公式编辑器
你说一下你要求解的微分方程和初值吧,不要用程序说,最好用公式编辑器
这个问题前几天我在百度知道里已经说过了,也有好几个好心人给了我答案,但我还是不太放心,所以今天又来到这里再补充一下.
该选文还是理?我怕报了理科后会跟不上,但是报文科的话又会后悔,马上就分科了,因为成绩不好烦了好几个星期,我想要更多的意见
求一些周记话题写了好几个星期的周记了,感觉没什么写的了,能提供些话题么?现在我十分郁闷啊!英语语文都可以·
谁知道我的 fine with me 写了好几个月了 别人总误解.弄的不好意思的!
我长大了 作文.各位好心人帮帮忙. 要求:童趣、好词好句多、600字. 在线等,快点啊!
中文摘要翻译成英文,哪位好心人帮帮忙,我的英语水平太差了
暑假我学会了什么的作文急需!好心人帮帮忙400-500字的
英语题目好心人帮帮忙教下我吧~谢谢了!在线等哦 !
我要关于“读书有益”的名人名言,请各位好心人帮帮忙,多谢了
想念老师的作文我们班老师生病了,好几个星期都没回来,别的班老师帮我班上的课.要写出对老师的思念
如何用skype加上外国人聊天,为什么我加了好几个外国人都没理我
问个数学题 好心人帮帮忙