南邮Matlab实验报告三

Matlab上机实验报告

实验名称：MATLAB的数值运算

班级：自动化二班

学号： B11050216

姓名：李鹏飞

南京邮电大学

2013年4月9日

一、实验目的

1．学习MATLAB 的基本矩阵运算； 2．学习MATLAB 的点运算；

3．学习复杂运算。

二、实验基本知识

1.基本矩阵运算；

2.多项式运算；

3.代数方程组；

4.数值微积分。

三、上机练习

1. 给a,b,c 赋如下数据:

]6,46,23,4,2,6,3,8,0,1[,356838241248

7,278744125431-=??????????--=??????????=c b a

1)求a+b,a*b,a.*b,a/b,a./b,a^2,a.^2,length(c)的结果.

2)求c 中所有元素的平均值、最大值.

3)求d=b(2:3,[1,3]).

程序代码：

（1）

>> a=[1,3,4;5,12,44;7,8,27]; >> b=[-7,8,4;12,24,38;68,-5,3]; >> c=[1,0,8,3,6,2,-4,23,46,6]; >> d=a+b d =

-6 11 8 17 36 82 75 3 30

>> e=a*b

e =

301 60 130

3101 108 608

1883 113 413 >> f=a.*b

f =

-7 24 16

60 288 1672

476 -40 81 >> g=a/b

g =

0.0966 0.0945 0.0080

-3.6125 1.5838 -0.5778

-1.9917 0.9414 -0.2682

>> h=a./b

h =

-0.1429 0.3750 1.0000

0.4167 0.5000 1.1579

0.1029 -1.6000 9.0000

>> i=a^2

i =

44 71 244

373 511 1736

236 333 1109 >> j=a.^2

j =

1 9 16 25 144 1936 49 64 729 >> k=length(c) k =

10

（2）

>> l=max(c) l =

46

>> m=mean(c) m =

9.1000

（3）

>> d=b(2:3,[1,3]) d =

12 38

68 3

2.求解方程???

?

??????=??????????001987654321X ，并对所得结果作出解释。

程序代码：

>> A=[1,2,3;4,5,6;7,8,9];

b=[1,0,-1]'; ra=rank(A) ra =

2

>> rab=rank([A,b])

rab =

2

>> xs=A\b;

Warning: Matrix is close to singular or badly scaled.

Results may be inaccurate. RCOND = 1.541976e-018. >> xg=null(A); c=rand(1);

ba=A*(xs+c*xg) ba =

1.0000 -0.0000 -1.0000

>> norm(ba-b)

ans =

1.7764e-015

3.已知有理分式)

()

()(x D x N x R =

，其中)5.0)(3()(33++=x x x x N ，

)125)(22()(232++-+=x x x x x D 。求该分式的商多项式和余多项式。

程序代码：

>> format rat

p1=conv([3,0,1,0],[1,0,0,0.5]); p2=[5 12 -6 -3 2 -2]; [q,r]=deconv(p1,p2);

cq='商多项式为 ';cr='余多项式为 ';

disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')]) 商多项式为 0.6 s - 1.44

余多项式为 21.88 s^4 - 5.34 s^3 - 5.52 s^2 + 4.58 s - 2.88 >> qp2=conv(q,p2); pp1=qp2+r; pp1==p1

ans =

1 1 0 0 1 1 1 4.求解Lorenz模型的状态方程，并图示.

初值在(0,1)中取值

程序代码：

M-文件：

function dx = myfun(t,x)

dx = zeros(3,1);

dx(1)=16*(x(2)-x(1));

dx(2)=45*x(1)-x(2)-x(1)*x(3);

dx(3)=x(1)*x(2)-4*x(3) ;

comand window中：

[T,Y] = ode45(@myfun,[0 1],[1 1 1]); plot(T,Y)

tspan=[0,1];

x0=[0.5;0.4;0.6];

[tt,xx]=ode45(@myfun,tspan,x0); subplot(3,1,1),

plot(tt,xx(:,1),'-b')

hold on

xlabel('t')

ylabel('x1(t)')

title('Lorenz?￡Díí?')

legend('\fontsize{14}\it x1(t)') hold on

subplot(3,1,2),

plot(tt,xx(:,2),'r*')

xlabel('t')

ylabel('x2(t)')

legend('\fontsize{14}\it x2(t)') subplot(3,1,3),

plot(tt,xx(:,3),'go')

xlabel('t')

ylabel('x3(t)')

legend('\fontsize{14}\it x3(t)')

0.1

0.2

0.3

0.4

0.50.6

0.7

0.8

0.9

1

-50050

t

x 1(t )

Lorenz 模型图

0.1

0.2

0.3

0.4

0.50.6

0.7

0.8

0.9

1

t

x 2(t )

t

x 3(t )

四、实验小结：

在实验三：MATLAB 的符号运算中，我们学会了用MATLAB 解决实际的数学和符号运算问题，但是在使用软件的过程中，我还是遇到了不少的问题，有很多的代码不熟悉，这需要我在日后的学习过程中不断加强深入学习，多增加动手操作，争取早日熟悉MATLAB 的基本操作要求。在实验第四题中，用到了M 文件的使用，同时我们也学习到了数值计算里的ode45指令，并熟练使用M 文件。

总之在这次MATLAB 实验中我们学到了解决数值函数问题的方法，初步掌握了数值运算，为日后再继续深入学习MATLAB 打下了基础。

B11050216

李鹏飞

《MATLAB与数值分析》第一次上机实验报告

电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名：李培睿 学号：2013020904026 指导教师：程建

一、实验名称 《MATLAB与数值分析》第一次上机实验 二、实验目的 1. 熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算 操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序） 2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号 转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序） 3. 掌握Matlab函数的编写规范。 4、掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、 三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释） 5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。 三、实验内容 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x， y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 2. 编程实现奥运5环图，允许用户输入环的直径。 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort 函数。 四、实验数据及结果分析 题目一： ①在Editor窗口编写函数代码如下：

MATLAB实验报告50059

实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1

例1-2 例1-3 例1-4

3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。

5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？

南京邮电大学

南京邮电大学 毕业设计(论文)任务书题目：OFDMA系统信道估计的研究 专业通信工程 学生姓名施祥威 班级学号10001417 指导教师(签字) 指导教师职称副教授 指导单位通信与信息学院 负责人签字 院(系)领导签字 日期2013.12

题目OFDMA系统信道估计的研究 课题类型理论研究■工程设计□ 产品开发□实验研究□ 课题类别 软件■硬件□ 软硬结合□其他□ 设计 内容与技术要求、成果形式 研究和学习无线移动通信原理，探讨数字通信系统的信号解调接收，针对信道对信号造成的失真，要进行信道估计，然后进行均衡（补偿）。具体到OFDMA系统，要力求准确估算多频率、时间点信道特性。 要求探讨无线信道的特点和模型，参考WiMAX系统标准，建立OFDMA系统模型。分析上、下行链路信道特点，分析不同导频符号安排下的信道估计算法。分析多种不同的信道估计算法，并作比较。在解调中用信道估计值进行补偿。 用Matlab仿真，绘出 ~,~ b MSE SNR BER E n曲线，分析其性能。提交程序清单，仿真数据。 设计进度11月中旬～12月初：查阅资料，学习无线通信和OFDMA的基础知识，消化课题内容，掌握仿真工具的使用方法，拟定系统总体实现方案。写开题报告。3月1日~3月底：依照总体方案，逐个编写程序，调试。联接程序，总体调试。 4月1日~5月底：调试程序，修改完善，测试系统数据、曲线，分析性能。撰写毕业论文。 6月1日~6月20日：完成毕业论文，毕业答辩。 参考资料[1]IEEE 802.16 IMT-Advanced Evaluation Group Coordination Meeting, Overview of IEEE P802.16m Technoledge and Candidate Rit for IMT-Advanced[S], IEEE L802.16-10/0002, IEEE,2010 [2](芬) 霍尔马 (Holma.H.) . (芬) 托斯卡拉 (Toskala.A.), UMTS中的LTE-基于OFDMA和SC-FDMA的无线接入, [M]机械工业出版社，2010 [3]Ramjee Prasad l Fernando J. Velez，WiMAX Networks: Techno-Economic Vision and Challenges[M], Springer press, 2010 [4]W-C.Huang,C-H Pan, C-P Li,H-J Li, Subspace-Based Semi-Blind Channel Estimation in Uplink OFDMA Systems[J], IEEE Trans. on Broadcasting, vol. 56, no. 1, March 2010 [5]沈嘉，下一代无线宽带通信的核心OFDM/OFDMA[J]，计算机世界，2007.1 [6]Yong Soo Cho , Jaekwon Kim, Won Young Yang , Chung G.Kang, MIMO-OFDM无线通信技术及MATLAB实现[M]，电子工业出版社，2013.6 [7]徐明远、邵玉斌，MATLAB仿真在通信与电子工程中的应用(第2版)[M]，西安电子科技大学出版社，2010

Matlab实验报告3

实验三函数的可视化与Matlab作图 一、按要求绘制如下曲线（面）： 1. 在[0，4pi]上画sin(x),cos(x)在同一图像中，其中cos(x)图像用红色小圆圈，并在函数图上标注“y=sin(x)”，”y=cos(x)”,X轴，Y轴，标题为“正弦余弦函数图像。”答：>> clear >> clf, x=linspace(0,4*pi,200);y1=sin(x);y2=cos(x); plot(x,y1,'k-',x,y2,'ro') >> title('正弦余弦函数图像。') >> legend('y=sin(x)','y=cos(x)') >> ylabel('\it{Y轴}'); >> xlabel('\it{X轴}'); 2.任意绘制彗星曲线图。 答：>> clf; >> x=[1:10]; y=[5 6 3 4 8 1 10 3 5 6]; >> z=0:0.1:100; x=sin(z);y=cos(z).*10; >> %三维彗星图 comet3(x,y,z) >> %二维彗星图

t = -pi:pi/200:pi; comet(t,tan(sin(t))-sin(tan(t)))

3.在多窗口中绘制y=sin(t)*sin(t);y1=sin(3*t+2.5);y2=sin(5*t+5)并加以标注。答：>> clf; t=0:0.1:4*pi; subplot(3,1,1),plot(sin(t).*sin(t)),legend('y=sin(t)*sin(t)') subplot(3,1,2),plot(sin(3*t+2.5)),legend('y1=sin(3*t+2.5)') subplot(3,1,3),plot(sin(5*t+5)),legend('y2=sin(5*t+5)') 4.自拟题目绘制三维线图。 绘制以下方程y1=sin(t),y2=cos(t),x=t在t=[0,2π] 对应的三维曲线。 >> clf; >> t=0:pi/10:2*pi; >> y1=sin(t);y2=cos(t); >> plot3(y1,y2,t);grid on; >> xlabel('Dependent Variable Y1'); >> ylabel('Dependent Variable Y2'); >> zlabel('Dependent Variable X'); >> title('Sin and Cos Curve');

南邮MATLAB软件设计

通信与信息工程学院 2015 / 2016 学年第 2 学期 软件设计实验报告 模块名称MATLAB软件设计 专业通信工程 学生班级B130102 学生学号B13010229 学生姓名郭智超 指导教师孟田甜

报告内容 一、要求练习的实验部分 1.题目：在时间区间 [0,10]中，绘出 t e y t 2cos 15.0--=曲线。 程序代码： t=0:0.02:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'b'); xlabel('time'); ylabel('amplitude'); title('figure 1'); grid 实验结果： 24 6810 00.51 1.5 time a m p l i t u d e figure 1 2.题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦波，要求它的负半波被置零，且在 2 3 处被削顶。

程序代码： t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点） y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(ytol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ? （3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ?

matlab实验报告3详解

实验四、LTI系统的响应 课程名称： MATLAB应用技术专业班级：通信1422 学生学号： 1430119231 学生姓名：周妍智 所属院部：电子信息工程系指导教师：徐树梅 2015 —— 2016 学年第二学期

实验项目名称： LTI 系统的响应 实验学时： 16 学生姓名： 周妍智 实验地点： 微机11 实验日期： 2016.4.17 实验成绩： 批改教师： 徐树梅 批改时间： 一、 实验目的 1. 熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法 2. 熟悉连续（离散）时间系统在任意信号激励下响应的求解方法 3. 熟悉应用MATLAB 实现求解系统响应的方法 二、 实验原理 1.连续时间系统 对于连续的LTI 系统，当系统输入为f (t )，输出为y (t )，则输入与输出之间满足如下的线性常系数微分方程： () ()0 ()()n m i j i j i j a y t b f t ===∑∑，当系统输入为单位冲激信号δ(t )时产生 的零状态响应称为系统的单位冲激响应，用h(t)表示。若输入为单位阶跃信号ε(t )时，系统产生的零状态响应则称为系统的单位阶跃响应，记为g(t)，如下图所示。 系统的单位冲激响应h (t )包含了系统的固有特性，它是由系统本身的结构及参数所决定的，与系统的输入无关。我们只要知道了系统的冲激响应，即可求得系统在不同激励下产生的响应。因此，求解系统的冲激响应h(t )对我们进行连续系统的分析具有非常重要的意义。 在MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应， 并绘制其时域波形的函数impulse( ) 和step( )。如果系统输入为f (t )，冲激响应为h(t)，系统的零状态响应为y (t )，则有：()()()y t h t f t =*。 若已知系统的输入信号及初始状态，我们便可以用微分方程的经典时域求解方法，求出系统的响应。但是对于高阶系统，手工计算这一问题的过程非常困难和繁琐。 在MATLAB 中，应用lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真，求出系统在任意激励信号作用下的响应。lsim( )函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解，而且还能同时绘制出系统响应的时域波形图。 以上各函数的调用格式如下： ⑴ impulse( ) 函数 函数impulse( )将绘制出由向量a 和b 所表示的连续系统在指定时间范围内的单位冲激响应h (t )的时域波形图，并能求出指定时间范围内冲激响应的数值解。

matlab实验报告

实验一小球做自由落体运动内容：一小球竖直方向做自由落体，并无损做往返运动。程序： theta=0:0.01:2*pi x=cos(theta) y=sin(theta) l=1 v=1 while l<10 for t=1:10 y=y+(-1)^l*v*t plot(x,y,[-1,1],[-56,2],'.') axis equal pause(0.1) end l=l+1 end 结果：

-50 -40 -30 -20 -10 收获：通过运用小球自由落体规律，及（-1）^n 来实现无损往 返运动！ 实验二 旋转五角星 内容：一个五角星在圆内匀速旋转 程序：x=[2 2 2 2 2 2] y=[0 4/5*pi 8/5*pi 2/5*pi 6/5*pi 0] y1=2*sin(y) x1=2*cos(y) theta=0:4/5*pi:4*pi

x2=2*cos(theta) y2=2*sin(theta) plot(x,y,x1,y1,x2,y2) axis equal theta1=theta+pi/10 x2=2*cos(theta1) y2=2*sin(theta1) plot(x2,y2) axis equal theta=0:4/5*pi:4*pi for rot=pi/10:pi/10:2*pi x=2*cos(theta+rot) y=2*sin(theta+rot) plot(x,y) pause(0.1) end 结果：

-2 -1.5-1-0.500.51 1.52 -2-1.5-1-0.500.511.5 2 收获：通过theta1=theta+pi/10，我们可以实现五角星在圆内匀速 旋转！ 实验三 转动的自行车 内容：一辆自行车在圆内匀速转动 程序：x=-4:0.08:4; y=sqrt(16-x.^2); theta1=-pi/2:0.01*pi:3*pi/2; x3=0.5*cos(theta1); y3=0.5*sin(theta1); theta=-pi/2+0.02*pi for k=1:100

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;

>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =

0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1

6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5

8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =

参考答案Matlab实验报告

实验一 Matlab基础知识 一、实验目的： 1.熟悉启动和退出Matlab的方法。 2.熟悉Matlab命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握Matlab各种表达式的书写规则以及常用函数的使 用。 二、实验内容： 1.求[100，999]之间能被21整除的数的个数。（rem） 2.建立一个字符串向量，删除其中的大写字母。（find） 3.输入矩阵，并找出其中大于或等于5的元素。（find） 4.不采用循环的形式求出和式 63 1 2i i= ∑ 的数值解。(sum) 三、实验步骤： ●求[100，199]之间能被21整除的数的个数。（rem） 1.开始→程序→Matlab 2.输入命令： ?m=100:999; ?p=rem(m,21); ?q=sum(p==0) ans=43 ●建立一个字符串向量，删除其中的大写字母。（find） 1.输入命令：

?k=input('’,’s’)； Eie48458DHUEI4778 ?f=find(k>=’A’&k<=’Z’)； f=9 10 11 12 13 ?k(f)=[ ] K=eie484584778 ●输入矩阵，并找出其中大于或等于5的元素。（find） 1.输入命令： ?h=[4 8 10;3 6 9; 5 7 3]； ?[i,j]=find(h>=5) i=3 j=1 1 2 2 2 3 2 1 3 2 3 ●不采用循环的形式求出和式的数值解。（sum） 1.输入命令： ?w=1:63； ?q=sum(2.^w) q=1.8447e+019

实验二 Matlab 基本程序 一、 实验目的： 1. 熟悉Matlab 的环境与工作空间。 2. 熟悉M 文件与M 函数的编写与应用。 3. 熟悉Matlab 的控制语句。 4. 掌握if,switch,for 等语句的使用。 二、 实验内容： 1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 2. 编程完成，对输入的函数的百分制成绩进行等绩转换，90~100为优，80~89为良，70~79为中，60~69为及格。 3. 编写M 函数文件表示函数 ，并分别求x=12和56时的函数值。 4. 编程求分段函数 2226;03 56;0532 1;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它，并求输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。 三、 实验步骤： 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 1. 打开Matlab ，新建M 文件 2. 输入命令： 51022-+x

南京邮电大学 数值代数实验

数值代数实验 数值线性代数实验一 一、实验名称：矩阵的LU分解. 二、实验目的：用不选主元的LU分解和列主元LU分解求解线性方程组Ax=b, 并比较这 两种方法. 三、实验内容与要求 （1）用所熟悉的计算机语言将不选主元和列主元LU分解编成通用的子程序，然后用编写的程序求解下面的84阶方程组 将计算结果与方程组的精确解进行比较，并就此谈谈你对Gauss消去法的看法. （2）写出追赶法求解三对角方程组的过程，并编写程序求该实验中的方程组 Gauss消去法： 用消去法解方程组的基本思想是用逐次消去未知数的方法把原来方程组Ax=b化为与其等价的三角方程组，而求解三角方程组就容易了。换句话说，上述过程就是用行的初等变换将原方程组系数矩阵化为简单形式，从而将求解原方程组的问题转化为求解简单方程组的问题。 利用Gauss消去法对线性方程组Ax=b进行求解。 用MATLAB建立m文件DelGauss.m,程序如下： function x=DelGauss(a,b) [n,m]=size(a); nb=length(b); det=1; x=zeros(n,1); for k=1:n-1 for i=k+1:n if a(k,k)==0 return end

m=a(i,k)/a(k,k); for j=k+1:n a(i,j)=a(i,j)-m*a(k,j); end b(i)=b(i)-m*b(k); end det=det*a(k,k); end det=det*a(n,n); for k=n:-1:1 for j=k+1:n b(k)=b(k)-a(k,j)*x(j); end x(k)=b(k)/a(k,k); End 在matlab中输入如下： 结果如下：

MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案

“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)

1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408

3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans =

matlab实验报告

Matlab实验报告 实验二图像处理 一、实验目的 （1）通过应用MA TLAB语言编程实现对图像的处理，进一步熟悉MATLAB软件的编程及应用； （2）通过实验进一步掌握图像处理的基本技术和方法。 二、实验内容及代码 ㈠.应用MA TLAB语言编写显示一幅灰度图像、二值图像、索引图像及彩色图像的程序，并进行相互之间的转换 首先，在matlab页面中的current directory下打开存放图像的文件夹。 1.显示各种图像 ⑴显示彩色图像： ①代码：>> mousetif=imread('tif.TIF'); >> image(mousetif) 显示截图： ②代码：>> mousetif=imread('tif.TIF'); >> imshow(mousetif) 显示截图：

③代码：mousetif=imread('tif.TIF'); subimage(mousetif) 显示截图： 显示截图：

⑵显示二值图像 ①代码：>> I=imread('单色bmp.bmp'); >> imagesc(I,[0 2]) 显示截图： ②代码：>> I=imread('单色bmp.bmp');

>> imshow(I,2) 显示截图： ③代码：>> I=imread('单色bmp.bmp'); >> subimage(I) 显示截图：

⑶显示灰度图像 ①代码：>> I1=imread('256bmp.bmp'); >> imagesc(I1,[0,256]) 显示截图： 代码：>> I1=imread('256bmp.bmp'); >> colormap(gray); >> subplot(1,2,1); >> imagesc(I1,[0,256]); >> title('灰度级为[0 256]的mouse.bmp图'); >> subplot(1,2,2); >> imagesc(I1,[0,64]); >> colormap(gray); >> title('灰度级为[0 64]的mouse.bmp图'); 显示截图：

南京邮电大学matlab软件设计(超详细).

邮电大学通信学院 软件课程设计 实验报告 模块名称：___MATLAB软件设计 专业班级：__通信工程 姓名：____ ____ 学号：___ 实验日期：2013年 6 月 17—28日实验报告日期： 2013 年 7 月 1 日

一、要求练习的实验部分 1．在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。 程序： t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 结果： 2. 写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦波，要求它的负半 波被置零，且在2 3处被削顶。 程序： t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(y

xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数') legend('y=sin(t)','z=f(t)',4),hold off 结果： 3. 令???? ??????=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ? （3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ? （5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ? （7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ? 程序： clear; clc; disp('设 A=[1,2,3;4,5,6;7,8,9]，得')

Matlab程序设计实验报告

实验七Matlab 程序设计 实验目的： 1、掌握建立和执行M 文件的方法； 2、掌握实现选择结构的方法； 3、掌握实现循环结构的方法。 实验内容： 1. 编写用 5 次多项式拟合函数y=sin(x), x [0, 2 ]的脚本M 文件，要求绘图观察拟合的效果。 function shiyan1 x=0:0.5:2*pi y=sin(x) p=polyfit(x,y,5) x1=0:0.2:2*pi y1=polyval(p,x1) plot(x,y, 'b' ,x1,y1, '*r' x =

Columns 1 through 9 0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 Columns 10 through 13 4.5000 5.0000 5.5000 6.0000 y = Columns 1 through 9 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411 -0.3508 -0.7568 Columns 10 through 13 -0.9775 -0.9589 -0.7055 -0.2794 p = -0.0056 0.0881 -0.3967 0.2671 0.8902 0.0029 x1 = Columns 1 through 10 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 Columns 11 through 20

2. 2.2000 2.4000 2.6000 2.8000 3.0000 3.2000 3.4000 3.6000 1.8001 Columns 21 through 30 4.0 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 5.6000 5.8000 Columns 31 through 32 6.0 6.2000 y1 = Columns 1 through 10 0.29 0.1886 0.3786 0.5585 0.7172 0.8461 0.9391 0.9926 1.0048 0.9761 Columns 11 through 20 0.9083 0.8048 0.6701 0.5098 0.3301 0.1381 -0.0590 -0.2538 -0.4389 -0.6073 Columns 21 through 30 -0.7524 -0.8685 -0.9505 -0.9949 -0.9991 -0.9626 -0.8863 -0.7732 -0.6288 -0.4606 Columns 31 through 32

广州大学学生实验报告1 matlab 程序设计

广州大学学生实验报告 开课学院及实验室：机械与电气工程学院计算机楼 301室2014 年10 月30 日

2、MATLAB指令窗的基本操作 MATLAB指令窗给用户提供了最直接的交互界面，可用于输入和执行指令、显示指令运行结果、调试MATLAB程序等常用的MATLAB仿真计算功能。本实验掌握以下在指令窗执行的基本操作，达到熟悉使用指令窗的目的： （1）最简单的计算器使用方法：在MATLAB指令窗中，可按计算器的方式进行一般的数学计算，MATLAB的运算符的含义大致与常见的运算规则一致； （2）在指令窗中输入和生成矩阵：与一般的计算器不同，在MATLAB中可直接输入和生成矩阵。实际上，矩阵是MATLAB工作的基本元素。 （3）数值表述方法：在MATLAB中的大部分数值的表述方式与平常是相同的，需要注意的是在表示比较大的数时，MATLAB默认采用科学计数法显示； （4）变量命名规则：对于MATLAB变量命名规则，需要注意以下几点： a、变量名、函数名对字母大小写敏感 b、变量名的第一个字母必须是英文字母，后续可以是字母、数字、下划线 c、变量的有效时限：在变量定义赋值之后，会作为内存变量保存并显示在Workspace Browser中。因此，凡是显示在Workspace Browser中的变量 都是“有效”的，其后可以被调用，否则不能被调用。 d、对于像 等常用的数学常量，MATLAB定义了预定义变量与其对应，在使用时需多加留意。 e、复数和复数矩阵的表示方法。 （5）其他操作的操作要旨和操作技巧的运用。 3、计算结果的图形表示 计算结果可视化是MATLAB的主要组成部分，借助图形表现数据是十分常用的“数据表达手段”，尤其当数据量相当庞大时，因为图形可以表现数据内在联系和宏观特征。关于MATLAB绘图的基本方法在后续章节中详细讲述，本实验主要通过示例了解MATLAB绘图的基本功能。 4、Current Directory、路径设置器和文件管理 理解当前目录Current Directory和搜索路径的作用是正确使用MATLAB的关键环节。当前目录指的是当前MA TLAB工作的目录，MATLAB运行指令需要打开或者保存的文件，都首先在目录中查找或保存。搜索路径则是MATLAB工作时，需查找相应的文件、函数或变量所在的相关文件夹所在的路径。 在理解当前目录Current Directory和搜索路径的作用的基础上，也要掌握当前目录Current Directory和搜索路径的设置方法，这是正确使用MA TLAB 的必要步骤。 为了理解MATLAB当前目录Current Directory和搜索路径的作用，可以大致了解一下当用户从指令窗送入一个名为cow的指令后，MATLAB的“运作次序”： （1）MATLAB在内存中检查，看cow是不是变量；如果不是，进行下一步； （2）检查cow是不是内建函数；如果不是进行下一步； （3）在当前目录下，检查是否有名为cow的M文件存在；如果不是，进行下一步； （4）在MA TLAB搜索路径的其他目录下，检查是否有名为cow的M文件存在。

南邮matlab软件设计实验报告

通信与信息工程学院 2014 /2015 学年第2学期 软件设计实验报告 模块名称MATLAB软件设计专业通信工程 学生班级通工2班 学生学号 学生姓名 指导教师孟田甜

报告内容 一、要求练习的实验部分 1． 题目：在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。 实现方法：plot 画图 程序：clc clear t=0:0.01:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); xlabel('t'); ylabel('y'); title('绘制曲线'); 结果及结 果分析： 1234 5 678910 00.51 1.5 t y 绘制曲线 对源程序的改进：1.加了坐标轴标示，标题 2.将步长细化 2． 题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦

波，要求它的负半波被置零，且在 2 3 处被削顶。 图 E2-1 题15所要表示的波形（逐段解析函数） 实现方法：plot 画图 程序 t=0:0.01*pi:10; y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(y

金融MATLAB实验报告三答案详解

安徽财经大学金融证券实验室实验报告 实验课程名称《金融》TLABMA 金融学院部课系开

级班 学号 姓名 师导指教日年月 1．

2 一、期权定价分析 1.black-scholes方程求解

例1：假设欧式股票期权，六个月后到期，执行价格90元，现价为102元，无股利支付， 股价年化波动率为55%，无风险利率为8%，计算期权价格。 解：clear Price=102; >>Strike=90; >>Rate=0.08; >>Time=6/12; >>V olatility=0.55; [CallDelta,PutDelta]=blsprice(Price,Strike,Rate,Time,V olatility) 计算结果： CallDelta= 23.5648 PutDelta= 8.0358 2.期权价格与波动率关系分析 Price=102; >>Strike=90; >>Rate=0.08; >>Time=6/12; V olatility=0.08:0.01:0.5; >>N=length(V olatility) Call=zeros(1,N); Put=zeros(1,N); for i=1:N [Call(i),Put(i)]=blsprice(Price,Strike,Rate,Time,V olatility(i)); N= 43 end plot(Call,'b--'); hold on plot(Put,'b'); xlabel('V olatility') ylabel('price') legend('Call','Put')

南邮MATLAB数学实验答案(全)

第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

matlab实验报告

MATLAB 数学实验报告 指导老师： 班级： 小组成员： 时间：201_/_/_

Matlab第二次实验报告 小组成员： 1题目：实验四，MATLAB选择结构与应用实验 目的：掌握if选择结构与程序流程控制，重点掌握break，return，pause语句的应用。 问题：问题1：验证“哥德巴赫猜想”，即：任何一个正偶数（n>=6）均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。 问题分析：由用户输入一个大于6的偶数，由input语句实现。由if判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数，再向屏幕输出两个质数即可。 编程：function [z1,z2]=gede(n); n=input('please input n')

if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n-i) break end

end end 结果分析： 如上图，用户输入了大于6的偶数返回两个质数5和31，通过不断试验，即可验证哥德巴赫猜想。 纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰，更快的解决问题。 2题目：实验四，MATLAB选择结构与应用实验 目的：用matlab联系生活实际，解决一些生活中常见的实际问题。