计算物理作业.docVIP

计算物理

姓名：杨永富

学号：109040669

专业：物理电子学

任课教师：沈中华

作业一：

x

0

1

2

3

4

5

6

7

y

7.82

7.93

7.98

7.59

7.92

7.91

7.80

7.71

试用最小二乘法寻求经验公式，并拟合以上数据

程序：

clearall,closeall,clc

%===============原始数据========================

X=[01234567];

Y=[7.827.937.987.997.927.917.807.71];

%===============================================

plot(X,Y,k*);%绘制原始数据散点图

A=polyfit(X,Y,2);%多项式拟合

x1=0:0.1:8;

z=polyval(A,x1);

holdon

plot(x1,z,r);%绘制拟合曲线

xlabel(x);

ylabel(y);

gridon

legend(试验数据散点图,拟合曲线);%图形标注

title(y-x曲线);

gtext(109040669杨永富)

vpa(poly2sym(A),3)%输出多项式的表达式，系数取3位有效数字

计算结果：

因此

在某次实验中，需要观察水分的渗透速度，测得时间t与水的重量W的数据如下：

t/s

1

2

4

8

16

32

64

W/g

4.22

4.02

3.85

4.59

3.44

3.02

2.59

设t与W的关系w=Ats试用最小二乘法确定A,s

程序：

clearall,closeall

t=[124163264];

W=[4.224.023.853.443.022.59];

s_lnt=log(t);s_lnW=log(W);plot(s_lnt,s_lnW,k+),holdon;

A=polyfit(s_lnt,s_lnW,1);%多项式拟合

z=polyval(A,s_lnt);plot(s_lnt,z,r);

title(lnW-lnt曲线);xlabel(lnt);ylabel(lnW);

gridon,gtext(109040669杨永富)

figure,plot(t,W,r*),holdon;

s=A(1);A=exp(A(2));t=0:0.5:70;W=A*t.^s;

plot(t,W),gridon;

legend(试验数据散点图,拟合曲线);%图形标注

xlabel(t/s);ylabel(W/g);

title(W-t曲线);gtext(109040669杨永富)

因此

作业二：

复化梯形公式计算定积分：

要求：递交算法说明过程、源程序及与实际结果的比拟

算法说明：

把积分区间分割成n个小区间，在每个小区间上用梯形积分公式:

再将n个小区间上的数值积分累加:

有误差估计：

把积分区间分割成2n个小区间，有误差估计：

即计算中构造序列：。

每次计算的积分值均要与前一次的积分值进行比拟，直到，停止计算，并将作为计算结果。

程序：

%复化梯形公式

clearall

closeall

clc

%**************************

%======在此修改积分限=====

I_low=0;

I_high=pi;

%======在此修改积分限=====

%**************************

epsilong=input(Pleaseinputtheerroryoucouldaccept:\n);

fprintf(\n\n);

n=2;

T_back=0;

T_front=0;

delta=5*epsilong;

tic

whiledelta=3*epsilong;

Cigema_f=0;

T_front=T_back;

n=2*n;

h=I_high/n;

x=0:h:I_high;

%**************************

%======在此修改被积函数=====

f=sin(x);

%======在此修改被积函数=====

%**************************

Cigema_f=sum(f(2:n-1));

T_back=h*(0.5*f(1)+Cigema_f+0.5*f(n))