Численные методы интегрирования
Материал из Department of Theoretical and Applied Mechanics
Версия от 19:51, 4 июня 2019; 89.223.47.218 (обсуждение)
Описание
Визуализация различных численных методов
Исполнитель: Дурнев Андрей
Группа 13632/1 Кафедра Теоретической механики
Файл:[[1]]
Код программы
Текст программы на языке JavaScript:
1 //Метод Монте-Карло
2
3 window.addEventListener("load", program_code, false);
4 function program_code() {
5
6 var ctx = canvas_example.getContext("2d");
7 var w = canvas_example.width;
8 var h = canvas_example.height;
9
10 ctx.translate(w/2, h/2);
11
12 var x = [];
13 var y = [];
14 var min = 1000;
15 var max = -1000;
16 var n = 0;
17 var points;
18 var cur_points;
19 var ok = 0;
20 var a;
21 var b ;
22 var scalex, scaley;
23 var dx;
24 var dy;
25 var func;
26 var time = 5000;
27 var dt;
28 var flag = false;
29
30
31 var result;
32
33
34
35 //функция
36 function f(x){
37 return (eval(func));
38 }
39
40 function search(){
41
42 func = document.getElementById("fun").value;
43
44 for (var i = a; i < b; i+=0.001){
45
46 //максимум
47 if (f(i) > max){
48 max = f(i);
49 }
50
51 //минимум
52 if (f(i) < min){
53 min = f(i);
54 }
55 }
56
57 //масштаб
58 scaley = 100/(Math.abs(max) + Math.abs(min));
59
60 if (Math.abs(b) >= Math.abs(a)){
61 scalex = 250/b * 0.9;
62 }
63
64 if (Math.abs(b) < Math.abs(a)){
65 scalex = 250/Math.abs(a) * 0.9;
66 }
67
68 //сдвиг
69 dx = (b + a)/2;
70 dy = (max + min)/2;
71
72 //точки
73 //points = 100;
74 }
75
76 function calculation(){
77
78 cur_points++;
79
80 //положительная функция
81 if (min * max >= 0 && min >= 0 && n < points){
82 x[n] = Math.random() * (b - a) + a;
83 y[n] = Math.random() * (max - 0) + 0;
84
85 if (y[n] < f(x[n])){
86 ok++;
87 }
88
89
90 result = (ok / cur_points) * (b - a) * max;
91 n++;
92
93 }
94
95 //отрицательная функция
96 if (min * max >= 0 && max <= 0 && n < points){
97 x[n] = Math.random() * (b - a) + a;
98 y[n] = Math.random() * (min - 0) + 0;
99
100 if (y[n] > f(x[n])){
101 ok++;
102 }
103
104 result = -(ok/cur_points) * (b - a) * (-min);
105 n++;
106 }
107
108 //знакопеременные функции
109 if (min * max < 0 && n < points){
110
111 x[n] = Math.random() * (b - a) + a;
112 y[n] = Math.random() * (max - min) + min;
113
114
115 if (y[n] > 0 && f(x[n]) > 0 && y[n] < f(x[n])){
116 ok++;
117 }
118
119 if (y[n] < 0 && f(x[n]) < 0 && y[n] > f(x[n])){
120 ok--;
121 }
122
123 result = (ok/cur_points) * (b - a) * (max - min);
124 n++;
125 }
126 }
127
128 function draw(){
129 //оси
130 ctx.beginPath();
131 ctx.moveTo(0 - dx * scalex, -h/2);
132 ctx.lineTo(0 - dx * scalex, h/2);
133 ctx.moveTo(w/2 , dy * scaley);
134 ctx.lineTo(-w/2 , dy * scaley);
135 ctx.stroke();
136
137 ctx.font = '10px Arial';
138 ctx.strokeText('1', (1 - dx) * scalex, (-0.1 + dy) * scaley);
139 ctx.strokeText('1', (0.1 -dx) * scalex, (-1 + dy) * scaley);
140
141
142 //кривая
143 ctx.beginPath();
144 for (var i = a; i < b; i+=0.01){
145 ctx.moveTo((i - dx) * scalex, (-f(i) + dy) * scaley);
146 ctx.lineTo((i+0.01 - dx) * scalex, (-f(i+0.01) + dy) * scaley);
147 }
148
149 ctx.stroke();
150
151
152 //точки
153 ctx.fillStyle = 'red';
154 ctx.beginPath();
155 ctx.arc((x[n-1] - dx) * scalex, (-y[n-1] + dy) * scaley, 1, 0, 2 * Math.PI);
156 ctx.fill();
157
158 }
159
160 function control(){
161 calculation();
162 draw();
163 result1.innerHTML = Math.round(result * 100)/100;
164 }
165
166 knopka.onclick = function(){
167
168 if (flag){
169 location.reload();
170 }
171
172 flag = true;
173 knopka.value = 'Очистить';
174 cur_points = 0;
175 points = parseFloat(document.getElementById("qual").value);
176 dt = time/points;
177 console.log(dt)
178 a = parseFloat(document.getElementById("a1").value);
179 b = parseFloat(document.getElementById("b1").value);
180 search();
181 var name = setInterval(control, dt);
182 }
183
184
185 if (n == points){
186 clearInterval(name);
187 console.log('a');
188 }
189 }
190
191 //Метод трапеций
192 window.addEventListener("load", program_code, false);
193 function program_code() {
194
195 var ctx = canvas_example.getContext("2d");
196 var w = canvas_example.width;
197 var h = canvas_example.height;
198
199 var flag = false;
200 var a;
201 var b;
202 var func;
203 var max = -100000;
204 var min = 100000;
205 var scale, scale1, scale2;
206 var x = [];
207 var y = [];
208 var result = 0;
209
210 ctx.translate(w/2, h/2);
211
212 //функция
213 function f(x){
214 return (eval(func));
215 }
216
217 function input(){
218
219 a = parseFloat(document.getElementById("a1").value);
220 b = parseFloat(document.getElementById("b1").value);
221 func = document.getElementById("fun").value;
222 number = parseFloat(document.getElementById('qual').value);
223
224 //минимум и максимум
225 for (var i = a; i < b; i+=0.001){
226
227 //максимум
228 if (f(i) > max){
229 max = f(i);
230 }
231
232 //минимум
233 if (f(i) < min){
234 min = f(i);
235 }
236 }
237
238 //масштаб
239 scaley = 100/(Math.abs(max) + Math.abs(min));
240
241 if (Math.abs(b) >= Math.abs(a)){
242 scalex = 250/b * 0.9;
243 }
244
245 if (Math.abs(b) < Math.abs(a)){
246 scalex = 250/Math.abs(a) * 0.9;
247 }
248
249 //сдвиг
250 dx = (b + a)/2;
251 dy = (max + min)/2;
252 dt = (b - a)/number;
253
254 //разбиение
255 for (var i = a; i < b + 1/2 * dt; i += dt){
256 x.push(i);
257 y.push(f(i));
258 }
259 }
260
261 function calculation(){
262 for (var i = 0; i < number; i++){
263
264 if (f(x[i]) >= 0 && f(x[i+1]) >= 0){
265 result += (x[i + 1] - x[i]) * (Math.abs(y[i]) + Math.abs(y[i + 1])) * 0.5;
266 }
267
268 if (f(x[i]) < 0 && f(x[i+1]) < 0){
269 result -= (x[i + 1] - x[i]) * (Math.abs(y[i]) + Math.abs(y[i + 1])) * 0.5;
270 }
271
272 else{
273 result += 0;
274 }
275
276 }
277 }
278
279 function draw(){
280
281 //оси
282 ctx.beginPath();
283 ctx.moveTo(0 - dx * scalex, -h/2);
284 ctx.lineTo(0 - dx * scalex, h/2);
285 ctx.moveTo(w/2 , 0 + dy * scaley);
286 ctx.lineTo(-w/2 , 0 + dy * scaley);
287 ctx.stroke();
288
289 ctx.font = '10px Arial';
290 ctx.strokeText('1', (1 - dx) * scalex, (-0.1 + dy) * scaley);
291 ctx.strokeText('1', (0.1 -dx) * scalex, (-1 + dy) * scaley);
292
293 //кривая
294 ctx.beginPath();
295 for (var i = a; i < b; i+=0.01){
296 ctx.moveTo((i - dx) * scalex, (-f(i) + dy) * scaley);
297 ctx.lineTo((i + 0.01 - dx) * scalex, (-f(i + 0.01) + dy) * scaley);
298 }
299 ctx.stroke();
300
301
302 ctx.strokeStyle = 'red';
303
304 //разбиение
305 ctx.beginPath();
306 for (var i = 0; i < number + 1; i++){
307 ctx.moveTo((x[i] - dx) * scalex, dy * scaley);
308 ctx.lineTo((x[i] - dx) * scalex, (-y[i] + dy) * scaley)
309 }
310 ctx.stroke();
311
312 ctx.beginPath();
313 for(var i = 0; i < number; i++){
314 ctx.moveTo((x[i] - dx) * scalex, (-y[i] + dy) * scaley);
315 ctx.lineTo((x[i+1] - dx) * scalex, (-y[i+1] + dy) * scaley);
316 }
317 ctx.stroke();
318
319 }
320
321
322 knopka.onclick = function(){
323
324 if (flag){
325 location.reload();
326 }
327
328 flag = true;
329 knopka.value = 'Очистить';
330 input();
331 calculation();
332 result1.innerHTML = Math.round(result * 100)/100;
333 draw();
334 }
335 }
336
337
338 //Метод прямоугольников
339 window.addEventListener("load", program_code, false);
340 function program_code() {
341
342 var ctx = canvas_example.getContext("2d");
343 var w = canvas_example.width;
344 var h = canvas_example.height;
345
346 var flag = false;
347 var a;
348 var b;
349 var func;
350 var max = -100000;
351 var min = 100000;
352 var scale, scale1, scale2;
353 var x = [];
354 var y = [];
355 var result = 0;
356 var dt;
357 var dx = 0;
358 var dy = 0;
359
360 ctx.translate(w/2, h/2);
361
362 //функция
363 function f(x){
364 return (eval(func));
365 }
366
367 function input(){
368
369 a = parseFloat(document.getElementById("a1").value);
370 b = parseFloat(document.getElementById("b1").value);
371 func = document.getElementById("fun").value;
372 var number = parseFloat(document.getElementById('qual').value);
373 dx = (b - a)/number;
374
375 //минимум и максимум
376 for (var i = a; i < b; i+=0.001){
377
378 //максимум
379 if (f(i) > max){
380 max = f(i);
381 }
382
383 //минимум
384 if (f(i) < min){
385 min = f(i);
386 }
387 }
388
389 //масштаб
390 scaley = 100/(Math.abs(max) + Math.abs(min));
391
392 if (Math.abs(b) >= Math.abs(a)){
393 scalex = 250/b * 0.9;
394 }
395
396 if (Math.abs(b) < Math.abs(a)){
397 scalex = 250/Math.abs(a) * 0.9;
398 }
399
400 //сдвиг
401 dx = (b + a)/2;
402 dy = (max + min)/2;
403 dt = (b - a)/number;
404
405 }
406
407 function calculation(){
408
409 if (menu.value == 'first'){
410 console.log('a');
411 for (var i = a; i < b; i += dt){
412
413 if (f(i) >= 0 & f(i + dt) >= 0){
414 result += Math.abs(f(i))*dt;
415 }
416
417 if (f(i) < 0 & f(i + dt) < 0){
418 result -= Math.abs(f(i))*dt;
419 }
420
421 else{
422 result += 0;
423 }
424
425 }
426 }
427
428 if (menu.value == 'second'){
429 for (var i = a + dt; i < b; i += dt){
430 if (f(i) >= 0 & f(i + dt) >= 0){
431 result += Math.abs(f(i))*dt;
432 }
433
434 if (f(i) < 0 & f(i + dt) < 0){
435 result -= Math.abs(f(i))*dt;
436 }
437
438 else{
439 result += 0;
440 }
441 }
442 }
443
444 if (menu.value == 'third'){
445 for (var i = a + dt/2 ; i < b; i += dt){
446 if (f(i) >= 0 & f(i + dt) >= 0){
447 result += Math.abs(f(i))*dt;
448 }
449
450 if (f(i) < 0 & f(i + dt) < 0){
451 result -= Math.abs(f(i))*dt;
452 }
453
454 else{
455 result += 0;
456 }
457 }
458 }
459
460 }
461
462 function draw(){
463
464 //оси
465 ctx.beginPath();
466 ctx.moveTo(0 - dx * scalex, -h/2);
467 ctx.lineTo(0 - dx * scalex, h/2);
468 ctx.moveTo(w/2 , 0 + dy * scaley);
469 ctx.lineTo(-w/2 , 0 + dy * scaley);
470 ctx.stroke();
471
472 ctx.font = '10px Arial';
473 ctx.strokeText('1', (1 - dx) * scalex, (-0.1 + dy) * scaley);
474 ctx.strokeText('1', (0.1 -dx) * scalex, (-1 + dy) * scaley);
475
476 //кривая
477 ctx.beginPath();
478 for (var i = a; i < b; i+=0.01){
479 ctx.moveTo((i - dx) * scalex, (-f(i) + dy) * scaley);
480 ctx.lineTo((i + 0.01 - dx) * scalex, (-f(i + 0.01) + dy) * scaley);
481 }
482 ctx.stroke();
483
484 ctx.strokeStyle = 'red';
485
486 if (menu.value == 'first'){
487 for (var i = a; i < b - dt/2; i += dt){
488 if(f(i) >= 0){
489 ctx.strokeRect((i - dx) * scalex, (-f(i) + dy) * scaley, dt * scalex, Math.abs(f(i)) * scaley);
490 }
491 else{
492 ctx.strokeRect((i - dx) * scalex, (dy) * scaley, dt * scalex, Math.abs(f(i)) * scaley);
493 }
494 }
495 }
496
497 if (menu.value == 'second'){
498 for (var i = a + dt; i < b + dt/2; i += dt){
499 if(f(i) >= 0){
500 ctx.strokeRect((i - dx) * scalex, (-f(i) + dy) * scaley, dt * scalex, Math.abs(f(i)) * scaley);
501 }
502 else{
503 ctx.strokeRect((i - dx) * scalex, (dy) * scaley, dt * scalex, Math.abs(f(i)) * scaley);
504 }
505 }
506 }
507
508 if (menu.value == 'third'){
509 for (var i = a + dt/2; i < b + dt/4; i += dt){
510 if(f(i) >= 0){
511 ctx.strokeRect((i - dt/2 - dx) * scalex, (-f(i) + dy) * scaley, dt * scalex, Math.abs(f(i)) * scaley);
512 }
513 else{
514 ctx.strokeRect((i - dt/2 - dx) * scalex, (dy) * scaley, dt * scalex, Math.abs(f(i)) * scaley);
515 }
516 }
517 }
518
519 }
520
521 knopka.onclick = function(){
522
523
524 if (flag){
525 location.reload();
526 }
527
528 flag = true;
529 knopka.value = 'Очистить';
530 input();
531 calculation();
532 draw();
533 result1.innerHTML = Math.round(result * 100)/100;
534
535 }
536 }
