数学软件实训2-MATLAB程序综合设计及应用

数学软件实训任务二

一 题目:MATLAB程序综合设计及应用

二 目的:熟练掌握MATLAB程序设计的基本方法,会根据MATLAB程序设计的

三 要求:

        1 熟练掌握控制流的基本语法结构。

   2 会熟练应用for循环,while循环及分支语句if …else和switch…case等进行混合编程,处理实际问题。

四 实训内容:

1 方程求根的二分法及程序设计

  (1)方程求根的二分法的基本原理

取[a,b]的中点$\xi \displaystyle _{1}=\frac{a+b}{2}$,计算 $f(\xi _{1})$.

如果$f(\xi _{1})=0$;那么$\xi=\xi_{1}$;

如果 $f(\xi _{1})$与$f(a)$同号,那么取$a_{1}=\xi_{1}$,$b_{1}=b$,由$f(a_{1})*f(b_{1})<0$,即知$a_{1}<\xi <b_{1}$,且$b_{1}-a_{1}=\frac{b-a}{2}$;

如果$f(\xi _{1})$与$f(b)$同号,那么取$a_{1}=a,b_{1}=\xi _{1}$,由$f(a_{1})*f(b_{1})<0$,即知$a_{1}<\xi <b_{1}$,且$b_{1}-a_{1}=\frac{b-a}{2}$;

总之,当$\xi \neq \xi_{1}$时,可求得$a_{1}<\xi<b_{1}$,且$b_{1}-a_{1}=\frac{b-a}{2}$.

以[a1,b1]作为新的隔离区间,重复上述做法,当$\xi\neq\xi_{2}=\frac{a_{1}+b_{1}}{2}$时,可求得

$a_{2}<\xi <b_{2}$,且$b_{2}-a_{2}=\frac{b-a}{2^{2}}$.

如此重复n次,可求得$a_{n}<\xi<b_{n} $,$b_{n}-a_{n}=\frac{b-a}{2^{n}}$.由此可知,如果以$a_{n}$或$b_{n}$作为$\xi$的近似值,那么其误差小于$\frac {b-a} {2^{n}}$.如图14-23是二分法的原理示意图.

 

图 14-23  二分法原理示意图

二分法的优点是对函数的要求低(只要求满足零点定理的条件),方法简便、可靠,程序设计容易,事先估计次数容易,收敛速度恒定。

   (2)二分法程序设计

   (3)绘出函数$f(x)=x^{3}+1.1x^{2}+0.9x-1.4$的图像,找出根的隔离区间,然后用所编写的二分法程序求方程$f(x)=x^{3}+1.1x^{2}+0.9x-1.4=0$实根的近似值,使误差不超过$10^{-3}$.

 

转载于:https://www.cnblogs.com/jiaxinwei/p/8151471.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
仿真软件应用实践1 MATLAB基础实验报告(代码总结心得)一、实验目的和要求: 1、掌握数组与矩阵的创建、矩阵运算与数组运算、基本元素群运算和向量与矩阵的特殊处理; 2、掌握二维图形、三维曲线和曲面图形的绘制以及标注; 二、实验内容: 1、利用“:”操作符和 linspace 函数生成数组 a=[0,6,12,18,24,30,36,42],并将数组 a 分别转化为 4*2 和 2*4 的矩阵。 2、在MATLAB中输入矩阵A=[■(3&6&8&1@-5&7&22&17@6&9&16&-12@15&13&-21&0)],并使用MATLAB回答以下问题: (1)创建一个有A中第2列到第4列所有元素组成的4*3数组B; (2)创建一个有A中第3行到第4行所有元素组成的2*4数组C; (3)创建一个有A中前两行和后三列所有元素组成的2*3数组D; (4)根据A利用单下标方法和双下标方法分别创建向量a=[-5 6 15],b=[6 8 1],并利用向量a和b生成矩阵E=[■(-5&6@6&8@15&1)];并利用“[]”删除矩阵A的第二行和第三列。 3、已知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]; (1)求a+b,a*b,a.*b,a/b,a./b,a^2,a.^2的结果,并观察分析运算结果; (2)求c中所有元素的平均值、最大值、中值,并将c中元素按从小到大顺序排列; (3)解释b(2:3,[1,3])的含义,并将生成的数组赋值给d。 4、已知A=[■(7&2&1&-2@9&15&3&-2@-2&-2&11&5@1&3&2&13)],用MATLAB求矩阵A的秩(rank)、行列式(determinant)、逆(inverse)、特征值及特征向量(eigenvalue and eigenvector)。 5、用MATLAB求方程组[■(7&2&1&-2@9&15&3&-2@-2&-2&11&3@1&3&2&13)]*[■(x_1@x_2@x_3@x_4 )]=[■(4@7@-1@0)]的解。 6、用不同线型和颜色在同一坐标内绘制曲线y=2e-0.5xsin(2πx)及其包络线,并为图形添加必要的图形标注。 7、在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线,并为图形添加必要的图形标注。 8、在xy平面内选择区域[-8,8]*[-8,8],绘制函数z=(sin√(x^2+y^2 ))/√(x^2+y^2 )的三种三维曲面图。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值