信号分析课程设计Word文件下载.docx
- 文档编号:6722646
- 上传时间:2023-05-07
- 格式:DOCX
- 页数:22
- 大小:158.61KB
信号分析课程设计Word文件下载.docx
《信号分析课程设计Word文件下载.docx》由会员分享,可在线阅读,更多相关《信号分析课程设计Word文件下载.docx(22页珍藏版)》请在冰点文库上搜索。
%录入参数
T2=t2*0.004;
ff1=20;
ff2=45;
w1=2*pi*ff1;
w2=2*pi*ff2;
wo=(w1+w2)/2;
ow=(w2-w1)/2;
h0=2*(w2-w1)/pi;
hn=2*cos(wo.*T2).*sin(ow.*T2)./(pi*T2);
h=[h0,hn];
%褶积运算
y=conv(x,h);
%用fft通过y(t)产生Y(k)
Y=fft(y,256);
IYI=abs(Y);
%出图
subplot(411);
plot(x);
title('
x(t)'
);
subplot(412);
plot(h);
h(t)'
subplot(413);
plot(y);
y(t)'
subplot(414);
plot(IYI);
|y(t)|'
(2)频域滤波:
%频域滤波,产生x(t)
%录入参数
t1=1:
T1=t1*0.004;
s=exp(-b*T1.*T1).*sin(2*pi*f1*T1)+exp(-b*T1.*T1).*sin(2*pi*f2*T1);
%分别表达出s(t),p(t)
p=0.5*sin(2*pi*f3*T1)+0.5*cos(2*pi*f4*T1);
%产生X(k)
X=fft(x,256);
%补齐FFT
IXI=abs(X);
%产生H(k)
df=1/(256*0.004);
m1=ceil(20/df);
m2=ceil(45/df);
m3=256-m2;
m4=256-m1;
h1=zeros(1,(m1-1));
h2=ones(1,(m2-m1+1));
h3=zeros(1,(m3-m2-1));
h4=ones(1,(m4-m3+1));
h5=zeros(1,256-m4);
H=[h1,h2,h3,h4,h5];
IHI=abs(H);
Y=X.*H;
subplot(311);
plot(IXI);
|X(k)|'
subplot(312);
plot(IHI);
|H(k)|'
subplot(313);
|Y(k)|'
(3)通过ifft反求y(t):
代码如下:
%产生s(t)
%产生s(t)
%产生p(t)
y=ifft(Y,256);
h=ifft(H,256);
y1=conv(x,h);
%通过IFFT产生y
%出图:
subplot(311)
plot(real(y));
plot(real(y1));
y1(t)'
图标如下:
第二题振幅谱补零
%产生x(n)
A1=1;
A2=0.8;
A3=1;
N=200;
f1=35;
f2=45;
f3=60;
N;
x=exp(-b*T.*T).*cos(2*pi*f1*T)+0.8*exp(-b*T.*T).*cos(2*pi*f2*T)+exp(-b*T.*T).*cos(2*pi*f3*T);
%表达出x(n)
subplot(211);
x(n)'
%通过FFT运算出X(k)
X=fft(x);
subplot(212);
%补零:
2的八次方,补56个零,产生X1(k)
zero=zeros(1,56);
x1=[x,zero];
plot(x1);
x1(n)'
X=fft(x1);
|X1(k)|'
2的九次方,补312个零,产生X2(k)
zero=zeros(1,312);
x2=[x,zero];
plot(x2);
x2(n)'
X=fft(x2);
|X2(k)|'
2的十次方,补824个零,产生X3(k)
zero=zeros(1,824);
x3=[x,zero];
plot(x3);
x3(n)'
X=fft(x3);
|X3(k)|'
第三题计算自互相关
%带入参数
m=1:
150;
M=m*0.004;
x=1*sin(2*pi*f1*M);
%出图x(n)
%产生y(n)
f2=50;
m1=1:
M1=m1*0.004;
y=0.8*exp(-100*M1.*M1).*cos(2*pi*f2*M1);
%出图y(n)
y(n)'
%产生rxy(τ)
rxy=xcorr(x,y);
%调用函数xcorr
subplot(4,1,1);
plot(rxy);
rxy(τ)'
%产生ryx(τ)
ryx=xcorr(y,x);
subplot(4,1,2);
plot(ryx);
ryx(τ)'
%产生rxy(-τ)
%产生x(-n)
m2=1:
M2=-m2*0.004;
x2=1*sin(2*pi*f1*abs(M2));
%表达出x(-n)
subplot(4,1,3);
x(-n)'
rx2y=conv(x2,y);
subplot(4,1,4);
plot(rx2y);
rxy(-τ)'
自相关:
%产生rxx(τ)
rxx=xcorr(x,x);
subplot(2,1,1);
plot(rxx);
rxx(τ)'
%产生ryy(τ)
ryy=xcorr(y,y);
subplot(2,1,2);
plot(ryy);
ryy(τ)'
选做题地震记录及信噪比分析
(1)
%产生b(t)
f=50;
%参数录入
N=1:
41;
T=N*0.004;
b=(1-2*pi*pi*f*f*T.*T).*exp(-pi*pi*f*f*T.*T);
B=fft(b,64);
IBI=abs(B);
U=angle(B);
;
plot(b);
title('
b(t)'
plot(IBI);
B(w)'
plot(U);
U(w)'
(2)
%产生函数x(t)
g=zeros(1,100);
%产生g(t)
g(35)=-0.5;
g(45)=0.4;
g(48)=-0.35;
g(65)=0.5;
s=conv(b,g);
n=0.4*rand(1,140)-0.2;
x=s+n;
p1=sum(s.*s);
%产生信噪比
p2=sum(n.*n);
p=p1/p2;
fprintf('
%d'
p);
plot(real(x));
title('
plot(n);
n(t)'
plot(g);
g(t)'
%产生函数x(t)
n=1*rand(1,140)-0.5;
1.631078e-002
学生姓名:
XXX
学号:
xxxxxxxxxxx
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 分析 课程设计