matlab实验程序例子Word文档下载推荐.docx
- 文档编号:7593721
- 上传时间:2023-05-08
- 格式:DOCX
- 页数:19
- 大小:727.22KB
matlab实验程序例子Word文档下载推荐.docx
《matlab实验程序例子Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《matlab实验程序例子Word文档下载推荐.docx(19页珍藏版)》请在冰点文库上搜索。
stem(n,imag(xn))
stem(n,angle(xn))
stem(n,xn)
实部:
虚部:
幅度:
相位角:
n=-10:
10;
xn=exp(-0.1+0.3*j)*n;
幅度:
相位:
X(n):
x=zeros(1,21);
x1=zeros(1,21);
x2=zeros(1,21);
forn=-10:
1:
10
switchn
case0
x(n+11)=2;
case2
x(n+11)=1;
case3
x(n+11)=-1;
case4
x(n+11)=3;
otherwise
x(n+11)=0;
end
forn=1:
21
ifn-3>
x1(n)=x(n-3);
x2(n+11)=x(-n+11);
n=-10:
subplot(1,3,1);
stem(n,x);
title('
x[n]'
);
xlabel('
n'
ylabel('
x'
axis([-1010-33]);
gridon;
subplot(1,3,2);
stem(n,x1);
x[n-3]'
x1'
subplot(1,3,3);
stem(n,x2);
x[-n]'
x2'
x=[1,2,3,4];
y=[1,1,1,1,1];
lenx=length(x);
leny=length(y);
n=length(x)+length(y)-1;
Xk=fft(x,N);
Xk=fft(x,n);
Yk=fft(y,n);
Tk=Xk.*Yk;
t=real(ifft(Tk));
plot(t);
k=0:
a=[1,4,2,1];
b=[1];
h=impz(b,a,k);
subplot(1,2,1)
stem(k,n)
x=sym('
cos(n)'
X=ztrans(x)
X=
(z-cos
(1))*z/(z^2-2*z*cos
(1)+1)
X=sym('
1/((1+z)*(2+z))'
x=iztrans(X)
x=
1/2*charfcn[0](n)-(-1)^n+1/2*(-2)^n
•零极点并画出零极点图;
•系统的单位脉冲响应h(n)和频率响应
b=[2,1];
a=[1,2,1];
[r,p,k]=tf2zp(b,a)
zplane(b,a)
r=-0.5000
p=-1
-1
k=2
num=[0,2,1];
den=[1,2,1];
h=impz(num,den);
figure
(1);
stem(h)
[H,w]=freqz(num,den);
figure
(2);
plot(w/pi,abs(H))
X=0;
w=-2*pi:
pi/100:
2*pi;
forn=0:
X=X+0.9*exp(j*pi/3*n)*exp(-j*w*n);
magX=abs(X);
angleX=angle(X);
plot(w,X);
\omega'
DTFT[x(n)]'
plot(w,magX);
|DTFT[x(n)]|'
figure(3);
plot(w,angleX);
angleDTFT[x(n)]'
运行结果如下(上图为DTFT,左图为DTFT幅频特性,右图为DTFT相频特性):
X=X+2^n*exp(-j*w*n);
50;
x=cos(0.82*pi*n)+2*sin(pi*n);
stem(n,real(dft(x,51)))
stem(n,angle(dft(x,51)))
其中dft函数:
functionXk=dft(x,N)
n=[0:
N-1];
k=n;
Wn=exp(-j*2*pi/N);
nk=n'
*k;
Wnnk=Wn.^nk;
Xk=x*Wnnk;
function[A]=run_time(m)
N=4^m;
n=0:
N-1;
xn=sin(0.2*pi*n)+randn(1,N);
dft_time=0;
fft_time=0;
t1=clock;
fft(xn);
fft_time=etime(clock,t1)
t2=clock;
dft(xn,N);
dft_time=etime(clock,t2)
run_time
(2)
fft_time=0
dft_time=0
run_time(3)
run_time(4)
dft_time=0.1400
4.研究高密度频谱和高分辨率频谱。
设有连续信号
•以采样频率32kHzsf=对信号x(t)采样,分析下列三种情况的幅频特性。
•采集数据长度N=16点,做N=16点的FFT,并画出幅频特性。
•采集数据长度N=16点,补零到256点,做N=256点的FFT,并画出幅频特性。
•采集数据长度N=256点,做N=256点的FFT,并画出幅频特性。
观察三种不同频率特性图,分析和比较它们的特点以及形成的原因。
Ts=1/32000;
x1(n+1)=cos(2*pi*6500*n*Ts)+cos(2*pi*7000*n*Ts)+cos(2*pi*9000*n*Ts);
x2(n+1)=cos(2*pi*6500*n*Ts)+cos(2*pi*7000*n*Ts)+cos(2*pi*9000*n*Ts);
forn=16:
255
x2(n+1)=0;
255;
x3(n+1)=cos(2*pi*6500*n*Ts)+cos(2*pi*7000*n*Ts)+cos(2*pi*9000*n*Ts);
Xk1=fft(x1,16);
stem(abs(Xk1));
采集数据长度N=16点的FFT幅频特性'
Xk2=fft(x2,256);
stem(abs(Xk2));
采集数据长度N=16点,补零到256点的FFT幅频特性'
Xk3=fft(x3,256);
stem(abs(Xk3));
采集数据长度N=256点的FFT幅频特性'
运行结果如下:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 实验 程序 例子