Geum.ru

Исследование переходных процессов в механической системе

Курсовой проект - Компьютеры, программирование

Другие курсовые по предмету Компьютеры, программирование

-0.024 -0.482 -0.062

.345 -0.053 -0.024 -0.366 -0.011

.352 -0.056 -0.024 -0.247 0.040

.358 -0.057 -0.024 -0.128 0.087

.365 -0.058 -0.023 -0.009 0.132

.372 -0.058 -0.022 0.105 0.173

.379 -0.058 -0.021 0.213 0.210

.385 -0.056 -0.020 0.314 0.241

.392 -0.054 -0.018 0.405 0.266

.399 -0.051 -0.016 0.485 0.285

.406 -0.048 -0.014 0.553 0.297

.413 -0.044 -0.012 0.607 0.303

.419 -0.040 -0.010 0.647 0.303

.426 -0.036 -0.008 0.673 0.297

.433 -0.031 -0.006 0.685 0.284

.440 -0.027 -0.004 0.682 0.267

.446 -0.022 -0.002 0.666 0.244

.453 -0.017 -0.001 0.637 0.218

.460 -0.013 0.001 0.596 0.188

.467 -0.009 0.002 0.546 0.155

.473 -0.005 0.003 0.487 0.121

.480 -0.002 0.004 0.421 0.087

.487 0.001 0.004 0.350 0.052

.494 0.003 0.005 0.277 0.019

.500 0.005 0.005 0.203 -0.013

.507 0.006 0.005 0.129 -0.042

.514 0.007 0.004 0.059 -0.067

.521 0.008 0.004 -0.006 -0.088

.527 0.008 0.003 -0.065 -0.105

.534 0.007 0.003 -0.115 -0.117

.541 0.006 0.002 -0.157 -0.124

.548 0.005 0.001 -0.188 -0.125

.555 0.004 0.000 -0.208 -0.121

.561 0.003 -0.001 -0.216 -0.112

.568 0.001 -0.001 -0.213 -0.098

.575 -0.000 -0.002 -0.199 -0.079

.582 -0.002 -0.003 -0.174 -0.057

.588 -0.003 -0.003 -0.138 -0.031

.595 -0.004 -0.003 -0.094 -0.002

.602 -0.004 -0.003 -0.042 0.028

.609 -0.005 -0.003 0.017 0.060

.615 -0.005 -0.003 0.080 0.091

.622 -0.004 -0.002 0.146 0.122

.629 -0.003 -0.001 0.213 0.151

.636 -0.002 -0.000 0.280 0.178

.642 0.000 0.001 0.343 0.201

.649 0.003 0.002 0.402 0.221

.656 0.005 0.004 0.454 0.236

.663 0.008 0.006 0.498 0.246

.669 0.012 0.007 0.532 0.250

.676 0.015 0.009 0.556 0.249

.683 0.019 0.011 0.535 0.242

.690 0.023 0.012 0.511 0.216

.697 0.026 0.014 0.485 0.177

.703 0.030 0.015 0.456 0.135

.710 0.033 0.016 0.424 0.096

.717 0.035 0.016 0.388 0.067

.724 0.038 0.017 0.350 0.049

.730 0.040 0.017 0.310 0.040

.737 0.043 0.018 0.268 0.038

.744 0.044 0.018 0.225 0.037

.751 0.046 0.018 0.183 0.034

.757 0.047 0.018 0.140 0.025

.764 0.048 0.018 0.097 0.011

.771 0.049 0.019 0.054 -0.009

.778 0.049 0.018 0.011 -0.033

.784 0.049 0.018 -0.031 -0.057

.791 0.049 0.018 -0.073 -0.079

.798 0.048 0.017 -0.115 -0.097

.805 0.048 0.017 -0.155 -0.112

.811 0.047 0.016 -0.195 -0.122

.818 0.045 0.015 -0.233 -0.130

.825 0.044 0.014 -0.269 -0.137

.832 0.042 0.013 -0.302 -0.143

.839 0.040 0.012 -0.334 -0.150

.845 0.038 0.011 -0.363 -0.157

.852 0.035 0.010 -0.390 -0.166

.859 0.033 0.009 -0.414 -0.174

.866 0.030 0.008 -0.435 -0.181

.872 0.027 0.007 -0.454 -0.188

.879 0.024 0.005 -0.470 -0.192

.886 0.021 0.004 -0.483 -0.194

.893 0.017 0.003 -0.493 -0.193

.899 0.014 0.002 -0.500 -0.191

.906 0.011 0.000 -0.504 -0.188

.913 0.007 -0.001 -0.505 -0.184

.920 0.004 -0.002 -0.503 -0.179

.926 0.000 -0.004 -0.498 -0.174

.933 -0.003 -0.005 -0.490 -0.169

.940 -0.006 -0.006 -0.479 -0.162

.947 -0.010 -0.007 -0.465 -0.155

.953 -0.013 -0.008 -0.449 -0.147

.960 -0.016 -0.009 -0.431 -0.138

.967 -0.019 -0.010 -0.410 -0.128

.974 -0.021 -0.011 -0.387 -0.118

.981 -0.024 -0.012 -0.362 -0.106

.987 -0.026 -0.012 -0.335 -0.094

.994 -0.029 -0.013 -0.307 -0.081

.001 -0.031 -0.013 -0.276 -0.068

.008 -0.033 -0.014 -0.245 -0.055

.014 -0.034 -0.014 -0.212 -0.042

.021 -0.036 -0.015 -0.179 -0.029

.028 -0.037 -0.015 -0.145 -0.016

.035 -0.038 -0.015 -0.110 -0.003

.041 -0.039 -0.015 -0.075 0.010

.048 -0.039 -0.015 -0.040 0.022

.055 -0.039 -0.015 -0.005 0.035

.062 -0.040 -0.014 0.029 0.047

.068 -0.039 -0.014 0.063 0.059

.075 -0.039 -0.014 0.096 0.071

.082 -0.038 -0.013 0.129 0.082

.089 -0.037 -0.013 0.160 0.092

.096 -0.036 -0.012 0.190 0.101

.102 -0.035 -0.011 0.218 0.110

.109 -0.034 -0.011 0.245 0.118

.116 -0.032 -0.010 0.270 0.126

.123 -0.030 -0.009 0.294 0.132

.129 -0.028 -0.008 0.315 0.138

.136 -0.026 -0.007 0.334 0.142

.143 -0.024 -0.006 0.351 0.146

.150 -0.021 -0.005 0.366 0.149

.156 -0.019 -0.004 0.379 0.152

.163 -0.016 -0.003 0.389 0.153

.170 -0.014 -0.002 0.397 0.153

.177 -0.011 -0.001 0.402 0.152

.183 -0.008 -0.000 0.405 0.151

.190 -0.006 0.001 0.405 0.148

.197 -0.003 0.002 0.403 0.145

.204 -0.000 0.003 0.399 0.141

.210 0.003 0.004 0.393 0.136

.217 0.005 0.005 0.384 0.131

.224 0.008 0.006 0.373 0.124

.231 0.010 0.007 0.360 0.117

.238 0.013 0.007 0.345 0.110

.244 0.015 0.008 0.328 0.102

.251 0.017 0.009 0.309 0.093

.258 0.020 0.009 0.289 0.084

.265 0.021 0.010 0.267 0.074

.271 0.023 0.010 0.244 0.064

.278 0.025 0.011 0.220 0.054

.285 0.026 0.011 0.194 0.044

.292 0.028 0.012 0.168 0.033

.298 0.029 0.012 0.141 0.023

.305 0.030 0.012 0.114 0.012

.312 0.031 0.012 0.086 0.002

.319 0.031 0.012 0.058 -0.009

.325 0.032 0.012 0.030 -0.019

.332 0.032 0.012 0.002 -0.029

.339 0.032 0.012 -0.026 -0.039

.346 0.032 0.011 -0.053 -0.049

maxy1 maxy2 maxV1 maxV2

.062 0.025 1.076 0.446

 

6.4 Программа Модифицированный метод Эйлера program ModEiler

 

typeodnmas = array [1..201] of real;

T2thenT0:=3*T1T0:=3*T2;:=T0/n;:=1;[i]:=0;[i]:=y10;y2[i]:=y20;[i]:=V10;V2[i]:=V20;(f2,ty1y2V1V2);[i+1]:=t[i]+ht;[i+1]:=y1[i]+ht*(V1[i]+V1[i]+ht*fv1(y1[i],y2[i],V1[i],V2[i],t[i]))/2;[i+1]:=V1[i]+ht*(fv1(y1[i],y2[i],V1[i],V2[i],t[i])+fv1(y1[i]+">var ht,y10,y20,V10,V20,C1,C2,k,m1,m2,,g,d1,d2,T1,T2,T0,di,maxy1,maxy2,maxV1,maxV2:real;,y1,y2,V1,V2:odnmas;,n:integer;,f2:text;fv1(y1,y2,V1,V2,t:real):real;tT2 then T0:=3*T1T0:=3*T2;:=T0/n;:=1;[i]:=0;[i]:=y10;y2[i]:=y20;[i]:=V10;V2[i]:=V20;(f2, t y1 y2 V1 V2);[i+1]:=t[i]+ht;[i+1]:=y1[i]+ht*(V1[i]+V1[i]+ht*fv1(y1[i],y2[i],V1[i],V2[i], t[i]))/2;[i+1]:=V1[i]+ht*(fv1(y1[i],y2[i],V1[i],V2[i],t[i])+fv1(y1[i]+

+ht*V1[i],y2[i]+ht*V2[i],V1[i]+ht*fv1(y1[i],y2[i],V1[i],V2[i],t[i]),V2[i]+ht*fv2(y1[i],y2[i],V1[i],V2[i],t[i]),t[i]+ht))/2;[i+1]:=y2[i]+ht*(V2[i]+V2[i]+ht*fv2(y1[i],y2[i],V1[i],V2[i],t[i]))/2;[i+1]:=V2[i]+ht*(fv2(y1[i],y2[i],V1[i],V2[i],t[i])+fv2(y1[i]+ht*V1[i],y2[i]+ht*V2[i],V1[i]+ht*fv1(y1[i],y2[i],V1[i],V2[i],t[i]),V2[i]+ht*fv2(y1[i],y2[i],V1[i],V2[i],t[i]),t[i]+ht))/2;

{di:=i/10;(frac(di)=0)or(frac(di)=0.1) then}(f2,t[i]:3:3,y1[i]:8:3,y2[i]:8:3,V1[i]:8:3,V2[i]:8:3);:=i+1;i>n+1;

:=y1[1];:=y2[1];:=V1[1];:=V2[1];

i:=2 to n doAbs(y1[i])>maxy1 then:=Abs(y1[i]);Abs(y2[i])>maxy2 then:=Abs(y2[i]);Abs(V1[i])>maxV1 then:=Abs(V1[i]);Abs(V2[i])>maxV2 then:=Abs(V2[i]);;(f2);(f2, maxy1 ,maxy2 ,maxV1 ,maxV2 );(f2,maxy1:8:3,maxy2:8:3,maxV1:8:3,maxV2:8:3);(f1);(f2);.

 

6.5 Файл вывода результата

y1 y2 V1 V2

.000 0.000 0.000 0.000 0.000

.007 0.000 0.000 0.005 0.000

.014 0.000 0.000 0.020 0.004

.020 0.000 0.000 0.043 0.014

.027 0.001 0.000 0.075 0.030

.034 0.001 0.000 0.113 0.052

.041 0.002 0.001 0.158 0.079

.047 0.003 0.002 0.208 0.107

.054 0.005 0.002 0.260 0.136

.061 0.007 0.003 0.314 0.164

.068 0.009 0.005 0.368 0.189

.074 0.012 0.006 0.419 0.211

.081 0.015 0.007 0.466 0.229

.088 0.018 0.009 0.507 0.242

.095 0.022 0.011 0.541 0.250

.101 0.026 0.012 0.566 0.253

.108 0.029 0.014 0.581 0.250

.115 0.033 0.016 0.585 0.242

.122 0.037 0.017 0.577 0.228

.128 0.041 0.019 0.556 0.209

.135 0.045 0.020 0.522 0.184

.142 0.048 0.021 0.477 0.155

.149 0.051 0.022 0.419 0.120

.156 0.054 0.023 0.349 0.082

.162 0.056 0.023 0.269 0.040

.169 0.058 0.024 0.180 -0.005

.176 0.058 0.023 0.084 -0.051

.183 0.059 0.023 -0.019 -0.099

.189 0.058 0.022 -0.126 -0.146

.196 0.057 0.021 -0.235 -0.192

.203 0.055 0.019 -0.344 -0.237

.210 0.052 0.018 -0.451 -0.278

.216 0.049 0.016 -0.554 -0.316

.223 0.045 0.013 -0.650 -0.349

.230 0.040 0.011 -0.738 -0.377

.237 0.035 0.008 -0.817 -0.399

.243 0.029 0.006 -0.883 -0.415

.250 0.023 0.003 -0.937 -0.423

.257 0.016 -0.000 -0.976 -0.425

.264 0.010 -0.003 -1.001 -0.420

.270 0.003 -0.006 -1.009 -0.408

.277 -0.004 -0.009 -1.002 -0.390

.284 -0.011 -0.011 -0.979 -0.364

.291 -0.017 -0.014 -0.941 -0.333

.298 -0.023 -0.016 -0.887 -0.297

.304 -0.029 -0.018 -0.821 -0.256

.311 -0.034 -0.019 -0.741 -0.211

.318 -0.039 -0.020 -0.651 -0.164

.325 -0.043 -0.021 -