MATLAB 实验7

MATLAB 实验报告

一、实验目的

熟悉和掌握MA TLAB 的m 文件程序和程序设计。

二、实验内容

1．for 循环结构

for i=1:10;

x(i)=i;

end;

x

这是通过“for ”循环产生从1到10的是个数字。执行结果如下图：

2．while 循环结构

Fibonaci 数组的元素满足Fibonacci 规则：);,2,1k (,a a a 1k k 2k =+=++且1a a 21==。现要求该数组中第一个大于10000的元素。输入命令后执行结果如下图：

a(1)=1;a(2)=1;i=2;。

while a(i)<=10000

a(i+1)=a(i-1)+a(i);

i=i+1;

end;

i ,a(i)

3.if-else-end 结构。这种结构就是对于多种可能的情况进行不同情况的划分，然后对

不同情况进行编辑执行不同的命令。输入命令执行结构如下：

cost=10;number=12;

if number>8

sums=number*0.95*cost;

end,sums

4.M 文件的一般结构

(A)编写一个画任意半径彩色线型的圆.(B)完整函数文件的基本结构.(C)函数文件各基本组成部分的作用.

如下是M 函数，只需要调用即可,执行结果如下：

function sa=functiondemo(r,s)

% CIRCLE

if nargin>2

error('输入变量太多。');

end

if nargin==1

s='b';

end;

clf;

t=0:pi/100:2*pi;

x=r*exp(i*t);

if nargin==2

plot(x,s);

else

sa=pi*r*r;

fill(real(x),imag(x),s)

end

axis('square')

比如输入“ functiondemo(4,'red')”即可得到右图。

三、思考题

1．编写函数文件，分别用for 和while 循环语句编写程序，求出 63654363

0i 2i 22222212K +++++++==∑= 之和，并比较算法的运算时间。

答：（输入一个n 值，就可以得到2的0次方一直加的2的n-1次方的累加值，只需要调用函数输入n 的值为63即可完成上题）M 函数如下：

for 循环：

function s=input(n)

tic

a=0;

for i=1:n+1

x(i)=i-1;

a=a+2^x(i);

end

t=toc

a

while 循环：

function inputs(n)

tic

b=0;

i=0;

while i<=n

b=b+2^i;

i=i+1;

end

t=toc

b

得出结果：

2.用MATLAB 语言实现下面分段函数：??

???-<-≤>==1,11||,1,1)(x x x x x f y

x 由键盘输入。

答：同上面做法一样，编辑M 函数实现，然后输入调用函数

即可）

M 函数如下：

function fenduan(n)

if n>1

y=1;

elseif n>=-1&n<=1

y=n;

elseif n<-1

y=-1;

end ;

y

3.用for-end 循环求出从100-200之间的所有素数。

答：编辑M 函数实现这一命令，如下：

function sushu

m=1;

for n=100:200

if rem(n,2)~=0&rem(n,3)~=0&rem(n,5)~=0&rem(n,7)~=0&rem(n,11)~=0 &rem(n,13)~=0

x(m)=n;

m=m+1;

end

end

x

输出结果如下：

4.用while-end循环求1-100的所有数值和。

答：while-end循环M函数如下：：（输入不同的n值就可以得到1~n的累加和。）function sb(n)

p=0;

i=0;

while i<=n

p=p+i;

i=i+1;

end,p

四、遇到的问题

1. 第一个问题是出现在运算时间的计算上，对于etime函数的使用，计算出来的时间总是0，经过反复的尝试和与同学的讨论，终于找到了好的途径，将“t1=clock”、“t1=clock”命令放在正确的地方才可以执行出正确的结果。

2. 在第二个思考题中的命令“-1<=n<=1”无法正常执行，这是一个小失误，只需要用到“&”罗技操作符，改变命令为“n>=-1&n<=1”就可以正常执行了。

五、体会

这次实习的目的就是熟悉和掌握MA TLAB的m文件程序和程序设计。在这里我们要先打开MATLAB，然后建立一个空的M文件程序，在空的M文件中编辑程序就可以了。

在这次实习中我们运用编辑了许多函数，比如“if-elseif-end”、“while-end”和“switch-case”等，用以实现不同条件下的问题解决。这些函数我们在学习C语言时也曾经接触过，用过，大概知道一些编辑注意，但是MATLAB中对于这些函数的使用条件、要求明显降低了很多，使得这些函数在使用时简便了很多。不仅如此，在练习这些函数解决实际问题的同时，我们并不是直接编辑输入命令，而是把这些函数编辑到M文件中，在需要的

时候直接调用就可以了，免去了每次重复输入命令，同时也练习了M函数的编辑与调用。

这次实习出现的问题就是一些符号的使用，比如“-1<=n<=1”就无法正常执行，需要用到的是逻辑操作符“&”——“n>=-1&n<=1”就可以正常执行了。

总之，在这次实习中，就是学会了比如“if-elseif-end”函数的使用方法以及M文件的编辑与调用，让我受益匪浅。

MATLAB数学实验第二版答案(胡良剑)

数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点 >> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500

实验三 MATLAB绘图(含实验报告)

实验三 MATLAB 绘图 一、实验目的 1．掌握二维图形的绘制。 2．掌握图形的标注 3．了解三维曲线和曲面图形的绘制。 二、实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 设计提示 1．Matlab 允许在一个图形中画多条曲线：plot(x1，y1，x2，y2，……) 指令绘制y 1 = f 1(x 1), y 2 = f 2 (x 2 )等多条曲线。Matlab 自动给这些曲线以不同颜色。标注可用text 函数。 2．绘图时可以考虑极坐标和直角坐标的转换。 3．三维曲线绘图函数为plot3，注意参考帮助中的示例。 三、实验内容 1．生成1×10 维的随机数向量a ，分别用红、黄、蓝、绿色绘出其连线图、 杆图、阶梯图和条形图，并分别标出标题“连线图”、“杆图”、“阶梯图”、“条形图”。 2、绘制函数曲线，要求写出程序代码。 (1) 在区间[0:2π]均匀的取50个点，构成向量t (2) 在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求y1曲 线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号。 (3) 分别在靠近相应的曲线处标注其函数表达式。 3.将图形窗口分成两个绘图区域，分别绘制出函数： ???+-=+=1 352221x x y x y 在[0,3]区间上的曲线，并利用axis 调整轴刻度纵坐标刻度，使1y 在[0,12] 区间上，2y 在[-2,1.5]区间上。 4．用mesh 或surf 函数，绘制下面方程所表示的三维空间曲面，x 和y 的

取值范围设为[-3，3]。 10102 2y x z +-= 思考题： 1. 编写一个mcircle(r)函数，调用该函数时，根据给定的半径r ，以原点 为圆心,画一个如图所示的红色空心圆。（图例半径r=5）；左图参考 polar 函数的用法，右图绘制圆形的参数方程为x=sin （t ），y=cos （t ）。其中，t 的区间为0~2*pi ，步长为0.1。 2．（1）绘一个圆柱螺旋线（形似弹簧）图。圆柱截面直径为10，高度为5， 每圈上升高度为1。如左图所示。 （2）利用（1）的结果，对程序做少许修改，得到如右图所示图形。

数学实验练习题(MATLAB)

注意:在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上）. 第一次练习题 1.求解下列各题： 1)30sin lim x mx mx x ->- 2)(4)cos ,1000.0=x mx y e y 求 3）21/2 0mx e dx ?(求近似值，可以先用inline 定义被积函数，然后用quad 命令) 4）4 224x dx m x +? 5 0x =展开(最高次幂为8). 2.对矩阵21102041A m -?? ?= ? ?-?? ，分别求逆矩阵，特征值，特征向量，行列式，并求矩阵,P D （D 是对角矩阵），使得1A PDP -=。 3. 已知2 1(),()2f x e x μσ=--分别在下列条件下画出)(x f 的图形： (1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). 4.画 （1）sin 020cos 02100x u t t y u t u t z m ??=≤≤?=?≤≤??=?

(2) sin()03,03z mxy x y =≤≤≤≤ （3）sin()(/100cos )02cos()(/100cos )02sin x t m u t y t m u u z u π π=+?≤≤?=+?≤≤?=? 的图（第4题只要写出程序）. 5．对于方程50.10200 m x x --=，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会. 第二次练习题 判断迭代收敛速度的程序 x0=1;stopc=1;eps=10^(-8);a=1;c=1;b=2*c;d=a;k=0; f=inline('(a*x+b)/(c*x+d)'); kmax=100; while stopc>eps&k

MATLAB实验报告

MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日

实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans =

30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans =

matlab实验(西北农林科技大学)

课程实验报告 学年学期2015-2016学年第2学期 课程名称MATLAB与科学计算 实验名称实验一 MATLAB环境与数值计算实验室测量测绘实验中心计算机室专业年级土木13级 学生姓名 学生学号 20130 提交时间 2016-06-02 成绩 任课教师王增红 水利与建筑工程学院

实验一 MATLAB环境与数值计算 1、目的和要求 （1）熟练掌握MATLAB的启动和退出、MATLAB的命令窗口；常用选单和工具栏，通过实例初步认识对MATLAB的命令和程序运行等。 （2）熟练掌握MATLAB变量的使用； （3）熟练掌握矩阵的创建、掌握MATLAB的矩阵和数组的运算； （4）熟悉MATLAB多项式的运用，掌握MATLAB的拟合和插值方法、积分运算和统计运算。2、内容和步骤 参见教材实验一、二。 3、实验报告提交要求 （1）用MALTAB命令和M程序两种方式求解三元一次方程组。 答：以求解三元一次方程组{2a?b+3c=5 3a+b?5c=5 4a?b+c=9 为例： （2）用MALTAB命令和M程序两种方式求解3×3复数矩阵的模和相角。答：

（3）创建一个4╳4矩阵，求逆矩阵、产生对角矩阵、求矩阵的秩等。 （5）创建一个一元四次多项式，进行多项式求值，求根计算。 432

（6）创建一个一元三次多项式，取十个不同自变量值，计算多项式的值，并分别进行一阶、二阶拟合，图示原始数据和拟合结果 （7）自建一个一元三次多项式，取十个不同自变量值，计算多项式的值，并用线性、相邻点和三次插值，图示其插值结果。 解： （8）给出多组温度的测定值（内含平均温度、最低、最高温度），并求其平均温度、最低、最高温度的平均值、标准差等。 解： >> a=fix(40*rand(6,3)+1)

浅析Matlab数学实验报告

数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)

cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);

matlab实验十--周期函数

实验十周期函数 【实验目的】 1.了解周期函数的基本概念。 2.了解周期函数经过四则运算、复合运算、求导运算、积分运算后的周期性。 3.学习、掌握MATLAB软件有关命令。 【实验内容】 从图形上观察六个三角函数的周期性。 【实验准备】 1.周期函数的基本概念 函数() f x是以T为周期的周期函数是指对任何x，有 += ()() f x T f x 使得上式成立的最小正数T称为函数的最小正周期。 2.周期函数的四则运算 3.周期函数的最小正周期 【实验重点】 1、周期函数的四则运算与复合 2、周期函数的求导与积分运算 【实验难点】 1、最小正周期的确定 【实验方法与步骤】 练习 1 图形上观察六个三角函数

sinx,cosx,tanx,cotx,secx,cscx的周期性。绘制正弦函数y=sinx 在区间[-6π，6π]的图形，相应的MATLAB代码为 >>x=-6*pi:2*pi/30:6*pi; >>y=sin(x); >>plot(x,y); >>xlabel('x');ylabel('y'); 运行结果见图10.1。 从图形中可以看出y=sinx为周期函数，最小正周期T≈6。实际上，最小正周期T=2π=6.28…。同样，可以画出余弦函数y=cosx的图形，见图10.2，其最小正周期也为T=2π。

画正切函数y=tanx 的图形时，要注意函数在,0,1,2,...2x k k π π=+=±±是不连续，所以我们只能分别绘出函数在区间(,),0,1,2,...22k k k ππ ππ-+=±±的图形。相应的MATLAB 代码为 >>x=-1.5:0.01:1.5; >>x1=x-pi;x2=x+pi; >>y=tan(x);y1=tan(x1);y2=tan(x2); >>plot(x,y,x1,y1,x2,y2); >>xlabel('x');ylabel('y'); 运行结果见图10.3。

(完整版)Matlab实验7图形绘制

实验7：图形绘制 一、实验目的 1、 掌握绘制二维图形的常用函数。 2、 掌握绘制三维图形的常用函数。 3、 掌握绘制图形的辅助操作。 二、实验内容 1、 已知2*13),2cos(2,12 y y y x y x y ===，完成下列操作： （1） 在同一坐标系下用不同的颜色和线型绘制三条曲线。 %homework_7_1_1.m x=0:pi/100:2*pi; y1=x.*x; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'r--',x,y2,'k:',x,y3,'b-.'); （2） 以子图形式绘制三条曲线。 %homework_7_1_2.m x=0:pi/100:2*pi; y1=x.*x; y2=cos(2*x); y3=y1.*y2; subplot(2, 2, 1); plot(x,y1,'r--'); subplot(2, 2, 2); plot(x,y2,'k:'); subplot(2, 2, 3); plot(x,y3,'b-.'); （3） 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。 %homework_7_1_3.m %ì?D?í??￠?×ìYí??￠??í?oíì?3?í? x=0:pi/100:2*pi; y1=x.*x; y2=cos(2*x); y3=y1.*y2; %μúò?DD subplot(4, 3, 1); bar(x,y1,'r'); subplot(4, 3, 2); bar(x,y2,'k'); subplot(4, 3, 3); bar(x,y3,'b'); %μú?tDD

stairs(x,y1,'r'); subplot(4, 3, 5); stairs(x,y2,'k'); subplot(4, 3, 6); stairs(x,y3,'b'); %μúèyDD subplot(4, 3, 7); stem(x,y1,'r'); subplot(4, 3, 8); stem(x,y2,'k'); subplot(4, 3, 9); stem(x,y3,'b'); %μú??DD subplot(4, 3, 10); fill(x,y1,'r'); subplot(4, 3, 11); fill(x,y2,'k'); subplot(4, 3, 12); fill(x,y3,'b'); 2、 绘制极坐标曲线)sin(θρn b a +=，并分析参数a ，b ，n 对曲线形状的影响。 %homework_7_2.m function homework_7_2(a,b,n) theta=0:0.01:2*pi; rho=a.*sin(b+n.*theta); polar(theta,rho,'k'); %homework_7_2_tiao.m % a μ?ó°?ì subplot(3, 4, 1); homework_7_2(1,1,1) subplot(3, 4, 2); homework_7_2(2,1,1) subplot(3, 4, 3); homework_7_2(3,1,1) subplot(3, 4, 4); homework_7_2(4,1,1) % b μ?ó°?ì subplot(3, 4, 1); homework_7_2(1,1,1) subplot(3, 4, 2); homework_7_2(1,2,1) subplot(3, 4, 3); homework_7_2(1,3,1)

南邮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实验三—迭代 二、实验目的和要求 a.熟悉MATLAB软件的用户环境，掌握其一般目的命令和MATLAB数组操作与 运算函数； b.掌握MATLAB软件的绘图命令，能够熟练应用循环和选择结构实现各种循环 选择功能； c.借助MATLAB软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜 想，发现进而证实其中的规律。 三、实验内容 问题一：将方程53 x x x +-+=改写成各种等价的形式进行迭代 5210 观察迭代是否收敛，并给出解释。 问题二：迭代以下函数，分析其收敛性。 4 f(x)=x-a 使用线性连接图、蛛网图或费根鲍姆图对参数a进行讨论和观察，会得到什么结论？ 问题一： （1）画图 x1=-6:0.01:6; x2=-3:0.01:3; x3=-1:0.01:1; x4=-0.8:0.01:-0.75; y1=x1.^5 +5*x1.^3-2*x1+1; y2=x2.^5 +5*x2.^3-2*x2+1; y3=x3.^5 +5*x3.^3-2*x3+1; y4=x4.^5 +5*x4.^3-2*x4+1; subplot(2,2,1),plot(x1,y1) ,title('图(1)') ,grid on, subplot(2,2,2),plot(x2,y2) ,title('图(2)'),grid on, subplot(2,2,3),plot(x3,y3) ,title('图(3)'),grid on, subplot(2,2,4),plot(x4,y4) ,title('图(4)') ,grid on,

Matlab实验第一次实验答案

实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 解：sqrt Square root Syntax B = sqrt(X) Description B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i

1.0000 1.4142 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B 解：A=[1 2;3 4 ]; B=[5 5;7 8 ]; A^2*B （2）矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B 解：A=[1 2 3;4 5 6;7 8 9 ]; B=[1 0 0;0 2 0;0 0 3 ]; A\B,A/B （3）矩阵的转置及共轭转置

MATLAB实验七

实验七、MATLAB的绘图操作（一） 1、>> x=0:0.1:10;y1=x.^2;y2=cos(2*x);y3=y1.*y2; 在同一坐标下用不同的颜色和线型绘制三条曲线 >> plot(x,y1,':'),hold on,plot(x,y2,'b'),hold on,plot(y3,'g') 以子图形式绘制3条曲线 >> subplot(2,2,1),plot(x,y1),subplot(2,2,2),plot(x,y2),subplot(2,2,3),plot(y3) 分别以条形图阶梯图直方图和填充图绘制曲线 条形 图：>>subplot(2,2,1),bar(x,y1),subplot(2,2,2),bar(x,y2),subplot(2,2,3),bar(y3)

阶梯图： >> subplot(2,2,1),stairs(x,y1),subplot(2,2,2),stairs(x,y2),subplot(2,2,3),stairs(y3) 填充图： >> subplot(2,2,1),fill(x,y1,'r'),subplot(2,2,2),fill(x,y2,'r'),subplot(2,2,3),fill(x,y3,'r') 2、绘制极坐标曲线p=a*sin(b+nθ) >> a=2;b=4;n=3;th=0:0.01:pi;polar(th,a*sin(b+n*th)) 3、绘制函数的曲面图 （1、）>> [x,y]=meshgrid(-1:0.01:1,-1:0.01:1);z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y);plot3(x,y,z,'y')

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联系生活实际，解决一些生活中常见的实际问 题。

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150

plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1];

for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图

2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？ 2.2.2问题分析 如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。

MATLAB)课后实验答案[1]

实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221ln(1)2 z x x =++，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 2242011 122123t t z t t t t 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基础实验

实验一 MATLAB基础知识 一、实验目的 1．熟练掌握Matlab的启动与退出； 2. 熟悉Matalb的命令窗口，常用命令和帮助系统； 3. 熟悉Matalb的数据类型； 4. 熟悉Matlab的基本矩阵操作，运算符和字符串处理 二、实验设备 1．方正电脑 2．MATLAB软件 三、实验内容 1. 已知矩阵 A=[3 4 -1 1 -9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;2 -4 5 -6 12 -8;-3 6 -7 8 -1 1;8 -4 9 1 3 0] 写出完成下列操作的命令: (1) 将矩阵A的第2-5行中第1，3，5列元素赋给矩阵B； (2) 删除矩阵A的第7号元素； (3) 将矩阵A的每个元素加上30； (4) 求矩阵A的大小和维数； (5) 将矩阵A的右下角3*2矩阵构成矩阵C； (6) 输出[-5,5]范围内的所有元素； 程序： A=[3 4 -1 1 -9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;2 -4 5 -6 12 -8;-3 6 -7 8 -1 1;8 -4 9 1 3 0] %(1) B=A(2:5,[1 3 5]) %(2) A2=A; A2(7)=[] %(3) A3=A+30 %(4) length(A) size(A) %(5) C=A(end-2:end,end-1:end) %(6) A6=A.*(A>=-5 & A<=5) 结果: A =

3 4 -1 1 -9 10 6 5 0 7 4 -16 1 -4 7 -1 6 -8 2 -4 5 -6 12 -8 -3 6 -7 8 -1 1 8 -4 9 1 3 0 B = 6 0 4 1 7 6 2 5 12 -3 -7 -1 A2 = Columns 1 through 7 3 6 1 2 -3 8 5 Columns 8 through 14 -4 -4 6 -4 -1 0 7 Columns 15 through 21 5 -7 9 1 7 -1 - 6 Columns 22 through 28 8 1 -9 4 6 12 -1 Columns 29 through 35 3 10 -16 -8 -8 1 0 A3 = 33 34 29 31 21 40 36 35 30 37 34 14

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

Matlab数学实验一2015(标准答案版)

Matｌab数学实验一——matlab初体验 一、实验目的及意义 ［1] 熟悉MAＴＬＡB软件的用户环境； ［2] 了解MAＴLAB软件的一般目的命令； ［３] 掌握MATLＡB数组操作与运算函数； 通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。 二、实验内容 １.认识ｍatｌab的界面和基本操作 2.了解mａtlab的数据输出方式（ｆormat) 3. MATLAＢ软件的数组(矩阵）操作及运算练习; 三、实验任务 根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式和要求填写实验报告 1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。 iｊeps inf nａn pi realｍaｘrealｍin 2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。 3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。 4. 写出在命令窗口中的计算步骤和运行结果。 （1)计算 1.22 10 (ln log) 81 e ππ +- ； >＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns = ０．０3４8 (2) >> x=2;y＝４； >> z=x^２+exｐ(x+y)-y*lｏg(x)-3 z ＝ 40１.65６2 （3)输入变量 13 5.3, 25 a b ?? ==?? ?? ，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0 １.ｍat”文件。测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3 ａ＝

实验十一,MATLAB模拟

实验十MATLAB模拟 一．实验目的 掌握应用matlab进行离散事件的蒙特卡罗方法模拟的一般原理和过程。 二．实验原理和方法 1．模拟的概念：模拟就是利用物理的、数学的模型来类比、模仿现实系统及其演变过程，以寻求过程规律的一种方法。模拟的基本思想是建立一个试验模型，这个模型包含所研究系统的主要特点．通过对这个实验模型的运行，获得所要研究系统的必要信息。 2．数学模拟：在一定的假设条件下，运用数学运算模拟系统的运行，称为数学模拟。现代的数学模拟都是在计算机上进行的，称为计算机模拟。 在实际问题中，面对一些带随机因素的复杂系统，用分析方法建模常常需要作许多简化假设，与面临的实际问题可能相差甚远，以致解答根本无法应用。这时，计算机模拟几乎成为唯一的选择。 蒙特卡洛（Monte Carlo）方法是一种应用随机数来进行计算机模拟的方法．此方法对研究的系统进行随机观察抽样，通过对样本值的观察统计，求得所研究系统的某些参数． 例1．在我方某前沿防守地域，敌人以一个炮排（含两门火炮）为单位对我方进行干扰和破坏．为躲避我方打击，敌方对其阵地进行了伪装并经常变换射击地点． 经过长期观察发现，我方指挥所对敌方目标的指示有50％是准确的，而我方火力单位，在指示正确时，有1/3的射击效果能毁伤敌人一门火炮，有1/6的射击效果能全部消灭敌人． 现在希望能用某种方式把我方将要对敌人实施的20次打击结果显现出来，确定有效射击的比率及毁伤敌方火炮的平均值。 分析：这是一个概率问题，可以通过理论计算得到相应的概率和期望值.但这样只能给出作战行动的最终静态结果，而显示不出作战行动的动态过程. 为了能显示我方20次射击的过程，现采用模拟的方式。 1．问题分析：需要模拟出以下两件事： [1] 观察所对目标的指示正确与否 模拟试验有两种结果，每一种结果出现的概率都是1/2． [2] 当指示正确时，我方火力单位的射击结果情况 模拟试验有三种结果：毁伤一门火炮的可能性为1/3(即2/6)，毁伤两门的可能性为1/6，没能毁伤敌火炮的可能性为1/2(即3/6)． 2. 符号假设 i：要模拟的打击次数；k1：没击中敌人火炮的射击总数； k2：击中敌人一门火炮的射击总数；k3：击中敌人两门火炮的射击总数．E：有效射击比率；E1：20次射击平均每次毁伤敌人的火炮数．

MATLAB数学实验A

clear; clc; a=1;b=1; ezplot(sprintf('x^2/%f-y^2/%f',a^2,b^2)); hold on; ezplot(sprintf('x^2/%f-y^2/%f-1',a^2,b^2)); ezsurf('sin(a)*cos(b)','sin(a)*sin(b)','cos(a)',[0,pi,0,2*pi],60); hold on; ezsurf('x^2+y^2',[-1,1,-1,1],60);

clear all; x=-8:0.1:8; y=-8:0.1:8; [X,Y]=meshgrid(x,y); Z=sin(sqrt(X.^2+Y.^2))./sqrt(X.^2+Y.^2+2); [X,Y,Z]=peaks(50); surf(X,Y,Z)

syms x y; y=2*x^3-6*x^2-18*x+7; solve(diff(y,x),x) x=-1;eval(y) x=3;eval(y)

syms x y; z='x*y'; dblquad(z,1,4,-1,2) 结果 ans = 11.2500 求函数1+x -exp(2*x)+5的原函数clear all syms x C; f=int(1+x -exp(2*x)+5,'x')+C syms x y; >> x=0:0.01:1; >> y=sin(sin(x)); >> trapz(x,y)

x=0:0.05:1; y=[1.97687 2.17002 2.34158 2.46389 2.71512 3.06045 3.27829 3.51992 3.8215 4.2435 4.55188 4.88753 5.15594 5.698 6.04606 6.42701 7.00342 7.50192 7.89178 8.49315 9.0938] cftool 解常微分方程y’=-0.9y/(1+2x)的数值解y(0)=1 从0到0. 1的数值解，取步长0.02 clear all x1=0; x2=0.1; h=0.02; y(1,1)=1;