Anonymous     Text only     09 Dec 2012    
  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
3.2--------------------------------------------------------
// 1 ????
// 2
a=0.8;
b=1.6;
h=0.02;
x=a:h:b;
y=(2*x^2+1)^(-0.5);
inttrap(x,y)
// 2a
integrate('(2*x^2+1)^(-0.5)','x',0.8,1.6)
// 3
function yd=f(t,y)
yd=0.133*(t^2+sin(2*t)+0.872*y)
endfunction
y0=0.25;
t0=0.2;
t=0.2:0.1:1.2;
y=ode(y0,t0,t,f)
plot(t,y)
4.1----------------------------------------------------------
function [P]=G(c,z)
P=z(2)-c(3)-c(2)*z(1)^2-c(1)*z(1)^3
endfunction
x=[0 1 1.5 2 2.5 3 3.5 4 4.5 5];
y=[12 10.1 11.58 17.4 30.68 53.6 87.78 136.9 202.5 287];
z=[x;y];
c=[0;0;0];
// Расчет суммы квадратов отклонений
// Решение
[a,err]=datafit(G,z,c)
// Графическое решение
plot2d(x,y,-4)
t=0:0.01:5;
Ptc=a(3)+a(2)*t^2+a(1)*t^3;
plot2d(t,Ptc);
// Коэффиценты регрессии
a=regress(x,y)
// Коэффицент кореляции
r=sum((x-mean(x)).*(y-mean(y)))/...
sqrt(sum((x-mean(x))^2)*sum((y-mean(y))^2))
4.2 ------------------------------------------------
x=[0.43 0.48 0.55 0.62 0.7 0.75];
y=[1.63597 1.73234 1.87686 2.03345 2.22846 2.35973];
plot2d(x,y,-4);
// График экспериментальных данных
koeff=splin(x,y);
X=[0.702 0.512 0.608];
// Значение функции в заданных точках
Y=interp(X,x,y,koeff)
// Нанесение точек на график
plot2d(X,Y,-3);
// Построение кубического сплайна
t=0.43:0.01:0.75;
ptd=interp(t,x,y,koeff);
plot2d(t,ptd);
4.3--------------------------------------------
function [f,g,ind]=fi(x,ind)
f=a*x.^4+b*x.^3+c*x.^2+d*x+e
g=4*a*x.^3+3*b*x.^2+2*c*x+d
endfunction
x=-5:0.1:1;
a=2;
b=3;
c=-15;
d=-5;
e=25;
y=a*x.^4+b*x.^3+c*x.^2+d*x+e;
plot(x,y);
xgrid();
y0=-3;
[ymin,xmin]=optim(fi,y0)
5.1-----------------------------------
// 1
function yd=f(t,y)
yd=1-sin(2*t+y)
endfunction
y0=0;
t0=0;
t=0:0.2:5;
y=ode(y0,t0,t,f)
plot(t,y)
5.2----------------------------------------
function dy=syst(t,y)
dy=zeros(2,1);
dy(1)=cos(y(1)+2*y(2))+2;
dy(2)=2/(t+2*y(1)^2)+t+1;
endfunction
y0=[1;0.05];t0=0;t=0:0.1:10;y=ode(y0,t0,t,syst);
plot(t,y)