当前位置：文档之家 › MATLAB实习题

MATLAB实习题

  • 格式：doc
  • 大小：363.00 KB
  • 文档页数：13

MATLAB实习题

1 矩阵运算

矩阵建立，转置、加、减、乘、求逆等运算。

(1) 矩阵的建立：

>> a=eye(5);

b=1*a;

b(1,2)=2;b(2,3)=3;b(3,4)=4;b(4,5)=5;

b

运行结果：

b =

1 2 0 0 0

0 1 3 0 0

0 0 1 4 0

0 0 0 1 5

0 0 0 0 1

（2） 矩阵的转置：

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

b=a'

运行结果：

a =

1 2 3

4 5 6

7 8 9

b =

1 4 7

2 5 8

3 6 9

（3）矩阵的逆：

>> a=[1 2 0;0 3 4;0 0 5];

inv(a)

运行结果：

ans =

1.0000 -0.6667 0.5333

0 0.3333 -0.2667

0 0 0.2000

（4）矩阵的加：

>> A=[1 1;1 1]

A =

1 1

1 1

>> B=[2 2;2 2]

B =

2 2

2 2

>> C=A+B

运行结果;

C =

3 3

3 3

(5)矩阵的减：

>> A=[1 1;1 1]

A =

1 1

1 1

>> B=[2 2;2 2]

B =

2 2

2 2

>> C=A-B

运行结果;

C =

-1 -1

-1 -1

(5)矩阵的乘：

>> A=[1 1;1 1]

A =

1 1

1 1

>> B=[2 2;2 2]

B =

2 2

2 2

>> C=A*B

运行结果;

C =

4 4

4 4

2 M文件的编写与运行

判断一年是否是闰年 function leapyear(year)

%使用该程序判断是否是闰年

sign=0;

if mod(year,4)~=0

sign=0;

elseif mod(year,400)==0

sign=1;

elseif mod(year,100)==0

sign=0;

else

sign=1;

end

if sign==1

fprintf('%4d year is a leap year!\n',year)

else

fprintf('%4d year is not a leap year!\n',year)

end

执行结果：

>> leapyear(1996)

1904 year is a leap year!

>> leapyear(1900)

1900 year is not a leap year!

3 数据文件的建立与读取

例：建立班级一门课成绩数据库，编程读取并排序、另存。

>> score=[80,82,81,65,77,67,55,76,90];

>> Fid=fopen('test.dat', 'w');

>> cnt=fwrite(Fid, score, 'float');

>> fclose(Fid);

>> Fid=fopen('test.dat', 'r')

>> [B,cnt]=fread(Fid, 'float');

>> B

>> sortScore=sort(B,1)%排序

结果：

score=[80,82,81,65,77,67,55,76,90];

Fid=fopen('test.dat', 'w');

cnt=fwrite(Fid, score, 'float');

fclose(Fid);

Fid=fopen('test.dat', 'r')

[B,cnt]=fread(Fid, 'float'); B

sortScore=sort(B,1)%排序

Fid =

3

B =

80

82

81

65

77

67

55

76

90

sortScore =

55

65

67

76

77

80

81

82

90

4 绘二维、三维图

(1)画双纽线

>> p=-2*pi:0.01:2*pi;

r=sqrt(cos(2*p));

x=r.*cos(p);

y=r.*sin(p);

plot(x,y)

运行结果;

(2)画螺线:

a=2;

t=0:.1:20;

r=a*t;

polar(t,r)

运行结果:

(3)画球面:

>> [x,y,z]=sphere(50);

surf(x,y,z)

运行结果：

5 多项式运算

>> p=[2,-1,0,3]

p =

2 -1 0 3

>> q=[2,1]

q =

2 1

>> k=conv(p,q)

运行结果：

k =

4 0 -1 6 3

6 曲线拟合

例：数据的多项式拟合。

多项式拟合

x=[0 30 60 90 120 150 180 210 240 270 300 330 360]

y=[-0.0167 -1.0927 -1.8725 -2.3586 -2.3061 -1.9576 -0.9574 -0.0080 0.8896

1.3877 1.1139 0.8517 -0.0167]

通过y=-0.52-1.90sin（2πt/360-0.079）进行二次多项式拟合

>> x=[0 30 60 90 120 150 180 210 240 270 300 330 360];

y=[-0.0167 -1.0927 -1.8725 -2.3586 -2.3061 -1.9576 -0.9574 -0.0080 0.8896

1.3877 1.1139 0.8517 -0.0167];

fun=@(a,t) a(1)+a(2)*sind(t+a(3))

a0=[-0.5 -1.9 -0.079];

a=nlinfit(x,y,fun,a0)

t=0:5:360;

yf=fun(a,t);

plot(x,y,'o',t,yf)

fun =

@(a,t) a(1)+a(2)*sind(t+a(3))

a =

-0.5239 -1.8995 -14.2382

运行结果:

7 线性规划

例：线性规划求解（如教材中）。

minz=32132xxx

0,,62382432121321xxxxxxxx

>> c=[2;3;1];

a=[1,4,2;3,2,0];

b=[8,6];

[x,y]=linprog(c,-a,-b,[],[],zeros(3,1))

Optimization terminated.

运行结果：

x =

0.8066

1.7900

0.0166 y =

7.0000

8 微分方程数值解

例：求微分方程yxy的通解

>> dsolve('Dy=x+y','x')

ans =

-x-1+exp(x)*C1

9 插值问题

例：样条插值。

给定函数 f(x)=1/(1+x^2)，取插值节点为 -5:1:5，利用编制的 Matlab 程序求出三次样条函数，并绘制出准确函数 f(x) 以及所求出的三次样条函数图像，进行比较。

x = -5:0.001:5;

f = 1./(1 + x.^2);

xi = -5:1:5;

fi = interp1(x, f, xi, 'spline');

figure;

plot(x, f); hold on; plot(xi, fi, 'r.-'); xlabel('x');ylabel('fx');

legend('Raw', 'Interp');

运行结果：

10综合练习：计算机模拟

例1追逐问题: 如图,正方形ABCD的四个顶点各有一人.在某一时刻,四人同时出发以匀速v=1米/秒按顺时针方向追逐下一人,如果他们始终保持对准目标,则最终按螺旋状曲线于中心点O.试求出这种情况下每个人的行进轨迹.

O BCDA

合集下载

信息与计算科学系Matlab软件实习(论文)——非线性方程求根

信息与计算科学系Matlab软件实习(论文)——非线性方程求根

Matlab软件实习论文

非线性方程求根

系 别 信息与计算科学

专 业 信息与计算科学

学 号

姓 名

指导教师

2008年 8 月 10 日

信息与计算科学系Matlab软件实习（论文） 第 I 页

非线性方程求根

摘 要

随着科学技术，生产力经济的发展，在科学与工程计算中存在着大量方程求根问题，例如贷款购房问题，工厂的最佳订货问题等都需要求解一类非线性方程的根，而本文就针对这些求根问题提出了解决方案，本文利用牛顿迭代法来结决方程的求根问题．

首先根据实际问题列出数学模型，确定变量，给出各个条件及相关函数；然后对建立的模型进行具体分析和研究，选择合适的求解方法；编写函数的程序，用计算机求出方程的解，通过所求解分析具体情况．

关键词：非线性方程，牛顿迭代法，Matlab

信息与计算科学系Matlab软件实习（论文） 第 II 页

目 录

摘 要 ........................................................................................................................................ I

1 绪 论 .................................................................................................................................... 1

1.1 非线性方程求根的背景 .................................................................................................. 1

matlab实训心得体会（通用23篇）

matlab实训心得体会（通用23篇）

matlab实训心得体会（通用23篇）

matlab实训 篇1

自己刚刚接触matlab有半个学期的时间，说实话我现在对MATLAB还是摸不着头脑，一方面是自己接触的时间太短，另一方面，就是自己在上机方面投入的时间有限，实践比较少。现在，我对MATLAB的印象仅仅在解决习题和绘制图形上，但是我很喜欢MATLAB的简单的语法,易于绘制图形,编程也非常容易, 并且具有功能强大的开放式的toolbox。 因此,尽管我一直没有这方面的应用,但是我还是对 它非常感兴趣，自己正打算暑假好好研究研究MATLAB。下面是我学习MATLAB在理论和实践方面的一点心得与体会，可能有些地方自己理解的不是很正确，但是随着学习的深入，我想我可以发现自己的错误所在。

首先我想说的是，在理论方面，在学习MATLAB过程中，我感觉到它和c语言有许多相似之处，他有c语言的特征，但是比c语言编程计算更加简单，适合于复杂的数学运算。但是MATLAB跟其他语言也有着很大的不同。现在用的比较多的编程语言，除了MATLAB就应该是c、c＋＋、VHDL，VB和Delphi也接触过，如果自己抱着“把其他语言的思想运用在MATLAB里面”的话，那么我想，即使程序运行不出错，也很难把握MATLAB的精髓，也就很难发挥MATLAB的作用了。 众所周知MATLAB是一个基于矩阵运算的软件，但是，真正在运用的时候，特别是在编程的时候，许多人往往没有注意到这个问题。在使用MATLAB时，受到了其他编程习惯的影响，特别是经常使用的C语言。因此，在MATLAB编程时，for循环（包括while循环）到处都是。.这不仅是没有发挥MATLAB所长，还浪费了宝贵的时间。我这里想说的一点是，往往在初始化矩阵的时候注意到这个问题，懂得了使用矩阵而不是循环来赋值，但是，在其他环节上，就很容易疏忽，或者说，仍然没有摆脱C＋＋、C的思想。

MATLAB博大精深，涉及的内容很多，所以，我认为不要试图掌握MATLAB的每一个功能，熟悉和你专业最相关的部分就可以了，这也是老师在课堂上经常说的。学MATLAB并不难，难的是学会怎么用，所以经常上机实践是很必要的。我自己感觉学习MATLAB和以前的编程能力没有太多的关系，所以不要担心自己编程能力差，自己一定用不好MATLAB，只要自己肯在这上面花费时间和精力，就一定能有所收获。在学习MATLAB的过程中，不要只问不学，并且学MATLAB要有耐心，要大胆的去试，哪怕只有一丁点儿可能，只有自己动手去实践了才能发现错误的所在，利用这个解决问题的方法要试着解决类似的问题，要举一反三，要学会变通。多读MATLAB高手写的程序，找到一个高手多向他请教这方面的问题，在几个大的论坛可以搜索出一大堆的帖子，然后慢慢去看吧，从中可以学到很多东西。善于总结，学习过的知识，看过好的经验介绍可以收藏起来，过段时间再复习一下，一段时间的积累，你会发现你的水平在慢慢提高（这属于实践方面的体会）还有一点比较重要多用help，see also， lookfor， get,

MATLAB)课后实验答案[1]

MATLAB)课后实验答案[1]

实验一 MATLAB运算基础

1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。

(1) 0122sin851ze

(2) 221ln(1)2zxx，其中2120.455ix

(3) 0.30.330.3sin(0.3)ln,3.0,2.9,,2.9,3.022aaeeazaa

(4) 2242011122123ttzttttt，其中t=0:0.5:2.5

解：

M文件:

z1=2*sin(85*pi/180)/(1+exp(2))

x=[2 1+2*i;-.45 5];

z2=1/2*log(x+sqrt(1+x^2))

a=-3.0:0.1:3.0;

z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)

t=0:0.5:2.5;

z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)

4. 完成下列操作：

(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量，删除其中的大写字母。

解：(1) 结果：

m=100:999;

n=find(mod(m,21)==0);

length(n)

ans =

43

(2). 建立一个字符串向量 例如：

ch='ABC123d4e56Fg9';则要求结果是：

ch='ABC123d4e56Fg9';

k=find(ch>='A'&ch<='Z');

ch(k)=[]

ch =

实验二 MATLAB矩阵分析与处理

1. 设有分块矩阵33322322ERAOS，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22ERRSAOS。

《数值计算方法实习》教学大纲

《数值计算方法实习》教学大纲

《数值计算方法实习》教学大纲

Numerical Computation Method Practice

适用本科四年制信息与计算科学专业（2周 2学分）

一、课程的目的和任务

本课程的授课对象是信息与计算科学专业本科生，属信息与计算科学专业公共基础课。

数值计算方法是一门专门研究各种数学问题近似解法的课程，它是一门与计算机应用密

切结合的实用性很强的数学课程。在数值计算方法课程中，讲授了各种数学问题的近似解法，这些近似解法的计算量很大，只有利用计算机计算，这些解法才具有实用意义。因而上机实习，掌握这些近似解法的计算机实现是数值计算方法课程学习的一个重要环节。

本课程实习的主要目的是通过科学计算语言MATLAB的学习，利用MATLAB求解各种数学问题的近似解，使学生对数值计算方法课程所学的各种近似解法能在计算机上实现，提高学生对数值计算方法课程讲授的各种数学问题近似解法的理解和掌握。

通过本实践环节，要求学生初步掌握MATLAB的使用方法，掌握利用MATLAB求解各种数学问题近似解的算法，通过上机实践，提高学生对各种数学问题近似解法的实际运用能力，并能应用所学的方法解决一些较简单的实际问题。

二、课程的基本要求和特点

本课程是一门既有系统理论又有较强实践性的技术基础课，学习本课程需坚持理论联系实际的学风，必须在学习数值计算方法课程讲授的各种数学问题近似解法的基础上，动手编写一些简单的MATLAB程序，利用MATLAB来实现求解各种数学问题的近似解；同时要注意数学软件的使用原理及使用方法。本课程是一门实用性很强的应用数学课程。

三、本课程与其它课程的联系

本课程实习是对前期《数值计算方法》课程的巩固，数值计算方法课程涉及面较宽，必须先修课程为《数学分析》、《高等代数》、《常微分方程》、《计算机应用基础》、《数值计算方法》。

四、课程的主要内容

1 数学软件MATLAB

教学要求：

了解：MATLAB的基本特点，MATLAB的启动方法和工作界面，MATLAB数值计算，