实验四matlab数值计算。
一、实验目的。
1.掌握线性方程组的求解方法。
2.掌握数值插值与曲线拟合的方法及应用3.掌握求数值导数和数值积分的方法4.掌握非线性方程组的求解方法。
二、实验内容。
1.求解线性方程组2.线性插值与曲线拟合3.求数值导数和数值积分4.非线性方程组的数值解。
三、实验过程。
1.求解线性方程组。
6x5y2z5u49xy4zu13
3x4y2z2u13x9y2u11
1）利用矩阵除法、逆矩阵、克莱姆法则。
a=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];b=[-4,13,1,11]';x1=a\bx2=inv(a)*bd1=[b,a(:,2:4)];d2=[a(:
,1:1),b,a(:,3:
4)];d3=[a(:,1:2),b,a(:
,4:4)];d4=[a(:,1:
3),b];dd=det(a);y1=det(d1)/dd;y2=det(d2)/dd;y3=det(d3)/dd;y4=det(d4)/dd;x3=[y1,y2,y3,y4]
2）使用命令linsolve(a,b)（3）lu分解[l，u]=lu(a)y=inv(u)*inv(l)*b
2.某气象观测站测得某日6：00~18：00之间每隔断小时的室内外温度如表所示：
时间室内温度室外温度。
试用三次样条插值分别求出该室内外6：30~17：30之间每隔断小时各点的近似值。
x=6:2:18;y=[18 20 22 25 30 28 24];z=[15 19 24 28 34 32 30];xi=6.
30:2:17.
30;yi=interp1(x,y,xi,’spline’);zi=interp1(x,z,xi,’spline’);plot(x,y,xi,yi)plot(x,z,xi,zi)
3.已知lg(x)在[1，101]区间11个整数采样点的函数值，试求lg(x)的次拟合多项式。x=1:10:101;
y=[0,1.0414,1.3222,1.
4914,1.6128,1.7076,1.
7853,1.8513,1.9085,1.
9590,2.0043];[p,s]=polyfit(x,y,5)4.求微分方程的数值解。
xd2ydy
5y02dxdx
y(0)0y(0)0
先将问题转化为一阶微分方程，然后编写函数文件，再调用ode45函数5.求i
ln(1x)
dx01x2
1）建立函数文件，然后再调用quad或者quad8命令求解2）或者使用内联函数。g=inline(‘log(1+x)/(1+x.^2)’)quad8(g,0,1)
6、求方程3xsinxex0在1.5附近的根。
编制函数文件fz:function f=f(x)f=3*x+sinx-exp(x);fxero(‘fz’,x0)
sinxy2lnz70
7、求3x2yz310在（1，1，1）附近的数值解。
xyz50编制方程组的函数文件。
在命令窗口中输入如下命令：x=fsolve('vvv','1,1,1]')
第4讲实验一
实验一项目开发的准备工作。实验学时 2实验类型 验证性。一 目的与任务。目的 确定课题，组织组员，合理分工，熟悉软件开发环境，培养团队精神。任务 学习软件开发小组的组织和管理，合理分工，将项目开发各阶段的任务明确，并熟悉相应的软件开发环境。二 参考题目及要求。1 学生信息管理系统。2 图书信息管理系...
第4讲概述
关于两个管理层次在工程建设各个阶段的咨询业务，表1 2进行了对比 项目管理层次阶段性工作和工程咨询业务关系表。2007年真题 项目执行管理层次委托的融资咨询是项目周期中 的工作。实施阶段 b 完工阶段 c 准备阶段 d 策划阶段。答案 c2007年真题 在项目实施阶段为项目决策管理层次提供的工程咨询...
第4讲概述
关于两个管理层次在工程建设各个阶段的咨询业务，表1 2进行了对比 项目管理层次阶段性工作和工程咨询业务关系表 重点每年都考 2007年真题 项目执行管理层次委托的融资咨询是项目周期中 的工作。实施阶段 b 完工阶段 c 准备阶段 d 策划阶段。答案 c2007年真题 在项目实施阶段为项目决策管理层次...