10.50-0.5-100.020.040.060.080.10.120.140.160.180.2200015001000500005001000150020002500300035004000
后来我在网上搜了一下关于正弦信号的频谱分析问题： 得到下面的结果： 程序：
t=-1:0.01:1; x=5*sin(2*pi*10*t); N=length(x); fx=fft(x); df=100/N; n=0:N/2; f=n*df; subplot(211); plot(t,x);grid on; subplot(212);
plot(f,abs(fx(n+1))*2/N); grid;
仿真结果：
50-5-15432100-0.8-0.6-0.4-0.200.20.40.60.815101520253035404550
这是对f=1/T=w/2?而言的，但是还是有点不懂。因为我们要求的是关于w而言的，所以我把上面的程序稍微改了一下，得到如下程序：
t=-1:0.01:1; x=5*sin(2*pi*10*t); N=length(x); fx=fft(x); df=100/N; n=0:N/2; w=n*df*2*pi; subplot(211); plot(t,x);grid on; subplot(212);
plot(w,abs(fx(n+1))*2/N); grid;
仿真结果为：
50-5-15432100-0.8-0.6-0.4-0.200.20.40.60.8150100150200250300350
这和我预想的结果不一样，我觉得应该在w=20的地方有冲击呀，但结果却差了好远。
【研讨内容】 题目2．分析实际物理系统的频率响应。
R25R21200kgnd2R22180k200kVCC8R26180kVCC8gnd6U1AU1BLF353UiViSINC15470pC163470pR19200kLF3531TP13C17C185470pR23200k7TP14470pR20130k4R24100kVEEgndgndgnd4VEE
【数学模型】 接入端。输出端。
通过分析：发现（1）和（2）的系统描述还有点相似，所以把它弄他们弄到一块讨论，（3）和（4）同点是（3）是对称的，而（4）中x[k-2]没有，且两边的数是相反数，所以也把他们弄到一块讨论。
此题考查的是离散系统的频率响应的求解，可根据描述方程的下列频率响应：
（1）H1(e(2) H1(e(3) H1(e(4) H1(ej?)?1.69?1.05e?j??0.45e?2j??0.45e?3j??1.05e?4j??1.69e?5j?; )?1.69?1.05e?j??0.45e?2j??0.45e?3j??1.05e?4j??1.69e?5j?; )?1.36?0.61e?j??1.08e?2j??0.61e?3j??1.36e?4j?; )?1.36?0.61e?j??0.61e?3j??1.36e?4j?;
j?j?j?再照着课本求离散系统频率响应。
【仿真程序】 （1） 和（2）
a1=[1]; a2=[1];
b1=[1.69 1.05 0.45 0.45 1.05 1.69]; b2=[1.69 1.05 0.45 -0.45 -1.05 -1.69]; [h1,w]=freqz(b1,a1,200); [h2,w]=freqz(b2,a2,200); h11=abs(h1);
h12=angle(h1); h21=abs(h2); h22=angle(h2); subplot(2,2,1); plot(w,h11); grid;
xlabel('角频率(w)'); ylabel('幅度'); title('(1)的相频特性'); subplot(2,2,2); plot(w,h12*180/pi); grid;
xlabel('角频率(w)'); ylabel('相位'); title('(1)的相频特性'); subplot(2,2,3); plot(w,h21); grid;
xlabel('角频率(w)'); ylabel('幅度'); title('(2)的相频特性'); subplot(2,2,4); plot(w,h22*180/pi); grid;
xlabel('角频率(w)'); ylabel('相位'); title('(2)相频特性');
（3）和（4） a3=[1]; a4=[1];
b3=[1.36 0.61 1.08 0.61 1.36]; b4=[1.36 0.61 0 -0.61 -1.36]; [h3,w]=freqz(b3,a3,200); [h4,w]=freqz(b4,a4,200); h31=abs(h3); h32=angle(h3); h41=abs(h4); h42=angle(h4); subplot(2,2,1); plot(w,h31); grid;
xlabel('角频率(w)'); ylabel('幅度');