MATLAB结课作业(2016.12)

计算应用软件

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

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

题目4：（1）编写一个函数式文件，使其产生如下函数段：

?2x?f(x)??5?0.5x?x?4?x??22?x??6

x?6（2）编写命令式文件，绘制x在[0,10]范围内f(x).f(x+2)的曲线。

（1）代码： function f=Untitled(x) if x<=2 f=2*x; elseif x>6 f=x-4; else

f=5-0.5*x; end

（2）代码： i=0;

for x=0:0.1:10 i=i+1;

y1(1,i)=Untitled(x); end i=0

for x=0:0.1:10 i=i+1;

y2(1,i)=Untitled(x+2); end y3=y1.*y2; x=0:0.1:10;

plot(x,y3,'LineWidth',5)

4 / 11