数值计算实验题目.doc

《数值计算实验题目.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验1Matlab环境一、实验目的熟悉Matlab操作环境和编程环境，Matlab数组变量的定义与使用二、实验内容与步骤<1）启动MATLAB软件，熟悉软件的界面和菜单<2）矩阵的运算1）矩阵的乘法已知A=[12。34]。B=[54。32]。观察A.*BA*BA.^BA’的计算结果2）矩阵除法已知A=[123。456。789]。B=[100。020。003]。观察A＼B,A/B的计算结果3）使用冒号选出指定元素已知：A=[123。456。789]。求A中第3列前2个元素；A中所有列第2，3行的元素；4）方括

2、号[]已知A=[123。456。789]。，删除该矩阵的第三列5>.输入下列特殊矩阵,并观察结果〉〉A=[]〉〉A=ones(5,10>>>A=zeros(5,10>6）已知X=0:pi/5:2*pi,观察sin(X>的执行结果。7>编写一个Matlab命令文件test.m,并将其保存在C:目录下。该文件内容如下：U1=10:20。U2=10:2:20。U3=linspace(10,20,6>。U4=[U2U3]。请在Matlab命令窗口执行该文件，并观察U1，U2，U3，U4的内容。注意文件的路径问题。

3、---------------------------------------------------------------------b5E2RGbCAP实验2Matlab编程与调试1、编写命令文件：计算1+2+…+n<2000时的最大n值；2、编写函数文件：分别用for和while循环结构编写程序，求---------------------------------------------------------------------p1EanqFDPw实验3Matlab函数图形绘制一实验目的<1）

4、学会用Matlab软件做常用函数的图形4/4<2）通过作图，进一步加深对函数的理解，观察函数的性质二预备知识Matlab拥有大量简单、灵活、易用的二维、三维图形命令，许多图形命令都在Matlab所带的演示程序中给出。DXDiTa9E3d（一）一元函数作图1.基本一元函数作图2.参数方程作图二、实验内容与步骤（一）一元函数作图1）已知X=[10962]。Y=[10081364]。请作图形表示出X与Y的关系。RTCrpUDGiT2）在同一图形窗口用红线画出sin(x>,绿圈画出cos(x>的曲线图，其中.在图形

5、上加格栅、图例和标注。5PCzVD7HxA3）定义分段函数，并绘出f(x>的图形。<二）参数方程作图画出星形线：的图形---------------------------------------------------------------------jLBHrnAILg实验4方程求根下图所示电路中，已知L=5H，C=F，开关合上后，电容器上t时刻的电荷，其中为电容器上的初始电荷，试用二分法或newton迭代法求得电阻R，使电容器上的电荷在开关合上后0.05s内衰减为原来值的1%,即，要求误差小于0.0

6、001%，并比较两种方法的收敛速度。<电阻初值取200）xHAQX74J0X---------------------------------------------------------------------LDAYtRyKfE实验题5插值方法4/4目的：对比各种插值方法题目：求出y=1./(1+x.^2>的lagrange插值、分段线性插值、曲线拟合结果，并把结果用图形显示在同一个图形窗口，并进行分析。Zzz6ZB2Ltk补充知识：实现分段线性插值不需另外编制函数程序，MATLAB自身提供了内部函数

7、interpl。调用格式如下yi=interpl(x，y，xi>。对一组节点(x，y>进行插值，计算插值点xi的函数值。x为节点向量值，y为对应的节点函数值。如果y为矩阵，则插值对y的每一列进行，若y的维数超出x或Xi的维数，则返回NaN。dvzfvkwMI1实现曲线拟合亦不需另外编制函数程序，MATLAB自身提供了内部函数polyfit。调用格式如下c=polyfit(x，y，n>。对一组节点(x，y>进行拟合，计算拟合系数c。x为节点向量值，y为对应的节点函数值。rqyn14ZNXI----------

8、-----------------------------------------------------------EmxvxOtOco实验题6曲线拟合的实际应用题目：用电压V=10伏的电池给电容器充电，电容器上t时刻的电压，其中是电容器的初始电压，是充电常数。实验测量了一组了t和n(t)的数据。请根据该数据表确定和的大小。SixE2yXPq5t(秒>0.51234579n(t)(伏>6.366.487.268