声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1732|回复: 4

求助,用ODE45来求解连续方程,但是程序运行出现错误

[复制链接]
发表于 2007-5-6 09:38 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我用ODE45来求解连续方程,但是程序运行出现错误,希望高手帮我查找一下!!
仿真主程序:clear all;
close all;
global F
x0=[0,0,0,0];
ts=0.001;
T=5.0;
TimeSet=[0:ts:T];
para=[];
[t,y]=ode45('chap_4eq',TimeSet,x0,[],para);
switch F
    case 1
        r1=pi/6;r2=pi/6;
        dr1=0;dr2=0;
    case 2
        r1=sin(3*t);r2=cos(3*t);
        dr1=3*cos(3*t);dr2=-3*sin(3*t);
    case 3
        r1=sign(sin(3*t));r2=sign(sin(3*t));
        dr1=0;dr2=0;
end
figure(1);
plot(t,r1,'r',t,y(:,1),'b');
xlable('time(s)');ylable('r1,x(1)');
figure(2);
plot(t,dr1,'r',t,y(:,2),'b');
xlable('time(s)');ylable('dr1,x(2)');
figure(3);
plot(t,r2,'r',t,y(:,3),'b');
xlable('time(s)');ylable('r3,x(3)');
figure(4);
plot(t,dr2,'r',t,y(:,4),'b');
xlable('time(s)');ylable('dr2,x(2)');
控制子程序:
function dx=PlantModel(t,x,flag,para)
global F
dx=zeros(4,1);
F=2;
switch F
    case 1
        r1=pi/6;r2=pi/6;
        dr1=0;dr2=0;
        ddr1=0;ddr2=0;
    case 2
        r1=sin(3*t);r2=cos(3*t);
        dr1=3*cos(3*t);dr2=3*sin(3*t);
        ddr1=-3*3*sin(3*t);ddr2=-3*3*cos(3*t);
    case 3
        r1=sign(sin(3*t));r2=sign(sin(3*t));
        dr1=0;dr2=0;
        ddr1=0;ddr2=0;
end
dr=[dr1;dr2];
ddr=[ddr1,ddr2];
e1=x(1)-r1;
de1=x(2)-dr1;
e2=x(3)-r2;
de2=x(4)-dr2;
kp1=3000;kv1=230;
kp2=2000;kv2=210;
D11=2.462;
D22=0.362;
D12=0.147;
D21=0.147;
M=[D11 D12*cos(x(1)-x(3));
   D21*cos(x(1)-x(3)) D22];
G=[D12*x(4)^2*sin(x(1)-x(3));
   -D12*x(2)^2*sin(x(1)-x(3))];
C=[-x(4)*sin(x(3)) -(x(2)+x(4))*sin(x(3));
    x(2)*sin(x(3)) 0];
S=3;
switch S
    case 1
        u1=-kp1*e1-kv1*de1;
        u2=-kp2*e2-kv2*de2;
    case 2
        u1=-kp1*e1-kv1*de1+M(1,:)*ddr+C(1,:)*dr+G(1);
        u2=-kp2*e2-kv2*de2+M(2,:)*ddr+C(2,:)*dr+G(2);
    case 3
        Md=[D11 D12*cos(r1-r2);
            D21*cos(r1-r2) D22];
        Gd=[D12*dr1^2*sin(r1-r2);
            -D12*dr2^2*sin(r1-r2)];
        Cd=[-dr2*sin(r2) -(dr1+dr2)*sin(r2);
            dr1*sin(r2)  0];
        u1=-kp1*e1-kv1*de1+Md(1,:)*ddr+Cd(1,:)*dr+Gd(1);
        u2=-kp2*e2-kv2*de2+Md(2,:)*ddr+Cd(2,:)*dr+Gd(2);
end
tol=[u1;u2];
Q=inv(M)*(tol-C*[x(2);x(4)]-G);
dx(1)=x(2);
dx(2)=Q(1);
dx(3)=x(4);
dx(4)=Q(2);

[ 本帖最后由 xinyuxf 于 2007-5-8 15:28 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-5-8 11:40 | 显示全部楼层
这么长的程序:)请说下是哪里错了,否则太长了...
 楼主| 发表于 2007-5-8 18:18 | 显示全部楼层

回复

后面的子程序里运行时说需要定义T和X,但是X是一个变量,不知道如何设置!!有知道的说下,谢谢!!
发表于 2007-5-10 19:19 | 显示全部楼层

这个问题现在解决了?
怎么解决的?
 楼主| 发表于 2007-5-11 09:52 | 显示全部楼层

回复 #4 cdwxg 的帖子

问题解决了!!我也没改它,在同学的电脑上运行不了,在自己的电脑上随便改了下就可以运行了!!我再看看有什么不同的地方吧!!!:@)
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 16:41 , Processed in 0.061103 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表