用最小二乘法计算同一多项式会得出不同的显示结果。
实验三
一、实验题目
复合梯形求积
对于函数fxsi
xx利用下表计算此函数在01区间的积分Xfx011809973978140989615838097672671209588510580936155634090885167808771925108414709
5
f二、程序功能
根据复合梯形求积公式,计算函数的积分
三、算法
步1:计算步长h。即hmaxxmi
x
1;步2:计算各个结点对应的函数值。即f1
y1
;步3:将第二步所得结果代入复合梯形公式,即If12sumf2
1f
h2;步4:计算结果。
四、重要标识符说明x函数中的自变量;xi每个结点处的自变量值;y每个自变量对应的函数值;h步长;
max是取最大值;mi
是取最小值;sum是求和函数;
五、程序运行实例
在matlab65环境中,运行程序,结果如图:
六、源程序
ComTrapzmclear
9x0181438125834781y1099739780989615809767267095885100936155609088516
6
f0877192508414709hmaxxmi
x
1f1
y1
If12sumf2
1f
h2dispI
七、个人实验总结
在matlab中用了很多矩阵和向量来解决问题,更加方便,不用再定义这些矩阵,直接输入就可以。掌握节点的个数和公式,了解一些内函数的相关知识,即可理解整个程序。
实验四
一、实验题目
二分法
求方程fxx3x10在区间10,15内的一个实根。
二、程序功能
利用二分法,求方程的根
三、算法
步1设区间为ab先把ab2代入fx,若为0则根xab2,若不为0看符号。步2:若fafab20则根在区间(aab2)内,否则根在区间(ab2b)内,继续在有根区间内二分,一直到达到精度。步3:循环步骤1和3,直到步骤2成立
7
f步4:取有根区间的中点为方程的根。
四、重要标识符说明ab:是区间的端点;emg:达到的精度;k:循环次数;x:所求近似值;Fa:将a代入方程的值;fab:将ab2代入方程的值;
feval:是求值函数;abs():是求绝对值函数;:函数调用;
五、程序运行实例
在matlab65环境中,运行程序,结果如图:
8
f六、源程序头文件fu
ctio
ffu
xfx3x1源程序二分法fu
ctio
xkdemimethodabfemgfafevalfafabfevalfab2k0whileabsbaemgiffab0xab2retur
elseiffafab0bab2elseaab2e
dfafevalfafabfevalfab2kk1e
dxab2七、个人实验总结
第一次用到了函数调用的方法,起初以为把它们写到一个文件里即可,后来才知道在程序中若要用到函r