代码搜索结果

close all;clear all;clc; t=0:0.01:4*pi; f=2; a=2; ut1=a*exp(j*2*pi*f*t); ut2=a*exp(-j*2*pi*f*t); ut3=a*exp(j*2*pi*f*t+pi/2); figure;subplot(3,1,1),plot(t,ut1); subplot(3,1,2),plot(t,ut2); sub

clear all close all clc M=32; alpha=(M-1)/2; magHk=[1,1,1,0.5,zeros(1,25),0.5,1,1]; k1=0:15; k2=16:M-1; angHk=[-alpha*2*pi/M*k1,alpha*2*pi/M*(M-k2)]; H=magHk.*exp(j*angHk); h=real(ifft(H,M))

% x = [1,1,1,1]; subplot(1,1,1); % % a) DTFT w = [0:1:500]*2*pi/500; [H] = freqz(x,1,w); magH = abs(H); phaH = angle(H); phaH(126)=-47.5841*pi/180; subplot(2,1,1); plot(w/pi,magH); grid xlabe

clear all; close all; clc; N=32; n=0:N-1; xn=cos(pi*n/6); k=0:N-1; WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^nk; Xk=xn*WNnk; figure(1) stem(n,xn) figure(2) stem(k,abs(Xk));

nfft = 1024; fs = 1000; t = 0:1/fs:1; xn = sin(2*pi*100*t)+sin(2*pi*300*t)+sin(2*pi*450*t)+randn(size(t)); [P, f] = pburg(xn, 18, nfft, fs); Pxx = 10 * log10(P); figure plot(f, Pxx); grid on;

randn('state', 1); n = 0:99; s = sin(pi/3*n) + 2*sin(pi/4*n) + randn(1, 100); X = corrmtx(s, 7, 'mod'); % 使用改进的协方差法估计互相关矩阵 [s, w, V, E] = peig(X, 4, 'whole'); peig(X, 4, 'whole');

Fs = 1000; % 采样频率 NFFT = 256; t = 0:1/Fs:1; 　% 采样时间 x = sin(2*pi*100*t) + sin(2*pi*200*t) + sin(2*pi*400*t) + randn(size(t)); w = hanning(NFFT)'; % 256个点的海宁窗，注意取转置 Pxx = (abs(fft(w.*x(1:NFFT

randn('state', 1); n = 0:99; s = exp(i*pi/2*n) + 2*exp(i*pi/4*n) + exp(i*pi/3*n) + randn(1, 100); X = corrmtx(s, 7, 'mod') % 使用改进的协方差法估计互相关矩阵 [s, w, V, E] = pmusic(X, 4, 'whole'); pmusic(X, 4, '

% 频率采样技术: 差分器 % M = 33; alpha = (M-1)/2; Dw = 2*pi/M; l = 0:M-1; wl = Dw*l; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; Hrs = [j*Dw*k1,-j*Dw*(M-k2)]; angH = [-alpha*Dw*k1, alpha*Dw*(M-k2

% 用PM算法进行的高通滤波器设计 % ws = 0.6*pi; wp = 0.75*pi; Rp = 0.5; As = 50; delta1 = (10^(Rp/20)-1)/(10^(Rp/20)+1); delta2 = (1+delta1)*(10^(-As/20)); deltaH = max(delta1,delta2); deltaL = min(delta1,de