发布时间 : 星期六 文章MATLAB结课作业(2016.12)更新完毕开始阅读6eb96eeebdeb19e8b8f67c1cfad6195f312be889
计算应用软件
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