第5次作业(2014.10.23)
班级： 学号：1431025008 姓名：兴丽娟 计算机编号：30 (练习后，把文件名”第5次作业.doc”，更名为你自己的名字“学号姓名5.doc”，并上传到指定服务器ftp://219.243.52.178的submitting子目录。
1. 某高分子材料在应力松弛实验中，应力与时间的关系如下：
t (h)????g/mm?? 0.2 4.094 0.6 2.744 1.0 1.839 1.4 1.233 1.8 0.826 2.2 0.554 利用线性回归方法计算出参数??和?，绘出?与t的函数图。 ???0exp???
???把程序及运行结果写在下面： t=[0.2,0.6,1.0,1.4,1.8,2.2]';
sigma=[4.094,2.744,1.839,1.233,0.826,0.554]'; Y=log(sigma);
X=ones(length(t),2); X(:,2)=t; B=X\\Y;
sigma0=exp(B(1)) tau=-1/B(2) plot(t,sigma,'o'), x=0:0.01:2.5;
y=sigma0*exp(-x/tau); hold on plot(x,y) xlabel('t'),
ylabel('\\sigma') grid on 结果：
sigma0 = 5.0002 tau =
0.9999
?t?6543?21000.51t1.522.5
2. 已知某物质的饱和蒸汽压和温度有关，并已测得下面一组数据。
温度 T/oC 压力 P/mmHg 利用模型：
0.534 0.101 0.822 0.348 1.110 0.647 1.398 0.841 1.686 0.925 1.972 0.958 2.342 0.978 2.712 0.988 3.206 0.995 P?a1?b1lnT?c1T3
计算方程中的待求参数，绘出P-T曲线。(c2的范围是-100 ~ 100)。 程序：
T=[0.534,0.822,1.110,1.398,1.686,1.972,2.342,2.712,3.206]'; P=[0.101,0.348,0.647,0.841,0.925,0.958,0.978,0.988,0.995]'; Y=P
X=ones(9,3); X(:,2)=log(T); X(:,3)=T.^3; B=X\\Y a1=B(1) b1=B(2) c1=B(3) plot(T,P) hold on xlabel('T') ylabel('P') grid on 结果：
a1 = 0.5680 b1 =0.7582 c1 = -0.0149
1.210.80.60.40.20-0.2-0.4-0.6-0.800.511.5T22.533.5P
把程序及运行结果写在下面：
3. 用非线性最小二乘法计算第1题中的参数??和?。
把程序及运行结果写在下面：
t=[0.2,0.6,1.0,1.4,1.8,2.2];
sigma=[4.094,2.744,1.839,1.233,0.826,0.554]; f=inline('a(1)*exp(-x/a(2))','a','x') a=lsqcurvefit(f,[1,1],t,sigma); sigma0=a(1) tau=a(2) 结果：
sigma0 = 5.0006 tau =
0.9998
4. 自由练习内容可拷贝在下面