MATLAB有关文献报告参考模板.docx
- 文档编号:15086053
- 上传时间:2023-06-30
- 格式:DOCX
- 页数:13
- 大小:183.64KB
MATLAB有关文献报告参考模板.docx
《MATLAB有关文献报告参考模板.docx》由会员分享,可在线阅读,更多相关《MATLAB有关文献报告参考模板.docx(13页珍藏版)》请在冰点文库上搜索。
MATLAB有关文献报告参考模板
论文学习报告
——对论文《MatlabSimulationinSignals&Systems》学习
摘要:
主要通过对学者的论文的学习和论文中仿真程序的运行和学习,掌握一些基本写论文的步骤和要点,同时加深对matlab的学习并对matlab在信号与系统中的应用有了基本的了解。
关键字:
matlab、信号与系统、应用、论文学习
一、介绍
MATLAB具有丰富和强大的功能,可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
尽管Matlab提供了许多主要信号和系统处理功能,如卷积,傅里叶转换、拉普拉斯变换,z变换等,它简化计算过程大大,同时一些相关的数学处理和物理意义被掩盖了。
我学习参考的文献《MatlabSimulationinSignals&Systems》主要通过Matlab在三个不同的水平实现和应用,运用函数模拟达到使学生在学习过程中对理论更深入全面的理解。
我对这篇文献的学习主要是
(1)了解论文的基本结构和框架;
(2)对文献内容进一步了解学习;(3)运用MATLAB仿真实现文献的内容并对出现的问题分析解决;(4)总结学习经验用论文的格式完成报告。
二、对论文的学习
论文一般由题名、作者、摘要、关键词、正文、参考文献和附录等部分组成,其中部分组成(例如附录)可有可无。
论文各组成的排序为:
题名、作者、摘要、关键词、英文题名、英文摘要、英文关键词、正文、参考文献、附录和致谢。
论文《MatlabSimulationinSignals&Systems》把MATLAB的应用分为三个阶段:
(1)第一阶段:
直接应用MATLAB强大的功能(FirstLevel-UsingMatlab’sPowerfulFunctionsDirectly);
(2)第二阶段:
运用MATLAB实现数学过程(SecondLevel—UsingMatlabtoRealizetheMathematicalProcess);(3)第一阶段:
开发新的优化算法(ThirdLevel—DevelopingnewOptimizingalgorithms);并通过介绍三个层面上应用MATLAB实现傅里叶变换,举例说明运用函数模拟达到使学生在学习过程中对理论更深入全面的理解。
三、文献的MATLAB仿真实现
——(程序中的问题直接在仿真过程中解决)
第一阶段、
symst;%Defineasignal.
F=fourier(heaviside(t))%CalculatetheFouriertransformofstepfunction.
ezplot(abs(F1),[-3,3]);%Drawmagnitude-frequencycurve.
F=
pi*dirac(w)-i/w
Undefinedfunctionorvariable'F1'.
改’F1’为‘F’结果正确且无误,图如下
第二阶段、
(1)
w=linspace(-6*pi,6*pi,512);%Assignthefrequencypointstobecalculated.
N=length(w);%Getthetotalnumberofthefrequencypoint.
X=zeros(1,N);%SetavectorXoflengthNtosavethecalculationresults.
fork=1:
N
X(k)=quad(@(t)qpulse(t,w(k)),-1,1);
end
%Foreachfrequencypointw(k),calculatethedefiniteintegral
.‘qpulse’isthe
squarewavefunctiondefinedinthefollowing
qpulse.mfile.
plot(w,real(X));gridon;%DrawFouriertransform.
The‘qpulse’functionquotedisdefinedasthefollowing.
functiony=qpulse(t,w);
y=(stepfun(t,-1)-stepfun(t,1)).*exp(-j*w*t);
(2)
T=2;%Definethetimedomainsampleintervallength.
N=200;%Definethenumberofthesamplepoints.
t=linspace(-T/2,T/2-T/N,N)';%Definethesamplepoints.
f=0*t;%Initializethetimefunction.
f(t>-1&t<1)=1;%Assignvaluetotimefunction.
OMG=16*pi;%Definethefrequencydomainsample
intervallength.
K=100;%Definethenumberofthesamplepoints.
w=linspace(-OMG/2,OMG/2-OMG/K,K)';
%Definethesamplepoints.
F=0*w;%Initializethespectrumfunction.
fork=1:
K%Calculatethespectrumforeachfrequency
pointinloop.
forn=1:
N
F(k)=F(k)+T/N*f(n)*exp(-j*w(k)*t(n));
end%Finishthesumin
(2)inloop.
end
k=1:
K;
plot(k,real(F));
用
fork=1:
K
F(k)=T/N*exp(-j*w(k)*t).'*f;
End
取代
fork=1:
K%Calculatethespectrumforeachfrequency
pointinloop.
forn=1:
N
F(k)=F(k)+T/N*f(n)*exp(-j*w(k)*t(n));
end%Finishthesumin
(2)inloop.
end
用
U=exp(-j*kron(w,t.’));%Calculatethetensorproductof
twovectors.
F=T/N*U*f;%Fouriertransformisfinished.
代替
fork=1:
K%Calculatethespectrumforeachfrequencypointinloop.
forn=1:
N
F(k)=F(k)+T/N*f(n)*exp(-j*w(k)*t(n));
end%Finishthesumin
(2)inloop.
end
k=1:
K;
(3)
function[w,tw]=cconv1(u,tu,v,tv)
%uandvaretwovectorstobeconvoluted,andtuandtv
aretheirsamplingtimesrespectively.
%wistheconvolutionresultandtwisw’ssamplingtime.
T=tu
(2)-tu
(1);%Getsampleinterval.
w=T*conv(u,v);%Calculate(6)viafunction‘conv’.
tw=tu
(1)+tv
(1)+T*[0:
length(u)+length(v)-2]';
%Distributesampletimeforw.
主程序:
t=[-2:
0.01:
2];%Generatesamplingtime.
e=(t>-1&t<1);%Defineasquarewavewithwidthof2.
h=(t>-0.5&t<0.5);
%Defineasquarewavewithwidthof1.
[r1,t1]=cconv1(e,t,h,t);
%Call‘cconv1’tofinishconvolutionintegral.
subplot(3,1,1);plot(t,e);axis([-22-0.21.2]);title('e(t)');
subplot(3,1,2);plot(t,h);axis([-22-0.21.2]);title('h(t)');
subplot(3,1,3);plot(t1,r1);axis([-44-0.21.2]);title('r(t)');
运行出现以下错误:
Undefinedfunction'conv2'forinputargumentsoftype'char'.
Errorinconv(line39)
c=conv2(a(:
),b(:
),shape);
Errorincconv1(line6)
w=T*conv(u,v);%Calculate(6)viafunction‘conv’.
修改:
创建一个新的函数myconv
functiona=myconv(b,c)
bs=size(b);
cs=size(c);
i=any(bs-cs);
ifi
error('error')
end
i=any(~(bs-1));
if~i
error('error')
end
ko=0;
ifbs
(1)>bs
(2)
b=b';
c=c';
ko=1;
end
bs=size(b);
cs=size(c);
ss=2*bs
(2)-1;
a=zeros(1,ss);
fori=1:
cs
(2)
q=zeros(1,i-1);
p=zeros(1,ss-cs
(2)+1-i);
ba=[q,c,p];
ma=b(i)*ba;
a=a+ma;
end
ifko
a=a';
end
end
cconv1中调用conv也改用mycovn得到以下结果:
说明运行无误。
function[w,tw]=cconv2(u,tu,v,tv)
T=tu
(2)-tu
(1);
L1=length(u);%Getthelengthofvectoru.
L2=length(v);%Getthelengthofvectorv.
forn=1:
L1+L2-1%Definethelengthofw.
w(n)=0;%Initializew.
form=max(1,n+1-L2):
mi(n,L1);
%Determinethestartingandendpointsofthesumming.
w(n)=w(n)+u(m)*v(n-m+1);%Calculatethesum.
end
w(n)=T*w(n);%Modifytheconvolutionsumto
convolutionintegral.
end
tw=tu
(1)+tv
(1)+T*[0:
L1+L2-2]';
用ccovn2代替ccovn1也得到以下结果:
至此文献中程序全部仿真完成。
参考文献:
(1)Xiao-hanGuan,Meng-mengZhang,YongZheng,《MatlabSimulationinSignals&Systems——UsingMatlabatdifferentlevels》,CollegeofElectromechanicalEngineeringNorthChinaUniversityofTechnology,Beijing,China,2009
(2)MATLAB数值计算方法\张德丰等著.—北京:
机械工业出版社,2010.1
致谢:
(1)王刚老师关于《MATLAB数值计算方法》的教授;
(2)王,蒋等同学的讨论;
(3)电子科技大学图书馆,XX知道,XX文库的大力支持
友情提示:
范文可能无法思考和涵盖全面,供参考!
最好找专业人士起草或审核后使用,感谢您的下载!
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 有关 文献 报告 参考 模板