代码搜索结果

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, '

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));

%截短对升余弦信号频谱的影响 M=128; Rb=1000; a=0.5; t=-5/Rb:0.0003:5/Rb; t1=cos(a*Rb*pi*t); t2=sin(pi*Rb*t)./(pi*Rb*t); t3=1-4*a*a*Rb*Rb*t.*t; P=t2.*t1./(t3); Y=fft(P,M); subplot(3,1,1),plot(t,P) ; grid

#include #include void main(void) { double pi = 3.14159265; printf("Tangent of pi is %f\n", tan(pi)); printf("Tangent of pi/4 is %f\n", tan(pi / 4.0)); }

% dbp.m - ideal bandpass FIR filter % % h = dbp(wa, wb, N) = row vector % % N = 2M+1 = filter length (odd) % wa, wb = cutoff frequencies in [rads/sample] function h = dbp(wa, wb, N) M = (N-1)/2; fo

%Funcion que evalua la matriz Ai. Se debe incresar en el orden %correspondiente los siguientes valores: %a=Angulo de giro alrededor del eje Z %b=Desplazamiento en el eje Z %c=Desplazamiento en el

#include #include void main(void) { double pi = 3.14159265; printf("Tangent of pi is %f\n", tan(pi)); printf("Tangent of pi/4 is %f\n", tan(pi / 4.0)); }

function [A0,A,B]=FZZ(func,T, n) syms t; func = subs(sym(func), findsym(sym(func)),sym('t')); A0=int(sym(func),t,-T/2,T/2)/T; for(k=1:n) A(k)=int(func*cos(2*pi*k*t/T), t,-T/2,T/2)*2/T; A