q609.m

来自「dsp相关的matlab程序」· M 代码 · 共 18 行

%《MATLAB在电子信息课程中的应用》第五章例5.9程序q509
%  非正弦信号的频谱分析和滤波
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear, format compact
L=0.4;C=10e-6;R=200;					% 输入元件参数
Um=100;w1=100*pi;
% 输入只有偶次谐波，故只分析偶次谐
N=input('需分析的谐波次数2N=（键入偶数）  ');
n=1:N/2;w=[eps,2*n*w1];					% 设定频率数组
Us = 4*Um/pi*[0.5,-1./(4*n.^2-1)];	% 输入信号频谱数组
z1=j*w*L;z2=1./(j*w*C);z3=R;		% 计算阻抗数组
z23=z2.*z3./(z2+z3);			% 注意为何用元素群运算符
Ur=Us.*z23./(z1+z23)			% 求Ur
disp('  谐波次数   谐波幅度  谐波相移(度)')
disp([2*[0,n]',abs(Ur)',angle(Ur)'*180/pi])