Green's Function Solution in Matlab
18 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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 comentarios
Respuestas (2)
Sousheel reddy Lachagari
el 4 de Abr. de 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 comentario
Walter Roberson
el 4 de Abr. de 2022
Please add comments explaining how this implements Green's function.
Ver también
Categorías
Más información sobre Symbolic Math Toolbox en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!