修改以下M文件程序并给出修改过后的代码要求：1 将显示图像代码进行修改显示出与原来代码不一样的图像2 将代码中kasier窗设计代码修改为矩形窗的设计代码3 将代码中Min-Max design 参数进行修改与原来代码的参数不一样代码如下：clc;clear all; close all;loadexam_tsmat;x=h_ts_data;x_1st=x1;x_3rd=x3;Ndot=lengt

clc; clear all; close all; load('exam_ts.mat'); x=h_ts_data; x_1st=x(1,:); x_3rd=x(3,:); Ndot=length(x_1st); fs=max(freq_axis)2; Nfft=512; %1 [PXX_1ST,F]=periodogram(x_1st,hamming(Ndot),Ndot,fs); [PXX_3RD,F]=periodogram(x_3rd,hamming(Ndot),Ndot,fs); %2 design a filter fcuts=[0.05fs/2 0.1fs/2]; mags=[0 1]; devs=[0.02 0.01]; %2.1 rectangular window Nfilter1=round(1.5(fs/fcuts(1))); h_rect=fir1(Nfilter1,fcuts/(fs/2),'bandpass'); %2.2 Min-Max design devs=[0.02 0.01]; [Nfilter2,fpts,mag,wt]=remezord(fcuts,mags,devs,fs); b=remez(Nfilter2,fpts,mag,wt); %3.filter y_1st_re=filter(h_rect,1,x_1st); y_1st_mm=filter(b,1,x_1st); y_3rd_re=filter(h_rect,1,x_3rd); y_3rd_mm=filter(b,1,x_3rd); %3.1 spectrum [PXX_Y_1ST_re,F]=periodogram(y_1st_re,hamming(Ndot),Ndot,fs); [PXX_Y_1ST_mm,F]=periodogram(y_1st_mm,hamming(Ndot),Ndot,fs); [PXX_Y_3RD_re,F]=periodogram(y_3rd_re,hamming(Ndot),Ndot,fs); [PXX_Y_3RD_mm,F]=periodogram(y_3rd_mm,hamming(Ndot),Ndot,fs); %3.2 TDFT FT_Y_1ST_re=fftshift(fft(y_1st_re,Nfft)); FT_Y_1ST_mm=fftshift(fft(y_1st_mm,Nfft)); FT_Y_3RD_re=fftshift(fft(y_3rd_re,Nfft)); FT_Y_3RD_mm=fftshift(fft(y_3rd_mm,Nfft));

figure(1) plot(freq_axis,10log10(fftshift(PXX_1ST)),'b',freq_axis,10log10(fftshift(PXX_3RD)),'r'); legend('spectra of 1st range samples','spectra of 3rd range samples') grid; xlabel('frequency(Hz)'); ylabel('power spectrum density(dB)'); velocity=150*nyquist_velocity/max(freq_axis)

figure(2) freqz(h_rect); title([num2str(Nfilter1),'th notch filter based on rectangular window'])

figure(3) freqz(b); title([num2str(Nfilter2),'th notch filter based on Min-max design'])

figure(4) subplot(411) plot(freq_axis,10log10(fftshift(PXX_Y_1ST_re)),'b',freq_axis,10log10(fftshift(PXX_Y_3RD_re)),'r'); legend('spectra of 1st range samples filterd','spectra of 3rd range samples filterd') title('spectra of 1st and 3rd range samples after rectangular filter') grid; xlabel('frequency(Hz)'); ylabel('PSD(dB)'); subplot(412) plot(freq_axis,10log10(fftshift(PXX_Y_1ST_mm)),'b',freq_axis,10log10(fftshift(PXX_Y_3RD_mm)),'r'); legend('spectra of 1st range samples filterd','spectra of 3rd range samples filterd') title('spectra of 1st and 3rd range samples after Min-max filter') grid; xlabel('frequency(Hz)'); ylabel('PSD(dB)');

subplot(413) plot(freq_axis2/fs,abs(FT_Y_1ST_re),'b',freq_axis2/fs,abs(FT_Y_3RD_re),'r'); legend('TDFT of 1st range samples filterd','TDFT of 3rd range samples filterd') title('TDFT of 1st and 3rd range samples after rectangular filter') grid; xlabel('\omega(pi/rad)'); ylabel('H(e^{j{\omega}})'); subplot(414) plot(freq_axis2/fs,abs(FT_Y_1ST_mm),'b',freq_axis2/fs,abs(FT_Y_3RD_mm),'r'); legend('TDFT of 1st range samples filterd','TDFT of 3rd range samples filterd') title('TDFT of 1st and 3rd range samples after Min-max filter') grid; xlabel('\omega(pi/rad)'); ylabel('H(e^{j{\omega}})')