[紧急求助]一个Matlab小程序

1.pdf (219.53 KB)

2010-4-24 17:07:49 上传

学习中遇到了一个小问题，需要编写一个Matlab小程序，由于本人不懂 Matlab，因此急切希望某位大侠能相助！
价钱还好商量，请各位仁兄帮忙！

、    

%第一题的主程序，ti1st.m
[t,Y] = ode45(@kalmon1,[0 1],[0 0 -1])
%把自变量取区间[0,1]。delta_x的初值直接取成了0，不知道是否合理
plot(t,Y,'r--')

%第一题的函数文件 kalmon1.m
function dy = kalmon1(t,y)
%y(1)即delta_x，y（2）即P,y(3)即x
dy = zeros(3,1);    % a column vector
q=0.5;
r=0.6;
dy(1) = -y(1)+3*(5-6*exp(-t))^2*y(2)*(-3*(5-6*exp(-t))^2*y(1))/r; %在这行的-3*(5-6*exp(-t))^2*y(1)前面加上delta_z（t）的表达式即可
dy(2) = -2*y(2)-9*(5-6*exp(-t))^4*y(2)^2/r+q;
dy(3) = 5-6*exp(-t)+y(1);

这种程序有难度，已经不是小程序了。牵涉到非线性方程组还有Kalman Filter，如果有人来做，估计没有几天时间不能做出来。

不是呀，有解题的方法了，不难的吧

希望有人会做

先看看，matlab我也懂一点，看能出上力不……

%第一题的主程序，ti2nd.m

[t,Y] = ode45(@rigid,[0 2],[0 1])
plot(t,Y,'r--')

%第二题的函数文件，rigid.m
function dy = rigid(t,y)
%y(1)即x，y（2）即P
dy = [-1,0]';    % a column vector
q=0.5;
r=0.6;
dy(1) = -y(1)+3*y(1)^2*y(2)*(-y(1)^3)/r+5;%在这行的-y（1）^3前面加上Z（t）的表达式即可
dy(2) = q-2*y(2)-9*y(1)^4*y(2)^2/r+q;

%=======已上是全部程序，由于某些符号我还不大懂，所以尤其是Z（t）不知为何物，就直接忽略掉了。请按照程序的注释修改再使用。