画分岔图，出现Vectors must be the same lengths.

lxx244lxx

function tanqingweiqi
hold on;box on;%生成空的坐标轴
N=100;%设置迭代次数
Tn=[];%初始化
kg2a=5.83;
options=odeset('RelTol',1e-9);%设置微分方程求解选项
for kg1a=5:0.01:8;%计算不同的kg1a值
    x0=[0,0,500];%设置初值
    tspan=[0:0.001:3]
   y=ode4(@ODE317_2,tspan,x0,options,kg1a,kg2a)
    x0=y(end,:);
    for n=1:N;
        y=ode4(@ODE317_2,[0:0.001:0.005],x0,options,kg1a,kg2a);%求解微分方程
        x0=y(end,:);%更新初值
        xd=y(:,3);%提取求解结果
        [M,IK]=max(xd);%计算最大值
        if 1<IK & IK<length(xd);%判断是否为最大值内部点
            Tn=[Tn;M];%记录最大值
        end
    end
    plot(kg1a,Tn,'k.','markersize',1);%绘图
    Tn=[];%清空
   % pause(0.01);
end
end

function dy=ODE317_2(t,y,flag,kg1a,kg2a)
dy=zeros(3,1);
r=0.5;pb=1175;Fco=1400;Fco2=400;F=20000;
r1x=(kg1a*16.84*(1-y(1))*533529*exp(-76857/8.314/y(3))*(16.84*(1-y(1)))^0.674)/(kg1a*16.84*(1-y(1))+533529*exp(-76857/8.314/y(3))*(16.84*(1-y(1)))^0.674);
r2x=(kg2a*4.81*(1-y(2))*296452638*exp(-105858/8.314/y(3))*(4.81*(1-y(2)))^0.008)/(kg2a*4.81*(1-y(2))+296452638*exp(-105858/8.314/y(3))*(4.81*(1-y(2)))^0.008);
dy(1)=(pi*r*r*pb/Fco)*r1x;
dy(2)=(pi*r*r*pb/Fco2)*r2x;
dy(3)=(pi*r*r*pb/F)*(r1x*1000*(188.1276+0.0724*y(3)-4.2782*10^(-5)*y(3).^2+7.9724*10^(-9)*y(3).^3)+r2x*1000*(146.4633+0.06996*y(3)-2.7313*10^(-5)*y(3).^2+1.8928*10^(-9)*y(3).^3))/(1-0.14*y(1)-0.04*y(2))/(26.2932+0.01354*y(3)+6.1860*10^(-6)*y(3).^2-4.1355*10^(-9)*y(3).^3);

ChaChing

Ref: 8F
[原创]写给学习 matlab 的新手们
http://forum.vibunion.com/forum/ ... p;extra=&page=1

lxx244lxx

回复 2 # ChaChing 的帖子

谢谢，我明白我的错误是什么意思，向量字节长度必须相同，但是看了好几遍感觉还是没问题啊，没找出错误在哪。。

gghhjj

lxx244lxx 发表于 2011-3-17 14:25

                               
登录/注册后可看大图

回复 2 # ChaChing 的帖子

谢谢，我明白我的错误是什么意思，向量字节长度必须相同，但是看了好几遍感觉还 ...

kg1a是一个数值
Tn=[Tn;M]是一个向量