%《MATLAB在電子信息課程中的應(yīng)用》第六章例6.10程序q610
% 調(diào)幅信號通過窄帶濾波器
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear
N=1000;
t=linspace(0,2*pi,N);dt=2*pi/(N+1);		% 信號周期為2*pi,分成1000份
w=[99,100,101];			% 輸入信號的三個頻率分量
U=[0.5,1,0.5];			% 三個頻率分量對應(yīng)的向量（虛部為零）
b=[2,0];a=[1,2,10001];	% 濾波器分子分母系數(shù)向量
u1=U*cos(w'*t+angle(U')*ones(1,N));	% 輸入信號的時間曲線
H=polyval(b,j*w)./polyval(a,j*w);	% 求濾波器在三個頻點上的頻率響應(yīng)，也可用H=freqs(b,a,w);
% 畫出濾波器的頻率響應(yīng)曲線，只用三個頻點，圖形不好看，讀者可修改程序得到其完整的頻率響應(yīng)
figure(1)
subplot(2,1,1),plot(w,abs(H)),grid			% 幅度
subplot(2,1,2),plot(w,angle(H)),grid		% 相位
%U2=U.*H;
%ewn=exp(w'*t*i);u2=U2*ewn;  % 求u2=idtft(U2,w,t);
u21=abs(U(1)*H(1))*cos(99*t+angle(U(1)*H(1)));	% 角頻率為99的分量	
u22=abs(U(2)*H(2))*cos(100*t+angle(U(2)*H(2)));	% 角頻率為100的分量
u23=abs(U(3)*H(3))*cos(101*t+angle(U(3)*H(3)));	% 角頻率為101的分量
u2=u21+u22+u23;								% 求和
% 巧妙地利用元素群運算和矩陣運算相結(jié)合可把四條語句合成一條語句如下
% u2=abs(U.*H)*cos(w'*t+angle((U.*H).')*ones(1,1001));
% 注意對復(fù)數(shù)矩陣(U.*H)，(U.*H)'為其共軛轉(zhuǎn)置，(U.*H).'為轉(zhuǎn)置而不共軛
figure(2)							% 畫出原信號和濾波后信號的波形作比較
subplot(2,1,1),plot(t,u1)
subplot(2,1,2),plot(t,u2)