matlab建模复习.docVIP

% ch3example1A.m 设计模拟低通滤波器 clear; f_p=2400; f_s=5000; R_p=3; R_s=25; [n, fn]=buttord(f_p,f_s,R_p,R_s, s); Wn=2*pi*fn; [b,a]=butter(n, Wn, s); f=0:100:10000; s=j*2*pi*f; H_s=polyval(b,s)./polyval(a,s); figure(1); subplot(2,1,1); plot(f, 20*log10(abs(H_s))); axis([0 10000 -40 1]); xlabel(频率 Hz);ylabel(幅度 dB); subplot(2,1,2); plot(f, angle(H_s)); xlabel(频率 Hz);ylabel(相角 rad); figure(2); freqs(b,a); % ch3example12prg1.m 计算频谱 Df=5; f_s=2*500; N=f_s/Df; t=0:1./f_s:(N-1)./f_s; freq=0:Df:(N-1)*Df; f_t=2*sin(2*pi*100*t)+cos(2*pi*180*t); F_f=1/f_s*fft(f_t,N); plot(freq-f_s/2,abs(fftshift(F_f))); xlabel(频率 Hz); ylabel(幅度谱×); % ch3example16prg4.m 用海明窗计算功率谱 fs=500; Df=1; N=floor(fs/Df)+1; t=0:1/fs:(N-1)/fs; F=0:Df:fs; xk=sin(2*pi*50*t)+2*sin(2*pi*130*t)+randn(1,length(t)); w=hamming(167); w=w*sqrt(167/sum(w.*w)); Pxx=(abs(fft(w.*xk(1:167))).^2+... abs(fft(w.*xk(83:249))).^2+... abs(fft(w.*xk(168:334))).^2+... abs(fft(w.*xk(250:416))).^2+... abs(fft(w.*xk(335:501))).^2)/5/((N/3)^2); Pav_timedomaim=sum(xk.^2)/N Pav_freqdomain=sum(Pxx) plot(0:3:fs,10*log10(Pxx));xlabel(freq Hz);ylabel(PSD dB); % ch5example2prg1.m 信噪比性能曲线 SNR_in_dB=-10:2:30; SNR_in=10.^(SNR_in_dB./10); m_a=0.3; P=0.5+(m_a^2)/4; for k=1:length(SNR_in) sigma2=P/SNR_in(k); sim(ch5example2.mdl); SNRdemod(k,:)=SNR_out; end plot(SNR_in_dB, SNRdemod); xlabel(输入信噪比 dB); ylabel(解调输出信噪比 dB); legend(包络检波,相干解调); Zero-Order Hold （采样信号） Quantizer（量化） Interger to Bit Converter（二进制编码） Frame StatusConversion （并串转换） Buffer Reshape （串并转换） Bit to Integer Conventer （D/A转换） Vector Scope