matlab用od45解微分方程组。求改错。

2013-06-10 51自学网

a、b、c的取值有问题，长度不为10，检查一下。

法一：在script和function里面申明abc为全局变量

global a b c

法二：使用匿名函数

clc
format long

t=0.001:0.001:20;
x0=[3 -4 2];
a=0.3:0.3:3;
b=-6:-1:-15;
c=-1:-1:-10;

for k=1:10
    f=@(t,x)([a(k)*x(1)-x(2)*x(3);b(k)*x(2)+x(1)*x(3);c(k)*x(3)+x(1)*x(2)]);
    [t,x]=ode45(f,t,x0);
    plot(t,x);
    hold on
end
xlabel('t'),legend('x','y','z');

为什么我的matlab不能自动调用?

matlab 如何把散点数据通过插值画出3维的图形

CAD教程网

51自学网

教师招聘

公招网

医院招聘网