MATLAB代做  | MATLAB代写  | MATLAB有偿编程
首 页
当前位置:MATLAB代做|MATLAB代写|MATLAB有偿编程matlab源码信号处理
matlab信号分析        [↓跳到下载地址]
官方主页:源码共享网站:www.buycode.cn

联系方式:客服QQ:1747812398 买代码源码软件站,欢迎咨询

运行环境:Win9X/2000/XP/2003/

源码语言:简体中文

源码类型:源码程序 - 信号处理

推荐星级:

更新时间:2088-08-04 22:32:33

源码简介

fn=1000;Ap=0.1;As=60;N=800;   %信号长度N为800

                

Fs=1000;T=1/Fs;Tp=N*T;   %采样频率Fs=10kHz,Tp为采样时间

t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;  

fc1=Fs/10;                %第1路调幅信号载波频率fc1=1000Hz

fm1=fc1/10;               %第1路调幅信号的调制信号频率fm1=100Hz

fc2=Fs/20;                %第2路调幅信号载波频率fc2=500Hz

fm2=fc2/10;               %第2路调幅信号的调制信号频率fm2=50Hz

fc3=Fs/40;                %第3路调幅信号载波频率fc3=250Hz

fm3=fc3/10;               %第3路调幅信号的调制信号频率fm3=25Hz

xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t);    %产生第1路调幅信号

xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t);    %产生第2路调幅信号

xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t);    %产生第3路调幅信号

st=xt1+xt2+xt3;        %三路信号相加,得到复合信号

fxt=fft(st,N);         %计算信号st的频谱

subplot(4,3,10);

plot(t,st);grid;xlabel('t/s');ylabel('s(t)');

axis([0,Tp,min(st),max(st)]);title('(a)s(t)的波形')

subplot(4,3,11);

stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b)s(t)的频谱')

axis([0,Fs/8,0,1.2]);

xlabel('f/Hz');ylabel('幅度');

低通

Wp=0.06*pi;

Ws=0.12*pi;

[E,wn]=buttord(Wp/pi,Ws/pi,Ap,As);              %计算巴特沃斯滤波器阶次和截止频率

[b,a]=butter(E,wn,'low');                         %频率变换法设计巴特沃斯高通滤波器

[db,mag,pha,grd,w]=freqz_m(b,a);                  %数字滤波器响应

subplot(4,3,1);grid;

plot(w/pi,mag);

xlabel('f');ylabel('s(t)');

axis([0,1,0,1.5]);

y1=filter(b,a,st);

subplot(4,3,2);grid;

plot(t,y1);grid;xlabel('t/s');ylabel('s(t)');

axis([0,Tp,min(y1),max(y1)]);

fxt=fft(y1,N);         %计算信号y1的频谱

subplot(4,3,3);

stem(f,abs(fxt)/max(abs(fxt)),'.');grid;

axis([0,Fs/8,0,1.2]);

xlabel('f/Hz');ylabel('幅度');

带通

Wp=[0.09*pi,0.13*pi];

Ws=[0.05*pi,0.2*pi];

[E,wn]=buttord(Wp/pi,Ws/pi,Ap,As);              %计算巴特沃斯滤波器阶次和截止频率

[b,a]=butter(E,wn,'bandpass');                     %频率变换法设计巴特沃斯带通滤波器

[db,mag,pha,grd,w]=freqz_m(b,a);                  %数字滤波器响应

y2= filter(b,a,st);

subplot(4,3,4);

plot(w/pi,mag);

xlabel('f');ylabel('s(t)');

axis([0,1,0,1.5]);

subplot(4,3,5);

plot(t,y2);grid;xlabel('t/s');ylabel('s(t)');

axis([0,Tp,min(st),max(st)]);title('(a)s(t)的波形');

fxt=fft(y2,N);         %计算信号st的频谱

subplot(4,3,6);

stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b)s(t)的频谱');

axis([0,Fs/8,0,1.2]);

xlabel('f/Hz');ylabel('幅度');

高通

Wp=0.15*pi;

Ws=0.1*pi;

[F,wn]=buttord(Wp/pi,Ws/pi,Ap,As);              %计算巴特沃斯滤波器阶次和截止频率

[d,c]=butter(F,wn,'high');                         %频率变换法设计巴特沃斯高通滤波器

y3=filter(d,c,st);

subplot(4,3,7);

[db,mag,pha,grd,w]=freqz_m(d,c);                  %数字滤波器响应

plot(w/pi,mag);

subplot(4,3,8);

plot(t,y3);grid;xlabel('t/s');ylabel('s(t)');

axis([0,Tp,min(st),max(st)]);title('(a)s(t)的波形');

fxt=fft(y3,N);         %计算信号st的频谱

subplot(4,3,9);

stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b)s(t)的频谱');

axis([0,Fs/8,0,1.2]);

xlabel('f/Hz');ylabel('幅度');

<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" /> 

第三节

[y,fs,nbits]=wavread('F:\a.wav');

%sound(y,fs,nbits);

N=length(y);

Y=fft(y,N);

subplot(3,2,1);

plot(y);

subplot(3,2,2);

plot(abs(Y));

t=0.02*randn(N,1);

x=y+t;

%sound(x,fs,nbits);

subplot(3,2,3);

plot(x);

S=fft(x);

subplot(3,2,4);

plot(abs(S));

Wp=0.25*pi;

Ws=0.3*pi;

Ap=0.1;

As=40;

[N,wn]=buttord(Wp/pi,Ws/pi,Ap,As);              %计算巴特沃斯滤波器阶次和截止频率

[b,a]=butter(N,wn,'low');                         %频率变换法设计巴特沃斯高通滤波器

So=filter(b,a,x);

subplot(3,2,5);

plot(x);

sf=fft(So);

subplot(3,2,6);

plot(abs(sf));

sound(So,fs,nbits);


提 示:本站所有源码只提供浏览,需要定制的朋友可以联系在线客服!
注意:价格根据项目的难易程度来定价格业务咨询

matlab信号分析 --下载地址

注意:为避免不必要的误会,本站项目里的源码只公开部分需要的联系在线客服

本站长期招聘程序代写高手,欢迎加入华南地区matlab团队

想创业却没有经验的人

无论你是否有过网上开店的经验,都可以随时联系在线客服,建立自己独立的网站
想开网店却不知道如何入手

淘宝创业成本低而且风险小,如果想开淘宝店的朋友可以联系在线客服。
想兼职创业,却不擅长交际与服务的人

在家创业月入5000元。网站程序+百套群发工具+网赚资料+域名+空间+本站终身代理资格,这样你网赚的条件全具备了。每天3小时管理、推广、收钱。
缺乏能快速赢利型产品的人

导入多种最新流行营销软件+网赚教程,让入驻者轻松加盟、复制有效成交技巧、快速赚钱。

源码评论评论内容只代表网友观点,与本站立场无关!

   评论摘要(共 0 条,得分 0 分,平均 0 分) 查看完整评论

浏览说明

* 本站所有源码全部公开,随时随地浏览!
* MATLAB软件如用于商业用途,请购买正版!
* 如果您发现下载链接错误,请点击报告错误谢谢!
* 站内提供的所有软件包含破解及注册码均是由网上搜集,若侵犯了你的版权利益,敬请来信通知我们!
Copyright © 2008-2014 www.buycode.cn. All Rights Reserved.
页面执行时间:218.75000 毫秒