计算物理MATLAB学习记录

第一次用MATLAB,折腾来折腾去码了三天终于码完了作业题目，放在这里留个纪念

实验题目1

work 1.1

Homework 1.1 生成矩阵
$\left[\begin{array}{ccccc}0& 1& 1& 1& 1\\ 1& 0& 1& 1& 1\\ 1& 1& 0& 1& 1\\ 1& 1& 1& 0& 1\\ 1& 1& 1& 1& 0\end{array}\right]$

A=ones(5,5);
B=eye(5,5);
C=A-B;
disp(C);

work 1.2

在$[0,2\pi ]$ 范围内绘制出$sin{x}^2$和$co{s}^{2}x$曲线，要求：步长为0.01，两条曲线同时在一个图中显示，分别用蓝色实线和红色点虚线表示，横纵坐标等长。

x=0:0.01:2*pi;
y1=sin(x.^2);
y2=cos(x).^2;
plot(x,y1,'b-',x,y2,'r:');

work 1.3

生成1行1000列的正态分布随机数组$x$，其中第5、15和20项置为100。求 $x$ 的平均值, 方差, 标准偏差, 并用拉依达方法剔除所测数据的错误值，结果$y$用直方图表示(默认)

%生成矩阵
x=randn(1,1000);
x(5)=100;
x(15)=100;
x(20)=100;
%计算均值、方差、标准差
Avg_x=mean(x);
Var_x=var(x);
Std_x=std(x);
%剔除错误值
y=[];
for i=1:length(x)
   d=abs(x(i)-Avg_x);
   if d<=3*Std_x
       y=[y,x(i)];
   end
end
%直方图
hist(y);

work 1.4

用肖维勒方法剔除第一题中$x$的错误值，结果用直方图表示

%生成矩阵,计算均值、方差、标准差
x=randn(1,1000);
x(5)=100;
x(15)=100;
x(20)=100;
Avg_x=mean(x);
Var_x=var(x);
Std_x=std(x);
%剔除错误值
n=length(x);
wn=1+0.4*log(n);
y=[];
for i=1:n
    d=abs(x(i)-Avg_x);
    if d<=wn*Std_x
        y=[y,x(i)];
    end
end
%绘制直方图
hist(y);

work 1.5

用差分法剔除$x=-0.9589,-0.7568,-0.5440,-5.0000,$