Green's Function Solution in Matlab
조회 수: 164 (최근 30일)
이전 댓글 표시
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
답변 (2개)
Sousheel reddy Lachagari
2022년 4월 4일
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
참고 항목
카테고리
Help Center 및 File Exchange에서 Symbolic Math Toolbox에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!