服务承诺





51Due提供Essay,Paper,Report,Assignment等学科作业的代写与辅导,同时涵盖Personal Statement,转学申请等留学文书代写。




私人订制你的未来职场 世界名企,高端行业岗位等 在新的起点上实现更高水平的发展




matlab信号分析命令代码解析--matlab编程代写
2017-02-22 来源: 51Due教员组 类别: 写作技巧
matlab是美国MathWorks公司旗下的一款基于模型的开发设计软件,也是让无数留学生们了为之头疼的一个软件,当国外大学教授们在布置matlab编程作业的时候,想必各位同学们的内心也开始走向崩溃,这也说明想要完成matlab作业并不是一件难事,我们需要的是学会其中的技巧,今天matlab编程代写平台要来为各位介绍的就是在matlab中那些常见的信号分析命令代码,让我们一起来看看吧。
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('幅度');
第三节
[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最常见的信号分析命令之后,各位留学生们日后再用matlab完成作业时,应该也已经学会该如何使用matlab信号分析了吧,当然,如果各位还是觉得matlab太难是实在难以完成,那么你可以找我们来帮忙,我们将为您提供最为优秀的matlab代写服务。
51due教育引领海外留学的风向标,打造最好的留学教育品牌,我们为自己的品牌代言,不负留学生的期待。提供全方位的服务,包括留学生作业代写、统计作业代写、essay代写、assignment代写服务。
更多关于写作技巧的内容,点开主页栏目 即可快速获取,有论文代写需求,咨询24小时客服QQ:800020041,为您详细解答。-xz
