资源描述：
《中国海洋大学声学基础实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、脉冲频率为1khz，周期10ms，N=1的正弦脉冲串。Fp=2500hz;Fs=3000hz;rp=1db,as=30db,前两个需要13阶，椭圆只需4阶。相关系数0.99，时间间隔104ms.开始50hz，结束500hz,采样时间500ms,采样率10000hz。clc;clearall;closeall;loadfish2;sound(fish2,FS);t=linspace(0,length(fish2)/FS,length(fish2));plot(t,fish2);xe1=fish2(9601:9601+0.07*FS);xe2
2、=fish2(14401:14401+0.07*FS);y1=abs(fft(xe1));y2=abs(fft(xe2));ff=linspace(0,FS,length(xe1));figure(1)plot(ff(1:350),y1(1:350),'r',ff(1:350),y2(1:350),'b');xlabel('f/hz'),ylabel('magnitude'),title('两次鱼叫的频谱');coe=corrcoef(y1,y2);y=xcorr(xe1,fish2);figure(2),plot(y);pausenum
3、1=find(y==max(y(0.5*10^4:0.8*10^4)));num2=find(y==max(y(1*10^4:1.2*10^4)));dt=(num2-num1)*1/FS%声学基础实验1clearall;closeall;clcf1=50;f2=20*10^3;fs=5*f2;angle1=pi/4;angle2=pi/3;t=[1/fs:1/fs:0.2];x=10*sin(2*pi*f1*t+angle1)+sin(2*pi*f2*t+angle2)+2;%去掉直流偏移x1=x-mean(x);figure(1);p
4、lot(t,x,'r',t,x1,'b');xlabel('time/s'),ylabel('magnitude'),legend('原始','去直流偏移'),title('signalx=10*sin(2*pi*f1*t+angle1)+sin(2*pi*f2*t+angle2)+2');%滤波器设计Fs=20*10^3;Fp=50;[N,wc]=buttord(2*Fp/fs,2*Fs/fs,1,60);[B,A]=butter(N,wc);[B2,A2]=butter(N+2,wc);[hf,ff2]=freqz(B,A,1024,
5、fs);[hf2,ff2]=freqz(B2,A2,1024,fs);[N2,wc2]=buttord(2*Fs/fs,2*Fp/fs,1,60);[B3,A3]=butter(N2,wc2,'high');[hf3,ff]=freqz(B3,A3,1024,fs);figure(2),plot(ff2,20*log10(abs(hf)/abs(hf(1))),'r',ff2,20*log10(abs(hf2)/abs(hf2(1))),'b');title('butterworth低通滤波器频响曲线');grid;xlabel('f/h
6、z');ylabel('magnitude/db');legend('2阶','4阶')axis([0,Fs,-60,5]);line([0,Fs],[-3,-3]);line([wc*fs/2,wc*fs/2],[-60,5]);%滤波作频谱x2=filter(B,A,x1);x3=filter(B3,A3,x1);figure(3)subplot(2,1,1),plot(t(1:4000),x2(1:4000),'r'),xlabel('t/s');ylabel('magnitude');title('低通滤波后的时域波形');sub
7、plot(2,1,2),plot(t(1:200),x3(1:200),'b'),xlabel('t/s');ylabel('magnitude');title('高通滤波后的时域波形');y=abs(fft(x1));y1=abs(fft(x2));y3=abs(fft(x3));ff=linspace(1,fs,length(x1));figure(4),subplot(3,1,1),plot(ff(1:4500),y(1:4500));axis([10,2.5*10^4,0,3*10^4]),xlabel('f/hz'),ylabe
8、l('magnitude'),title('原信号频谱');subplot(3,1,2),plot(ff(1:4500),y1(1:4500));axis([10,2.5*10^4,0,3*1
