|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我是新手,刚学matlab,想利用以下程序实现Y(:,1)的广义积分,下限是0,上限是无穷。望指教
clc,clear,close all;
global N
t0=0; tf=20; y1=0.66e-6; y2=0; a=1;
[T,Y]=ode45('li',[t0,tf],[y1,y2]);
b(a,2)=int(Y(:,1),0,inf); b(a,1)=N; a=a+1;
end
plot(b(:,1),b(:,2),'r-')
结果画的不对,应该是一条平滑的曲线,纵轴的范围应该是0到17左右。
li.m文件是
function dy=li(t,y)
global N
dy=zeros(2,1);
dy(1)=y(1)*N*exp(-y(2))-y(1);
dy(2)=y(1);
[ 本帖最后由 ChaChing 于 2010-1-14 10:09 编辑 ] |
|