求系统幅频响应和相频响应3.txt

求系统幅频响应和相频响应 求系统零极点图 二阶系统级联式

>> num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528];
>> den=[1 -1.8707 2.4947 -1.8801 0.9537 -0.2336];
>> freqz(num,den);   %   0~π中抽样，抽样点缺省（512点） ζ


>> num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528];
>> den=[1 -1.8707 2.4947 -1.8801 0.9537 -0.2336];
>> w=[0 pi/8 pi/4 pi*3/8 pi/2 pi*5/8 pi*3/4];  %自己定8个点  Θ
>> h=freqz(num,den,w);
>> subplot(2,2,1);
>> stem(w/pi,abs(h));
>> title('幅度谱')
>> xlabel('数字频率');
>> ylabel('振幅');

>> [h,w]=freqz(num,den,8);          % 系统在0~π之间均分8份，与“Θ”处效果一样
>> w
>> subplot(2,2,2);
>> stem(w/pi,abs(h));
>> title('幅度谱')
>> xlabel('数字频率');
>> ylabel('振幅'); 

>> h=freqz(num,den);          % 系统在0~π之间均分512份，与“ζ”处效果一样
>> subplot(2,2,3);
>> z=10*log(abs(h))
>> plot(z);                 % 与“ζ”处幅度谱效果一样
>> title('分贝幅度谱')
>> xlabel('数字频率');
>> ylabel('振幅'); 


>> num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528];
>> den=[1 -1.8707 2.4947 -1.8801 0.9537 -0.2336];
>> [z,p,k]= tf2zp(num,den);   % 求零极点
>> z  %零点
>> p  %极点
>> subplot(2,2,1);
>> zplane(z,p);   %  zplane (num,den)也可以
>> [r,p,k]=residuez(num,den);  % 部分分式展开
>> r
>> p
>> k
>> [sos,g]=zp2sos(z,p,k);   %二阶系统分解
>> sos
>> g
--------------------------------------------------------------------------------------