求助：LabVIEW利用MATLAB script如何编写PID控制程序（P和PI的都行了）

是做毕业设计，以下是我在MATLAB script节点里编的P和PI的控制（都是对的），传递函数也在图中。由于PID的编不来，所以网上找了
Function PI（PI的）
G=tf(1,[1,3,3,1]);
Kp=1;Ti=[0.7:0.1:1.5];
for i=1:length(Ti)
Gc=tf(Kp*[1,1/Ti(i)],[1,0]); G_c=feedback(G*Gc,1);
step(G_c),hold on
end
axis([0,20,0,2])

Function PID（PID的）
N=[100,1000,10000,1:10];
G=tf(1,[1,3,3,1]);
Kp=1;Ti=1;Td=1;
Gc=tf(Kp*[Ti*Td,Ti,1]/Ti,[1,0]);
G_c=feedback(G*Gc,1);step(G_c), hold on
for i=1:length(N)
mn=Kp*([Ti*Td,0,0]+conv([Ti,1],[Td/N(i),1]))/Ti;
cd=[Td/N(i),1,0]; Gc=tf(mn,cd);
G_c=feedback(G*Gc,1);
step(G_c)
end
axis([0,20,0,2])
我按照着这些，然后类比下，自己尝试编了下，但是运行错误，大家能否帮我看看，或者会的，自己曾今做的提点我一下，帮我下。

注：在编PI时，K取的是1；在编PID时，K取1，Ti取1，N取10

希望大家帮帮我，挺急的，一直卡在这里

最后一幅图中，我发现我自己弄错了，把K改为Td后可以运行了，但是波形不对

wdyjz

可以查看下matlab的PID工具箱

朋友，这个程序怎么编啊,毕业论文比较重要
y=log(2+z(1)*x*y)+log(2+z(2)*x*y)+...+log(2+z(i)*x*y)+...+log(2+z(1000)*x*y),
其中，z(i)为已知的，x,y为变量，i从1到1000,
E-mail  yuxiyazhu@163.com