如何实现MATLAB中的循环运算

yoyoo

我碰到一个比较棘手的问题,希望大家帮忙解答一下.
自己编写的程序,运行起来就有问题;还有,关于函数调用上不太明白,希望大家指点下.

自己编写的程序如下:
%
function f=yo_11(a,b)
l=0.5;m=0.1;mc=1;
f=(9.8*sin(a)-m*l*b^2*cos(a)*sin(a)/(mc+m))/l/(4/3-m*cos(a)^2/(mc+m));
%
function g=yo_12(a)
l=0.5;m=0.1;mc=1;
g=cos(a)/(mc+m)/l/(4/3-m*cos(a)^2/(mc+m));
%
function [A;A1]=yo_02()
h=0.01;A[1 1]=0.3491;A1[1 1]=0;
for(i=1:1000)
   x=A[1 i]+A1[1 i]*h;
   u=-25*A[1 i]-A1[1 i]+0.2*cos((2*pi/50)*i);
   f=yo_11(A[1,i],A1[1,1]);
   g=yo_12(A[1,i]);
   y=A1[1 i]+(f+g*u)*h;
   A=[A [x]];
   A1=[A1 [y]];
end
%
function C=yo_01()
k=0.95;D[1,i]=[A[1,i];A1[1,i];u];
F=[1,0,0;0,1,0;0,0,1];
C=[0,0;0,0;0,0];
for(i=1:1000)
   e=[A[1,i];A1[1,i]]-C`*D[1,i]`;
   F=1/k*(F-F*D[1,i]`*D[1,i]*F/(k+D[1,i]*F*D[1,i]`));
   C=C+F*D[1,i]`*e;
end
%
怎么才能使我编写的函数能运行起来,然后得到矩阵C.
麻烦大家了

eight

原帖由 yoyoo 于 2007-4-10 11:47 发表
我碰到一个比较棘手的问题,希望大家帮忙解答一下.
自己编写的程序,运行起来就有问题;还有,关于函数调用上不太明白,希望大家指点下.

自己编写的程序如下:
%
function f=yo_11(a,b)
l=0.5;m=0.1;mc=1;
f= ...

请阅读matlab的基础书，看看矩阵的定义和调用方式等