




已阅读5页,还剩8页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
DSP考试实验报告 姓名学号专业班级电信提交日期2014年12月18日 一、频分复用1 选择3个不用频段的信号,画出其频谱%1.150Hz,300Hz450Hz t=0:0.001:1;x1=cos(2*pi*150*t);x2=cos(2*pi*300*t);x3=cos(2*pi*450*t);y1=fft(x1,256);y2=fft(x2,256);y3=fft(x3,256);f=1000*(0:128)/256;figure(1);subplot(3,1,1);plot(f,abs(y1(1:129);xlabel(频率);ylabel(幅度);subplot(3,1,2);plot(f,abs(y2(1:129);xlabel(频率);ylabel(幅度);subplot(3,1,3);plot(f,abs(y3(1:129);2 将三路信号合成一路信号,画出合成信号的时域和频域波形%2.3路信号合成1路信号figure(2);x=x1+x2+x3;stem(t,x);subplot(2,1,1);stem(t,x);plot(x(1:128);xlabel(t);ylabel(y);title(时域);subplot(2,1,2);y=fft(x,256);plot(f,abs(y(1:129);xlabel(频率);ylabel(幅度);title(频域);3将合成信号分别通过设计好的数字滤波器,分离出原来的三路信号,分析得到的三路信号的时域波形和频谱,与原始的三路信号进行比较说明频分复用的特点4将所用滤波器的频谱画出来%低通t=0:0.001:1;x1=cos(2*pi*150*t);x2=cos(2*pi*300*t);x3=cos(2*pi*450*t);x=x1+x2+x3;y1=fft(x,256);fp=200;fs=1200;ft=4000;rp=1;rs=40; T=1/ft; wp=2*pi*fp/ft; ws=2*pi*fs/ft; omegap=(2/T)*tan(wp/2); omegas=(2/T)*tan(ws/2); N,wn=buttord(omegap,omegas,rp,rs,s); num,den=butter(N,wn,low,s); numz,denz=bilinear(num,den,ft); figure(3) freqz(numz,denz) y=filter(numz,denz,x);y2=fft(y,256);f=1000*(0:128)/256;subplot(4,1,3);plot(f,abs(y1(1:129);xlabel(频率);ylabel(幅度);title(原始信号频域);subplot(4,1,4);plot(f,abs(y2(1:129);xlabel(频率);ylabel(幅度);title(滤波后频域);%带通 t=0:0.001:1;x1=cos(2*pi*150*t);x2=cos(2*pi*300*t);x3=cos(2*pi*450*t);x=x1+x2+x3;y1=fft(x,256);fp1=900;fp2=1300;fs1=600;fs2=1500;ft=4000;rp=1;rs=40; T=1/ft; wp=2*pi*fp1 fp2/ft; ws=2*pi*fs1 fs2/ft; omegap=(2/T)*tan(wp/2); omegas=(2/T)*tan(ws/2); N,wn=buttord(omegap,omegas,rp,rs,s); num,den=butter(N,wn,bandpass,s); numz,denz=bilinear(num,den,ft); figure(5) freqz(numz,denz)y=filter(numz,denz,x);y2=fft(y,256);f=1000*(0:128)/256;subplot(4,1,3);plot(f,abs(y1(1:129);xlabel(频率);ylabel(幅度);title(原始信号频域);subplot(4,1,4);plot(f,abs(y2(1:129);xlabel(频率);ylabel(幅度);title(滤波后频域);%高通t=0:0.001:1;x1=cos(2*pi*150*t);x2=cos(2*pi*300*t);x3=cos(2*pi*450*t);x=x1+x2+x3;y1=fft(x,256);fp=1400;fs=1500;ft=4000;rp=1;rs=40; T=1/ft; wp=2*pi*fp/ft; ws=2*pi*fs/ft; omegap=(2/T)*tan(wp/2); omegas=(2/T)*tan(ws/2); N,wn=buttord(omegap,omegas,rp,rs,s); num,den=butter(N,wn,high,s); numz,denz=bilinear(num,den,ft); figure(4) freqz(numz,denz) y=filter(numz,denz,x);y2=fft(y,256);f=1000*(0:128)/256;subplot(4,1,3);plot(f,abs(y1(1:129);xlabel(频率);ylabel(幅度);title(原始信号频域);subplot(4,1,4);plot(f,abs(y2(1:129);xlabel(频率);ylabel(幅度);title(滤波后频域);频率复用系统的优点:信道复用率高,允许复用的路数多,同时它的分路也很方便。二、音频信号分析1读出所给音频信号,确定该音频信号的采样率=44100Hz和量化比特数=16.画出该音频信号的时域波形和频谱图2对该音频信号以间隔M=2进行采样,回放音频信号,感觉和原来有什么不同?N1=100000;N2=600000;x,fs,bits=wavread(C:UserslcDesktoplaughing.wav,N1,N2);wavwrite(x,44100,C:UserslcDesktoplaughingnew.wav);Nx=N1:N2;figure;subplot(2,1,1);plot(Nx,x(:,1);title(C 1);subplot(2,1,2);plot(Nx,x(:,2);title(C 2);figure;subplot(2,1,1);X1=fft(x(:,1);plot(Nx,fftshift(X1);title(C 1);subplot(2,1,2);X2=fft(x(:,2);plot(Nx,fftshift(X2);title(C 2);x1c=x(1:2:length(x),1);x2c=x(1:2:length(x),2);for i=1:length(x)/2xnew(i,1)=x1c(i);xnew(i,2)=x2c(i);endwavwrite(xnew,40000,C:UserslcDesktoplaughingnewoo.wav);采样后的音频信号比原来的速度更快,而且音调更尖锐。3你想怎么解决?为什么?在采样后的每个点中间插入相邻两点的平均值以恢复原信号。4画出解决前后音频信号的波形及频谱for i=1:2:length(x)-10xori(i,1)=x1c (i+1)/2);xori(i,2)=x2c (i+1)/2);xori(i+1,1)=0.5*(x1c (i+1)/2)+x1c (i+3)/2);xori(i+1,2)=0.5*(x2c (i+1)/2)+x2c (i+3)/2);endwavwrite(xori,40000,laughing_ori.wav);%4figure;subplot(3,1,1);plot(1:length(x)-1,x(1:length(x)-1,:);title(origin);subplot(3,1,2);plot(1:2:length(x)-1,xnew);title(Down Sampling);subplot(3,1,3);plot(1:length(x)-9,xori);title(Regain);figure;subplot(3,1,1);plot(1:length(x)-1,abs(fftshift(fft(x(1:length(x)-1,:);title(origi
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 01 第7讲 牛顿第一定律、牛顿第二定律 【答案】作业手册
- 保姆教育专业知识培训课件
- 保姆基本知识培训内容
- 人教版地理七下8.1《中东》教学设计
- 保健胸部专业知识培训课件
- 2025年10月“江南十校”2026届新高三第一次综合素质检测 语文试卷(含答案详解)
- 2025年中国茶道艺术NFT茶艺师专业技能提升课程模拟题及答案解析
- 2025年中国铁路总公司信息技术岗位招聘考试指南及模拟试题
- 2025年农村合作经济背景下的财务岗位应聘指南与模拟题解析
- 2025年市场营销策略制定实战模拟题集与解析
- TSG D2002-2006燃气用聚乙烯管道焊接技术规则
- GB/T 320-2025工业用合成盐酸
- 2024年公路水运工程助理试验检测师《水运结构与地基》考前必刷必练题库500题(含真题、必会题)
- 2025年社工招聘考试试题及答案
- 病理检查报告审核制度
- 《数据库系统概论》全套课件(南京农业大学)
- 山西建投集团考试真题
- JT-T-325-2018营运客运类型划分及等级评定
- 湘教版版八年级上册地理知识点复习总结
- 2069-3-3101-002WKB产品判定准则-外发
- 美国发布2024版《关键和新兴技术清单》(英)
评论
0/150
提交评论