代码搜索结果

% example 9 x=linspace(0,2*pi,30); y=sin(x); z=cos(x); plot(x,y,'g:',x,z,'r-',x,y,'bo',x,y,'c+')

n=0:100;x=cos(pi*n/2); k=-100:100;w=(pi/100)*k; X=x*(exp(-j*pi/100)).^(n'*k); y=exp(j*pi*n/4).*x; Y=y*(exp(-j*pi/100)).^(n'*k); subplot(2,2,1);plot(w/pi,abs(X)); xlabel('frequency');ylabel('|X|'

Wp=0.2*pi;Ws=0.3*pi;Rp=1;As=16; Ripple=10^(-Rp/20);Attn=10^(-As/20); %Analog filter design: [b,a]=afd_chb1(Wp,Ws,Rp,As) %Chebyshev-1 Filter Order=4 % Calculation of second-order sections: [C,B,A

% eg 1-6 n=0:50; x=3*cos(0.1*pi*n+pi/3); close all stem(n,x); title('余弦序列'); xlabel('n');ylabel('x(n)');

x=rand(1,11);n=0:10; k=0:500; w=(pi/500)*k; X=x*(exp(-j*pi/500)).^(n'*k); [y,m]=sigshift(x,n,2); Y=y*(exp(-j*pi/500)).^(m'*k); error=max(abs(Y-exp(-j*2*w).*X)); disp(error);

function hd=ideal_lp(wc,M) alpha=(M-1)/2; n=[0:1:(M-1)]; m=n-alpha+eps; hd=sin(wc*m)./(pi*m);

n=[-1:3]; x=1:5; k=0:500;w=(pi/500)*k; X=x*exp(-j*pi/500.*(n'*k)); magX=abs(X);angX=angle(X); realX=real(X);imagX=imag(X); subplot(2,2,1);plot(w/pi,magX);grid; xlabel('Frequency');ylabel('ampli

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

wp=0.2*pi;ws=0.3*pi; tr_width=ws-wp; M=ceil(6.6*pi/tr_width) n=[0:1:M-1]; wc=(ws+wp)/2; hd=ideal_lp(wc,M); w_ham=(hamming(M))'; h=hd.*w_ham; [db,mag,pha,grd,w]=freqz_m(h,[1]); delta_w=2*pi/

% example 12 x=linspace(0,2*pi,30); y=sin(x); z=cos(x); a=2*sin(x).*cos(x); b=sin(x)./(cos(x)+eps); subplot(2,2,1)%选取4个子图中的左上 plot(x,y);axis([0 2*pi -1 1]),title('sin(x)'); subplot(2,2,2)%选取4个