计算应用软件

MATLAB程序设计及绘图

班级： 学号： 姓名：

2016.12

计算应用软件. 工程力学专业

题目1：分别利用For循环和While循环求下值：K??2i?1?2?22???263。

063For循环代码： k=0; for i=0:1:63 k=k+2^i; end k

While循环代码: k=0; i=0; while (i<64) k=k+2^i; i=i+1; end k

结果：k = 1.8447e+019

?3x?4y?7z?12w?4?5x?7y?4z?2w??3?题目2：已知方程组：?，求方程组的解。

?8z?5w?9?x???6x?5y?2z?10w??8代码：

A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; B=[4;-3;9;-8]; X=A\\B

2 / 11

计算应用软件. 工程力学专业

结果：X = -1.4841 -0.6816 0.5337 -1.2429

题目3：某试验试件在受火后，其弹性模量发生变化，数据如下： 受火温度T=[20, 100, 200, 300, 400, 500, 600, 700, 800, 900];

受火后弹性模E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.32]*1011; 用合适的多项式拟合弹性模量随受火温度变化的规律。 代码：

T=[20,100,200,300,400,500,600,700,800,900];

E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.28]*10^11; plot(T,E,'LineWidth',4) 结果：

3 / 11