代码搜索结果

/* * Example program * * Copyright (c) 1988-1997 Shamus Software Ltd. */ #include #include "miracl.h" int main() { /* Brents example program */ flash x,pi; mira

wp=0.2*pi;Rp=1;ws=0.3*pi;As=15; [n,wn]=buttord(wp,ws,Rp,As,'s') [z,p,k]=buttap(n); [bap,aap]=zp2tf(z,p,k) [b,a]=lp2lp(bap,aap,wn); [db,mag,pha,w]=freqs_m(b,a,0.5*pi); [ha,x,t]=impulse(b,a); fig

wp=0.2*pi;Rp=1;ws=0.3*pi;As=15; [n,wn]=buttord(wp,ws,Rp,As) [b,a]=butter(n,wn); [db,mag,pha,grd,w]=freqz_m(b,a,0.5*pi); figure(3); subplot(2,2,1),plot(w/pi,mag);title('幅度'); subplot(2,2,2),plot(

clc ws1=0.2*pi;wp1=0.35*pi; wp2=0.65*pi;ws2=0.8*pi; wd=min((wp1-ws1),(ws2-wp2)); wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2; w_bla=(blackman(N+1))'; N=ceil(11*pi/wd); hd=ideal_lp(wc2,N+1)-ideal_lp(wc1,N+1

N=40;T1=0.5925;T2=0.1099; ws1=0.2*pi;wp1=0.35*pi;ws2=0.8*pi;wp2=0.65*pi; alpha=(N-1)/2; L=0:N-1;w1=(2*pi/N)*L; hrs=[zeros(1,5),T2,T1,ones(1,7),T1,T2,zeros(1,9),T2,T1,ones(1,7),T1,T2,zeros(1,4)];

ws1=0.2*pi;wp1=0.35*pi; wp2=0.65*pi;ws2=0.8*pi; wd=min((wp1-ws1),(ws2-wp2)); wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2; w_bla=(blackman(N+1))'; N=ceil(11*pi/wd); hd=ideal_lp(wc2,N+1)-ideal_lp(wc1,N+1); h

N=40;T1=0.5925;T2=0.1099; ws1=0.2*pi;wp1=0.35*pi;ws2=0.8*pi;wp2=0.65*pi; alpha=(N-1)/2; L=0:N-1;w1=(2*pi/N)*L; hrs=[zeros(1,5),T2,T1,ones(1,7),T1,T2,zeros(1,4)]; k1=[0:floor((N-1)/2)]; k2=floor(

clc; ws1=0.2*pi;wp1=0.35*pi;wp2=0.65*pi;ws2=0.8*pi; wd=min((wp1-ws1),(ws2-wp2)); wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2; N=ceil(11*pi/wd); wn=[wc1,wc2]/pi; b=fir1(N,wn,blackman(N+1)); freqz(b,1)

function [Xk]=dfs(xn,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^nk Xk=xn*WNnk;

n=0:0.1:2*pi; x1=2*cos(2*pi*10*n/50); stem(n,x1);