代码搜索结果

%Samp2_8 x=-10:0.1:10; %给出自变量序列值 y=sin(pi*x+eps)./(pi*x+eps); %给出sinc函数的离散序列值，eps是MATLAB系统的精度，这里防止被零除 plot(x,y) %绘图 xlabel('t') ylabel('sinc(t)')

%Samp6_8 Fs=1000; %采样频率 wp=300*2/Fs;ws=200*2/Fs; %根据采样频率将滤波器边界频率进行转换（6-19式） Rp=1;Rs=20; %通带波纹和阻带衰减 Nn=128; %显示滤波器频率特性的数据长度 [N,Wn]=buttord(wp,ws,Rp,Rs); %求得数字滤波器的最小阶数和截止频率（归一化频率） [b,a

%Samp6_5 Wp=2000*2*pi;Ws=3000*2*pi; %滤波器截止频率 Rp=3;Rs=15; %通带波纹和阻带衰减 Fs=10000; %采样频率 Nn=128; %调用freqz所用的频率点数 [N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); %模拟滤波器的最小阶数 [z,p,k]=buttap(N); %设计模拟低

%Samp10_2 Fs=100;dt=1/Fs; %给定模拟输入信号采样间隔 f1=6;f2=10;f3=19; %模拟输入信号的频率成分 N=500; %数据点数 t=[0:N-1]*dt; %模拟输入信号的时间序列 x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t)+0.5*sin(2*pi*f3*t); %模拟输入信号 X

%Samp5_10 wp=200*pi;ws=300*pi;Rp=1;Rs=16;%滤波器的通带阻带边界频率、通带波纹和阻带衰减 N=ceil(log10((10^(Rp/10)-1)/(10^(Rs/10)-1))/(2*log10(wp/ws))) % （5-22）式 Wcp=wp/((10^(Rp/10)-1)^(1/(2*N))) %（5-23）式 Wcs=ws/((10^(R

%Samp1_4 N=400;dt=0.01;f1=10;f2=9; %采样点数、间隔和两个振动的频率 n=0:N-1;t=n*dt; %定义时间离散值 x1=.3*sin(2*pi*f1*t+1); %第一个振动 x2=.3 *sin(2*pi*f2*t+1); %第二个振动 subplot(3,1,1),plot(t,x1); title('第一个振动') %第

%Samp5_15 wp=[1000 2000]*2*pi;ws=[500 2500]*2*pi;Rp=1;Rs=100; %滤波器设计参数，对于给定Hz应乘以2 [N,Wn]=buttord(wp,ws,Rp,Rs,'s'); %求得滤波器的最小阶数和截止频率 [b,a]=butter(N,Wn,'s'); %设计模拟Butterworth滤波器 w=linspace(1,3

%Samp5_19 N=5; %滤波器阶数 Wn=1000; %边界频率 [b,a]=besself(N,Wn); %设计Bessel滤波器 figure(1) [H,w]=freqs(b,a,512); %绘出滤波器的复数频率特性 magH=abs(H);phaH=unwrap(angle(H)); %求幅频响应和相频响应 subplot(2,1,1),plot(w/(2

%Samp7_5 wp=0.5*pi;ws=0.66*pi; %滤波器的边界频率 wdelta=ws-wp; %过渡带宽度 N=ceil(8* pi/wdelta);%求解滤波器的最小阶数，根据Hanning 窗主瓣宽 Wn=(0.5+0.66)*pi/2;%截止频率取通带和阻带边界频率的中点 b=fir1(N,Wn/pi,hanning(N+1));%

%Samp6_6 wp=0.2*pi;ws=0.3*pi;Rp=1;Rs=15; %数字滤波器截止频率、通带波纹和阻带衰减 T=0.01;Nn=128; %采样间隔 Wp=wp/T;Ws=ws/T; %得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式 [N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); %计算模拟滤波器的最小阶数 [z,p,k]=but