fgp711.m

数字信号处理教程MATLAB释义与实现里面大量的例程

%《数字信号处理教程——MATLAB释义与实现》第七章演示程序fgp711
% 滤波器绝对指标与相对指标之间的变换
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
delta1=[0:0.01:0.3];                        % 设定delta1作为自变量数组
for delta2=[0.001,0.003,0.01,0.03,0.1]      % 对若干个delta2循环计算
    Rp=-20*log10((1-delta1)./(1+delta1));   % 算Rp
    As=-20*log10(delta2./(1+delta1));       % 算As
    subplot(1,2,1),plot(delta1,Rp),         % 画Rp图
    hold on,grid on
    subplot(1,2,2),plot(delta1,As),hold on  % 画As图
end
set(gcf,'color','w'),grid on                % 置图形背景色为白
% 以下为按用户输入的指标值进行换算的程序段
disp('本程序先出了由delta1,delta2查Rp,As的曲线族')
disp('如想要精确结果，本程序将依次要求用户输入Rp,As，求delta1,delta2；')
disp('然后要求用户输入delta1,delta2；求Rp,As的精确值。')
Rp=input('Rp=  '),delta1=(1-10^(-Rp/20))./(1+10^(-Rp/20)) 	%	由Rp求delta1的语句	
As=input('As=  '),delta2=(1+delta1).*10^(-As/20) 			%	由As及delta1求delta2的语句
delta1=input('delta1=  '),
Rp=-20*log10((1-delta1)./(1+delta1));		% 	由delta1及delta2求Rp的语句
delta2=input('delta2=  '),
As=-20*log10(delta2./(1+delta1)); 		%	由delta1及delta2求As的语句