Редактирование: Теплопроводность с краевым условием

==Авторы==
[[Федоров_Евгений|Федоров Евгений]]

==Цель работы==
Визуализация процесса теплопередачи в неравномерно нагретой пластине. Постановка задачи в приложенном файле.

https://drive.google.com/file/d/1bhrHtb-4N6wRwiDy_Z7zCG0OoFWn8fS9/view?usp=sharing
==Визуализация==
*Визуализация решения первого уравнения
{{#widget: Iframe | url=http://tm.spbstu.ru/htmlets/js2020/Fedorov/newpj4.html | width=1000 | height=1000 | border = 0 }}
*Визуализация решения второго уравнения
{{#widget: Iframe | url=http://tm.spbstu.ru/htmlets/js2020/Fedorov/newpj6.html | width=1000 | height=1000 | border = 0 }}

== Код программы первого решения ==
<div class="mw-collapsible mw-collapsed">
'''JS:''' <div class="mw-collapsible-content">
<syntaxhighlight lang="javascript" line start="1" enclose="div">
function main() {
    var ctx = canvas_example.getContext('2d');
	var s = 400;
	var p = 400;
	var a = 1;
	var b = 4;
	function fe(x,y) {
		return Math.sin(b*x*Math.PI/p)*Math.sin(b*y*Math.PI/s);
	}
	var i = [];
	var m = [];
	for (var k = 1; k<1+b; k++) {
		i[k] = [];
		m[k] = [];
		for (var n = 1; n<1+b; n++) {
			var t = 0;
			var L = 0;
			for (var x = 1; x<1+p; x++) {
				for (var y=1; y<1+s; y++) {
					t = t + fe(x,y)*Math.sin(x*k*Math.PI/p)*Math.sin(y*n*Math.PI/s);
				}			
			}
			L = Math.PI*Math.PI*k*k/p/p+Math.PI*Math.PI*n*n/s/s;
			i[k][n] = t*4/p/s;
			if (Math.abs(i[k][n]) < 0.001) {
				i[k][n] = 0;
			}
			var j = 1/Math.exp(a*a*L);
			m[k][n] = j;
		}
	}
	var dc = 255/2;
	var r = 1;
	var xn = 100;
	var yn = 100;
    var vr = 0;
	function control() {
			
	for (x=1; x<1+p; x++) {
		for (y=1; y<1+s; y++) {
	        var ot=0;
			var q=0;
			for (k = 1; k<1+b; k++) {
				for (n = 1; n<1+b; n++) {
			     ot = ot + i[k][n]*Math.sin(x*k*Math.PI/p)*Math.sin(y*n*Math.PI/s)*Math.pow(m[k][n],vr);
				}
			}
			ot = (ot+1)*dc;
			c = `rgb(${ot},${ot/2},${0})`;
        ctx.fillStyle = c;
        ctx.beginPath();
        ctx.arc(x+xn, y+yn, r, 0, 0.5 * Math.PI);
        ctx.fill();
		}
	}
	vr = vr + 100;
	}
	setInterval(control, 300)
}
</syntaxhighlight>
</div>

== Код программы второго решения==
<div class="mw-collapsible mw-collapsed">
'''JS:''' <div class="mw-collapsible-content">
<syntaxhighlight lang="javascript" line start="1" enclose="div">
 
window.addEventListener('load', main, false);

function main() {
    var ctx = canvas_example.getContext('2d');
	var s = 400;
	var p = 400;
	var a = 1;
	var b = 3;
	function fe(x,y) {
		return x*(x-s)*y*(y-p);
	}
	var i = [];
	var m = [];
	for (var k = 1; k<1+b; k++) {
		i[k] = [];
		m[k] = [];
		for (var n = 1; n<1+b; n++) {
			var t = 0;
			var L = 0;
			for (var x = 1; x<1+p; x++) {
				for (var y=1; y<1+s; y++) {
					t = t + fe(x,y)*Math.sin(x*k*Math.PI/p)*Math.sin(y*n*Math.PI/s);
				}			
			}
			L = Math.PI*Math.PI*k*k/p/p+Math.PI*Math.PI*n*n/s/s;
			i[k][n] = t*4/p/s;
			if (Math.abs(i[k][n]) < 0.001) {
				i[k][n] = 0;
			}
			var j = 1/Math.exp(a*a*L);
			m[k][n] = j;
		}
	}
	var max = s*s*p*p/16;
	var dc = 255/max;
	var r = 1;
	var xn = 100;
	var yn = 100;
    var vr = 0;
	function control() {
			
	for (x=1; x<1+p; x++) {
		for (y=1; y<1+s; y++) {
	        var ot=0;
			var q=0;
			for (k = 1; k<1+b; k++) {
				for (n = 1; n<1+b; n++) {
			     ot = ot + i[k][n]*Math.sin(x*k*Math.PI/p)*Math.sin(y*n*Math.PI/s)*Math.pow(m[k][n],vr);
				}
			}
			ot = (ot+1)*dc;
			c = `rgb(${ot},${ot/2},${0})`;
        ctx.fillStyle = c;
        ctx.beginPath();
        ctx.arc(x+xn, y+yn, r, 0, 0.5 * Math.PI);
        ctx.fill();
		}
	}
	vr = vr + 1000;
	}
	setInterval(control, 300)
}
 
</syntaxhighlight>
</div>