代码搜索结果

function [Hr,w,a,L] = Hr_Type1(h); % 计算第一种低通滤波器设计的振幅响应Hr(w) % ----------------------------------------------------------- % [Hr,w,a,L] = Hr_Type1(h) % Hr = 振幅响应 % w = 在[0 pi] 区间内计算Hr 的500个频率点

function [Hr,w,d,L] = Hr_Type4(h); % 计算第四种低通滤波器设计的振幅响应Hr(w) % --------------------------------------------------- % [Hr,w,d,L] = Hr_Type4(h) % Hr = 振幅响应 % w = 在[0 pi] 区间内计算Hr 的500个频率点 % d =

function [Hr,w,c,L] = Hr_Type3(h); % 计算第三种低通滤波器设计的振幅响应Hr(w) % ----------------------------------------------------------- % [Hr,w,c,L] = Hr_Type3(h) % Hr = 振幅响应 % w = 在[0 pi] 区间内计算Hr 的频率点 %

wp=0.2;ws=0.4;Rp=1;As=20; % 定义滤波器参数(可修改) [n,wc]=buttord(wp,ws,Rp,As); % 计算滤波器阶次和3dB截止频率 [b,a]=butter(n,wc); % 设计巴特沃斯滤波器 w=[0:pi/200:pi];

function y=sinc(x) % 函数 y=sinc(x)=sin(pi*x)./(pi*x) % ----------------------- x = x + eps; y = sin(pi*x) ./ (pi*x);

% 生成三个正弦函数并计算其幅度谱 N=120; % 采样点数 fs=400; % 采样频率 n=1:N; % 离散时间 % 三个正弦函数和 x1=sin(2*pi*50*n/fs);x2=sin(2*pi*100*n/fs);x3=sin(2*pi*150*n/fs);x=x1+x2+x3; % 时域波形显示 figure(1);plot(n/fs,x);xlabe

%a21.m wp=0.2*pi;wr=0.3*pi;Ap=1;Ar=15; T=1; Omegap=(2/T)*tan(wp/2);Omegar=(2/T)*tan(wr/2); ep=sqrt(10^(Ap/10)-1); Ripple=sqrt(1/(1+ep*ep)); Attn=1/(10^(Ar/20)); %双线性变换法设计Chebshev 1型数字低通滤波器 [cs

function[Hr,w,d,L]=hr_type4(h); M=length(h); L=M/2; d=2*[h(L:-1:1)]; n=[1:1:L];n=n-0.5; w=[0:1:500]'*2*pi/500; Hr=sin(w*n)*d';

%a09.m n=[0:1:99];x=cos(0.48*pi*n)+cos(0.52*pi*n); n1=[0:1:9];y1=x(1:1:10); subplot(3,2,1);stem(n1,y1);title('信号x(n),0

%a25.m M=21;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l; Hrs=[1,1,1,1,zeros(1,14),1,1,1]; Hdr=[1,1,0,0];wdl=[0,0.35,0.35,1]; k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1; angH=[-alpha*(2*pi)/M*k1,alpha*(2