Question

MATLAB的问题

1. .已知传递函数
G(s)=(6s^2+1)/(s^3+3s^2+3s+1)
H(s)=[(s+1)(s+2)]/[(s+2j)(s-2j)(s+3)]
用MATLAB求出G（S）的零极点、H（S）的多项式形式，及G（S）H（S）的零极点图。
2. 编写 s=1+4+9+16...+n^2的M函数文件，并计算当n=20 时s 的值。

答案对的我还可以加分

Answer 1

%第1题
%由于H＝（s^2+3*s+2)/(s^3+3*s^2+4*s+12)，可用下面式子：
num1=[6 0 1];den1=[1 3 3 1];
num2=[1 3 2];den2=[1 3 4 12];
G=tf(num1,den1);H=tf(num2,den2);
GH=feedback(G,H);
[z1,p1]=tf2zp(num1,den1) %求出G的零极点
[z3,p3]=tf2zp(GH.num{1},GH.den{1}) %求出G的零极点
pzmap(GH);%画出零极点图
grid on;
%第2题
n=20;
s=0;
for i=1:n;
s=s+i^2;
end
s

>>s =

2870
另外：
1+4+9+16...+n^2＝
n*(n+1)*(2*n+1)/6