Green's Function Solution in Matlab
162 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
I wrote a code for stationary temperature distribution which is;
f = 0.09;
b = 0.0044;
q = 3.73e-9;
L = 1;
Tw = 250;
Tam = 27;
syms c x g
T = 2*c*cosh(x*((f-b*g)/q)^0.5)+g/(f-b*q);
c = solve(subs(T,'x',L/2)==0,c);
z = simplify(int(subs(T,'c',c),x,-L/2,L/2));
g = solve(z==L*(Tw-Tam),g)
T
Now I try to use perturbation method and find T1 temperature distribution. My equation is;
q*T1''(x)-T1(x)*(f-b*g+i*w*p)=T*(f1-b*g1)-g1
And my prof. say that can be solved by Green's function G(x,y), where the G(x,y) is soluton of;
q*G''(x,y)-G(x,y)*(f-b*g+i*w*p)=Diract(x-y)
And
G(L/2)=0 , G(-L/2)=0
if I can find the G(x,y), I will get the solution of T1.
Can anyone help me please?
Thank you for everything
Yusuf
0 Kommentare
Antworten (2)
Sousheel reddy Lachagari
am 4 Apr. 2022
clc
clear all
syms x y
N = input('the value of N = ');
M = input('M=')
My=diff(M,y)
Nx=diff(N,x)
x1=input('limit of x');
x2=input(' x u');
y1=input('limit of y');
y2=input('limit of y up');
I=int(int(Nx-My,y,y1,y2),x,x1,x2);
1 Kommentar
Siehe auch
Kategorien
Mehr zu Symbolic Math Toolbox finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!