资源描述：
《现代测控分析软件与应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、成绩中国矿业大学2016级硕士研究生课程考试试卷考试科目：现代测控分析软件与应用考试时间：2016年5月13日学生姓名：******学号：**************所在院系：机电工程学院任课教师：***中国矿业大学研究生院培养管理处印制一，（1）Labview系统时间显示*****图1.程序前面板图2.程序框图2.具体程序：cur_time=fix(clock);str=sprintf('%d%s%d%s%d%s%d%s%d%s%.2d%s',cur_time(1),’年’,cur_time(2),’月’,cur_time(3),’日’,cur_tim
2、e(4),’时’,cur_time(5),’分’,cur_time(6),’秒’);msgbox(str);图3.matlab程序和显示二信号的产生和读取1.用附件中的录音机记录一段语音信号，命名为“wangyangyang.wav”的波形文件保存在D盘的input_data文件夹下，添加该路径，以便matlab识别此段语音。filename='wangyangyang.wav';[y,fs,nbits]=wavread(filename);%读取语音信号N=length(y);%求出语音信号长度time1=(1:N)/fs;subplot(2,1,1);p
3、lot(time1,y)%绘制原始语音信号图图4.原始信号图A1=y(4001:5000);%截取一段语音信号，幅值给A1N2=length(A1);%求出A1的信号长度subplot(2,1,2);time2=(4001:5000)/fs;plot(time2,A1)%绘制A1图形xlswrite('d:input_datawangyangyangtest.xlsx',A1,'A1:A1000');%保存成*.xlsx格式A1=A1';%将A1转置save('test.txt','A1','-ascii');%保存成-ascii格式图5.A1信号2.
4、图6.程序框图图7.前面板3.Fs=1000;%采样频率T=1/Fs;%采样时间L=1000;%信号长度t=(0:L-1)*T;%时间轴A2=0.3*sin(2*pi*17*t)+0.4*sin(2*pi*155*t)+0.5*sin(2*pi*430*t);A3=0.25+0.25*randn(size(A2));S=A2+A3%混合信号figureplot(t,S)%绘制混合信号title('A1和A2的混合信号');%标题xlabel('时间t');%横轴标签ylabel('幅值
5、S(f)
6、');%纵轴标签图8.混合信号figureplot(t,A2)
7、%绘制信号A2title('信号A2');%标题xlabel('时间t');%横轴标签ylabel('幅值
8、S(f)
9、');%纵轴标签图9.信号A2三信号处理1.解：已知系统的传递函数为：又，则因此，系统的差分方程为：num=[10.5];den=[10.25];sys=filt(num,den);impulse(sys)其脉冲响应如图：图10.脉冲响应2.matlab程序subplot(2,2,1);%绘制原始信号S波形图plot(t,S);title(’原始信号的波形图’);xlabel('时间t');ylabel('幅值
10、S(f)
11、');gridon%
12、设计IIR椭圆低通滤波器Wp=80/500;%通带截止频率Ws=110/500;%阻带截止频率Rp=1;%通带波纹度Rp=1dbRs=60;%阻带衰减率Rs=60db[n,wn]=ellipord(Wp,Ws,Rp,Rs);[a,b]=ellip(n,Rp,Rs,wn);subplot(2,2,2);[H,w]=freqz(a,b,215,Fs);%绘制滤波器的频率响应曲线Hf=abs(H);Hx=angle(H);plot(w,Hf,w,Hx,'-.');legend('幅频特性','相频特性');title('滤波器的相频响应曲线');xlabel('频
13、率（Hz）');ylabel('幅值/相位');gridonTS=filter(a,b,S);%对原始信号进行滤波，并绘制滤波器后信TS的波形图subplot(2,2,3);plot(t,TS)title('滤波后信号的波形图');xlabel('时间t');ylabel('
14、TS(f)
15、');gridonNFFT=2^nextpow2(L);%对滤波后信号进行快速傅里叶变换，并绘制其频谱图SF=fft(TS,NFFT)/L;%进行fft变换f=Fs/2*linspace(0,1,NFFT/2+1);%频率轴subplot(2,2,4);plot(f,2*a
16、bs(SF(1:NFFT/2+1)))%画出滤波后信
