




已阅读5页,还剩2页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
MATLAB机考样题:(1)Generate and plot sequence , with.n1=-20:20;x1=2*cos(pi/8*n1);n2=n1-4;x2=2*cos(pi/8*n2);subplot(2,1,1);plot(n1,x1);subplot(212);plot(n2,x2);(2)Write a MATLAB program to compute and plot the impulse response of a causal finite-dimensional discrete-time system characterized by a difference equation of the following form: N=input(请输入你要求的点数N=);num=1.8 0.34 -1.32 -0.86;den=1 0.3 0.5 -0.72;x=1 zeros(1,N-1);(单位冲击)y=filter(num,den,x);plot(0:N-1,y);(3)Write a MATLAB program to compute and display the poles and zeros, to compute and display the second-order factored form, and to generate the pole-zero plot of a z-transform that is a ratio of two polynomials in z-1. Using this program to analyze the following G(z):num=8.1 6.93 -23.82 10.5;den=1 1.52 0.18 -0.1768;sos=tf2sos(num,den) %tf2sos表示为1/z的升幂zplane(num,den)(4)Try to give a program to evaluate the following DTFT in the range :%由于用freqz计算频点至少是2个,所以至少输入两个频点w1=input(请输入你要计算的频点w1=);w2=input(请输入你要计算的频点w2=);w=w1 w2;num=2 5 9 5 3;den=5 45 2 1 1;h=freqz(num,den,w)(6)Write a MATLAB program to compute and plot the magnitude response of a causal LTI discrete-time system with a transfer function given bynum=0.15*1 0 -1;den=1 -0.5 0.7;h,w=freqz(num,den,512);plot(w/pi,abs(h);(7)Consider the following FIR transfer function:Using MATLAB to determine its zero locations and plot its magnitude and phase response.h=1 0.6 .49 -0.48 -0.14 -0.12 0.09;figure(1)zplane(h,1);H,w=freqz(h,1,512);figure(2)plot(w/pi,abs(H);figure(3)plot(w/pi,angle(H);(8)Given a signal , when using a sampling frequency fT= 20KHz, plot the magnitude and phase spectrum of the sampled sequence(given length-64).fs=2e4;n=0:63;x=4*n/fs+cos(0.1*pi*n/fs);h=fft(x,1024);figure(1);plot(0:2/1023:2,abs(h);figure(2);plot(0:2/1023:2,angle(h);(9)design an IIR butterworth digital lowpass filter with the following specifications: sampling rate of 40kHz, passband edge frequency of 4kHz, stopband edge frequency of 8kHz, passband ripple of 0.5dB, and a minimum stopband attenuation of 40dB,plot frequency-magnitude and check if your design fits the specification.fs=40;wp=4*2/fs; %wp1,没有乘piws=8*2/fs; %ws1,没有乘piap=0.5;as=40;n,wn=buttord(wp,ws,ap,as);num,den=butter(n,wn);h,w=freqz(num,den,512);figure(1);plot(w/pi,20*log10(abs(h);axis(0 1 -50 0);figure(2);subplot(2,1,1);plot(w/pi,20*log10(abs(h)axis(0 wp -0.5 0);title(通带纹波);subplot(2,1,2);plot(w/pi,20*log10(abs(h);axis(ws 1 -50 -30);title(阻带纹波); (10)Design a Hanning FIR lowpass filter meeting the following specifications: passband edge frequency=2kHz, stopband edge frequency=2.5kHz, passband ripple p=0.005, stopband ripples=0.005, and sampling rate of 10kHz.Plot its gain and phase responses and check if it meets the specifications?ft=10;fp=2;fs=2.5;wp=2*pi*fp/ft;ws=2*pi*fs/ft;ds=0.005;ap=20*log10(1-ds)as=20*log10(ds)wc=(wp+ws)/2;dw=ws-wp;M=ceil(3.11*pi/dw);N=2*M;b=fir1(N,wc/pi,hann(N+1);h,w=freqz(b,1,512);figure(1);plot(w/pi,20*log10(abs(h);axis(0 1 -50 0);title(magitude response);figure(2);plot(w/pi,unwrap(angle(h);title(phase response);figure(3);subplot(211);plot(w/pi,20*log10(abs(h);axis(0 wp/pi ap 0);title(通带纹波)subplot(212);plot(w/pi,20*log10(abs(h);axis(ws/pi 1 as 0);title(阻带纹波);%从图中可以看出,通带和阻带中纹波都不满足要求,所以不满足指标%as= -46.020643.9 所以不能用hanning窗设计%应当用hamming或blackman窗设计(11)Writing a MATLAB program to compute 128-point DFT of the following sequence, you must firstly use DFT definition (directly computing DFT) to compute and use MATLAB function to test the result. Plot the two results in one figure.n=0:31;k=1:128;x=sin(pi*n/4);Xk1=zeros(1,128); for t=1:128 for m=1:32 Xk1(t)=Xk1(t)+x(m)*exp(-1i*2*pi*(t-1).*(m-1)/128); endendsubplot(2,1,1);plot(k,Xk1);Xk2=fft(x,128);subplot(2,1,2);plot(k,Xk2);(12) Using the function fir1 and window of Kaiser, design a linear-phase FIR lowpass filter meeting the following specifications: passband edge frequency=2kHz, stopband edge frequency=2.5kHz, passband ripple p=0.005, stopband ripples=0.005, and sampling rate of 10kHz.Plot its gain and phase responses and check if it meets the specifications?Design a Type 1 Chebyshev IIR lowpass filter meeting the specifications as below: sampling rate of 12kHz, passband edge frequency of 2.1kHz, stopband edge frequency of 2.7kHz, passband ripple of 0.6dB, and a minimum stopband attenuation of 45dB. Write down the exact expression for the transfer function generated. Does your design meet the specifications? Fp=2100;Fs=2700;Ft=12000;Rp=0.6;Rs=45;Wp=2*Fp/Ft;Ws=2*Fs/Ft;N,Wn=cheb1ord(Wp,Ws,Rp,Rs)B,A=cheby1(N,Rp,Wn)h,w=freqz(B,A,512);figure(1);plot(w/pi,20*log10(abs(h);axis(0 1 -50 0);figure(2);subplot(2,1,1);plot(w/pi,20*log10(abs(h)axis(0 Wp -0.6 0);title(通带纹波);subplot(2,1,2);plot(w/pi,20*log10(abs(h);axis(Ws 1 -50 -30);title(阻带纹波);(13)Using the function fir1 and window of Kaiser, design a linear-phase FIR lowpass filter meeting the following specifications: passband edge frequency=2kHz, stopband edge frequency=2.5kHz, passband ripple p=0.005, stopband ripples=0.005, and sampling rate of 10kHz.Plot its gain and phase responses and check if it meets the specifications?n,wn,beta,typ=kaiserord(2000 2500,1 0,0.005 0.005,10000);b=fir1(n,wn,kaiser(n+1,beta),noscale);h,omega=freqz(b,1,256);subplot(2,1,1);plot(omega/pi,20*log10(abs(h);xlabel(omega/pi); ylabel(Gain, dB);subplot(2,1,2)plot(omega/pi,angle(h);gridtitle(Phase Spectrum);xlabel(omega/pi); ylabel(Phase, radians);(14)Given a signal when using a sampling frequency fT= 20KHz, plot the magnitude and phase spectrum of the sampled sequence(given length-64). ft=20000;n=1:64;x=sin(0.1*pi*n/ft)+2*cos(0.3*pi*5*n/ft)+3*sin(0.5*pi*n/ft);subplot(2,1,1);h,w=freqz(x,1,256);plot(w/pi,abs(h);title(Magnitude spectrum of the sampled samples);xlabel(omega/pi); ylabel(Magnitude);subplot(2,1,2);plot(w/pi,angle(h);title(Phase spectrum of the sampled samples);xlabel(omega/pi); ylabel(Phase);(15) Write a MATLAB program to compute the first L samples of the inverse of rational z-transforms where the value of L is provided by the user through the command input. Using this program to compute and plot the first 50 samples of the inverse of following G(z). Use the command stem for plotting the sequence gener
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 成长记录袋小学生课件
- 2025年环保厕所项目合作计划书
- 2025年重铬酸钠项目建议书
- 加强网络信息安全保障条例
- 公司股份制实施方案
- 金融投资顾问投资风险提示书
- 小王子电影故事解读
- StA-IFN-1-生命科学试剂-MCE
- 石油库区员工年终总结
- 2025年太阳能热发电系统项目合作计划书
- 图书外借服务计划
- 软考系统集成项目管理工程师教程完整版
- GB/T 45091-2024塑料再生塑料限用物质限量要求
- 人教版八年级上册地理 2024-2025学年八年级上册地理期中测试卷(二)(含答案)
- 危险性较大的分部分项工程清单和安全管理措施范文
- 2024-2025年江苏专转本英语历年真题(含答案)
- 投标废标培训
- 脑卒中课件完整版本
- 药房保洁流程规范
- 装修合同违约解除通知书
- (新版)六西格玛绿带认证考试复习题库(含答案)
评论
0/150
提交评论