代码搜索结果

%Samp7_6 wp=[0.35 0.65];N=48; %通带边界频率（归一化频率）和滤波器阶数 Fs=50; b=fir1(N,wp); %设计FIR带通滤波器 figure(1) [H,f]=freqz(b,1,512,Fs); %以50Hz为采样频率求出滤波器频率响应 subplot(2,1,1),plot(f,20*log10(abs

% Samp3_17 N1=10;dt=1; %采样点数和采样间隔 n1=0:N1-1;t1=n1*dt; %时间序列 x=cos(2*pi*0.24*t1)+cos(2*pi*.26*t1); %原始信号 subplot(2,3,1);stem(t1,x); hold on; plot([0 10],[0 0]); %绘制横轴 title('信号x(n), n=0

close all;figure(1) ws=[1000 2000]*2*pi;wp=[500 2500]*2*pi;Rp=1;Rs=50; %滤波器设计参数，乘以2 转换为弧度/s [N,Wn]=cheb1ord(wp,ws,Rp,Rs,'s'); %求模拟滤波器的最小阶数和截止频率 w=linspace(1,3000,1000)*2*pi; %设置计算复数频率响应的频率点 [

%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