Методы обработки экспериментальных данных
Автор: albina002 • Октябрь 28, 2021 • Реферат • 1,191 Слов (5 Страниц) • 284 Просмотры
Әл-Фараби атындағы Қазақ Ұлттық Университеті
[pic 1]
СӨЖ
Тақырыбы: МЕТОДЫ ОБРАБОТКИ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ.
Орындаған: Кадерова Альбина СиС349
Тексерген: Алдиярова А.Н
Алматы 2021
Задание 1
x=[1.20;1.57;1.94;2.31;2.68;3.05;3.42;3.79]
x =
1.2000
1.5700
1.9400
2.3100
2.6800
3.0500
3.4200
3.7900
y=[2.56;1.06;1.58;1.25;0.91;0.66;0.38;0.21]
y =
2.5600
1.0600
1.5800
1.2500
0.9100
0.6600
0.3800
0.2100
x2=[x.^2];
xy=[x.*y];
Mx=1/6*sum(x);
My=1/6*sum(y);
Mx2=1/6*sum(x2);
Mxy=1/6*sum(xy);
M(1,1)=Mx2;
M(1,2)=Mx;
M(2,1)=Mx;
M(2,2)=1;
d(1,1)=Mxy;
d(2,1)=My;
Coeff=M^-1*d
Coeff =
1.0513
-2.0622
a=(Mxy-Mx*My)/(Mx2-Mx^2)
a =
1.0513
b=(Mx2*My-Mx*Mxy)/(Mx2-Mx^2)
b =
-2.0622
y1=a*x+b;
e2=(y-y1).^2;
S=sum(e2)
S =
21.3291
plot(x,y1,x,y,'*')
[pic 2]
b)
x=[1.20;1.57;1.94;2.31;2.68;3.05;3.42;3.79]
x =
1.2000
1.5700
1.9400
2.3100
2.6800
3.0500
3.4200
3.7900
u=(1./x);
y=[2.56;1.06;1.58;1.25;0.91;0.66;0.38;0.21]
y =
2.5600
1.0600
1.5800
1.2500
0.9100
0.6600
0.3800
0.2100
u2=[u.^2];
uy=[u.*y];
Mu=1/6*sum(u);
My=1/6*sum(y);
Mu2=1/6*sum(u2);
Muy=1/6*sum(uy);
M(1,1)=Mu2;
M(1,2)=Mu;
M(2,1)=Mu;
M(2,2)=1;
d(1,1)=Muy;
d(2,1)=My;
Coeff=M^-1*d
Coeff =
1.3080
0.6337
a=(Muy-Mu*My)/(Mu2-Mu^2)
a =
1.3080
b=(Mu2*My-Mu*Muy)/(Mu2-Mu^2)
b =
0.6337
y2=a./x+b;
e2=(y-y2).^2;
S=sum(e2)
S =
2.1441
plot(x,y2,x,y,'*')
[pic 3]
Задание 3 А
>> x=[-4.38;-3.84;-3.23;-2.76;-2.22;-1.67;-1.13;-0.60];% задание исходных данных
>> y=[2.25;2.83;3.44;4.51;5.29;6.55;8.01;10.04];
>> x2=[x.^2];
>> xy=[x.*y];
>> Mx=1/6*sum(x); % вычисление элементов матриц М и d
>> My=1/6*sum(y);
>> Mx2=1/6*sum(x2);
>> Mxy=1/6*sum(xy);
>> M(1,1)=Mx2; % задание матрицы М
The desktop configuration was not saved successfully
>> M(1,2)=Mx;
>> M(2,1)=Mx;
>> M(2,2)=1;
>> d(1,1)=Mxy; % задание матрицы d
>> d(2,1)=My;
>> Coeff=M^-1*d % решение системы линейных уравнений (11.1)
Coeff =
-14.2293
-39.8744
>> a=(Mxy-Mx*My)/(Mx2-Mx^2) % вычисление коэффициентов с помощью формул (11.3)
a =
-14.2293
ans =
11.3000
>> b=(Mx2*My-Mx*Mxy)/(Mx2-Mx^2)
b =
-39.8744
% вычисление суммы квадратов отклонений
>> y1=a*x+b;
>> e2=(y-y1).^2;
>> S=sum(e2)
S =
4.0051e+03
>> % построение графика полученной функции и исходных данных
>> plot(x,y1,x,y,'*')
[pic 4]
Б)
>> x=[-4.38;-3.84;-3.23;-2.76;-2.22;-1.67;-1.13;-0.60];
>> u=(1./x); % введение замены переменной х
>> y=[2.25;2.83;3.44;4.51;5.29;6.55;8.01;10.04];
>> u2=[u.^2];
>> uy=[u.*y];
>> Mu=1/6*sum(u);
>> My=1/6*sum(y);
>> Mu2=1/6*sum(u2);
>> Muy=1/6*sum(uy);
>> M(1,1)=Mu2;
>> M(1,2)=Mu;
>> M(2,1)=Mu;
>> M(2,2)=1;
>> d(1,1)=Muy;
>> d(2,1)=My;
>> Coeff=M^-1*d
Coeff =
-0.0230
7.1351
>> a=(Muy-Mu*My)/(Mu2-Mu^2)
a =
-0.0230
>> b=(Mu2*My-Mu*Muy)/(Mu2-Mu^2)
b =
7.1351
>> y2=a./x+b;
>> e2=(y-y2).^2;
>> S=sum(e2)
S =
75.8907
>> plot(x,y2,x,y,'*')
[pic 5]
Задание 5
x=[5.89;3.84;6.19;9.22;7.87;6.29;4.43;8.91]
x =
5.8900
3.8400
6.1900
9.2200
7.8700
6.2900
4.4300
8.9100
...