[MATLAB毕业设计] 光耦合方程用ode45编程如下，请高手看看那里错了？ - MATLAB爱好者论坛-LabFans.com

wumaobo · 2011-04-19, 21:10

目标函数
function dy=rigid(t,y)
global g A as vp vs ap;
dy=zeros(2,1);
dy(1)=g*y(1)*y(2)/A-as*y(1);
dy(2)=-vp*g*y(1)*y(2)/(vs*A)-ap*y(2);
主函数
global g A as vp vs ap ;
g=0.609*10^(-13);A=8*10^(-13);as=-25;ap=4*10^8;vp=1550;vs=1450;
[T,Y]=ode45(@rigid,[0 40],-25);
plot(T,Y)
请问那里出错

2011-04-19, 21:10	#1
wumaobo 初级会员注册日期: 2011-04-07 帖子: 2 声望力: 0	光耦合方程用ode45编程如下，请高手看看那里错了？目标函数 function dy=rigid(t,y) global g A as vp vs ap; dy=zeros(2,1); dy(1)=gy(1)y(2)/A-asy(1); dy(2)=-vpgy(1)y(2)/(vsA)-apy(2); 主函数 global g A as vp vs ap ; g=0.60910^(-13);A=810^(-13);as=-25;ap=4*10^8;vp=1550;vs=1450; [T,Y]=ode45(@rigid,[0 40],-25); plot(T,Y) 请问那里出错

主题工具
显示可打印版本邮寄本页给好友
显示模式
切换到平板模式切换到混合模式树形模式