Green's Function Solution in Matlab

50 次查看（过去 30 天）
yusuf
yusuf 2014-2-4
评论： Walter Roberson 2022-4-4
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

请先登录，再回答此问题。

回答（2 个）

Ted Shultz
Ted Shultz 2019-10-16
The person who asked this question has some (slow) sample code that will solve Green's function.
Sousheel reddy Lachagari
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);

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by