资源描述：
《MATLAB 选修课作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、课程：计算方法学院：信息工程学院班级：电信0906姓名：黄星星学号：0120909310615第一个图clearcloseallb=1000;a=[11000];w=[0:1000*2*pi];[hf,w]=freqs(b,a,w);subplot(1,3,1)plot(w/2/pi,abs(hf));grid;xlabel('f/Hz');ylabel('幅度');fs0=[1000,500];form=1:2fs=fs0(m);[d,c]=impinvar(b,a,fs);wd=[0:512]*pi/512;hw1=freqz(d,c,wd);subplot(1,3,2);plot(wd/
2、pi,abs(hw1)/abs(hw1(1)));holdon;endgrid;xlabel('f/(Hz)');text(0.52,0.88,'t=0.002s');text(0.12,0.54,'T=0.001s');form=1:2fs=fs0(m);[f,e]=bilinear(b,a,fs);wd=[0:512]*pi/512;hw2=freqz(f,e,wd);subplot(1,3,3)plot(wd/pi,abs(hw2)/abs(hw2(1)));holdon;endgrid;xlabel('f/(Hz)');text(0.5,0.74,'T=0.002S');text(0.
3、12,0.34,'t=0.001s');得到波形如图（1）Figure1图（1）Figure1第二个图clcclearcloseallN=21;wc=pi/4;n=0:N-1;r=(N-1)/2;hdn=sin(wc*(n-r))/pi./(n-r);ifrem(N,2)~=0hdn(r+1)=wc/pi;endwn1=boxcar(N);hn1=hdn.*wn1';wn2=hamming(N)hn2=hdn.*wn2';subplot(2,2,1)stem(n,hn1,'.')line([020],[00]);title('矩形窗设计的h(n)');xlabel('n');ylabel('h
4、(n)');subplot(2,2,3)stem(n,hn2,'.')line([0,20],[00])title('hamming窗设计的h(n)');xlabel('n');ylabel('h(n)');hn11=fft(hn1,512);w=2*[0:511]/512;subplot(2,2,2)plot(w,20*log10(abs(hn11)))gridaxis([02-805]);title('幅度特性');xlabel('w/pi');ylabel('幅度（dB)');hn22=fft(hn2,512);subplot(2,2,4)plot(w,20*log10(abs(hn22
5、)))gridaxis([02-805]);title('幅度特性');xlabel('w/pi');ylabel('幅度（dB)');得到波形如图（2）Figure2图（2）Figure2第三个图clcclear[t,omg,ft,ift]=prefourier([0,5],1000,[-250,250],1000);g=3*cos(10*t)+2*cos(20*t);%由定义生成基带信号f=g.*cos(100*t);%调制g0=f.*cos(100*t);%解调G0=ft*g0;%解调输出的频谱h=(omg<30&omg>-30);%定义低通滤波器G1=G0.*h;%在频域进行低通滤波g
6、1=ift*G1;%逆变换得到时域输出G=ft*g;F=ft*f;subplot(4,1,1)plot(t,g)xlabel('t')ylabel('g(t)')subplot(4,1,2)plot(t,f)xlabel('t')ylabel('f(t)')subplot(4,1,3)plot(t,g0)xlabel('t')ylabel('g0(t)')subplot(4,1,4)plot(t,g1)xlabel('t')ylabel('g1(t)')figuresubplot(5,1,1)plot(omg,G)axis([-250250-28])xlabel('w')ylabel('G(w
7、)')subplot(5,1,2)plot(omg,F)axis([-250250-24])xlabel('w')ylabel('F(w)')subplot(5,1,3)plot(omg,G0)axis([-250250-24])xlabel('w')ylabel('G0(w)')subplot(5,1,4)plot(omg,h)axis([-250250-0.51.5])xlabel('w')y
