روش کار متلب برای حل معادلات پاره ای رامی توانید در آدرس زیر ببنید

http://www.mathworks.com/access/helpdesk/help/toolbox/pde/ug/bqivs1t-1.html

اگر ارل بالا را باز کنید می‌بینید که برای حل pde ابتدا باید ماتریس هندسه مسئله (Geometry Description matrix)  را بسازیم. بعد  با استفاده از دستور decsg این ماتریس را تجزیه می‌کنیم که به ماتریس حاصل ماتریس هندسه تجزیه شده می‌گویند(Decomposed Geometry matrix). این ماتریس ناحیهای (domain) که معادله بر روی آن حل می شود را مشخص میکند

مثال زیر روش بدست آوردن این دو ماتریس را نشان می‌دهد.

مثال)

s = tf('s');

sys = 1/(s+1)^2;

nyquist(sys)

می بینید که نایکویست این تابع شبیه قلب است؛ و همانطور که می‌دانید نمودار نایکویست در صفحه اعداد مختلط رسم می‌شود. برای حل  یک معادله پاره ای را بر روی این شکل باید آن را از فضای اعداد مختلط به فضای اعداد حقیقی تبدیل کنیم(جزییات کار به عهده خودتان).

w=-10:.1:10;

x=(1-w.^2)./(1+w.^2).^2;

y=2*w./(1+w.^2).^2;

plot(x,y)

بد نشد ولی یک مقدار تمیز کاری می‌خواهد.

w = [linspace(-5,-2,7) linspace(-1.9,-.2,20) -.15:.05:.15...

    fliplr(-linspace(-1.9,-.2,20)) fliplr(-linspace(-5,-2,7))];

x = (1-w.^2)./(1+w.^2).^2;

y = 2*w./(1+w.^2).^2;

حالا یک مقدار چرخش

xy = [0 x ;0 y ];

XY = [0 1;-1 0]*xy;

اگر به معادلات x,y دقت کنید می بینید هنگامی که w به سمت بینهایت میرود x,y برابر با صفر می شوند ولی در عمل نمی توان w را به سمت بی نهایت میل داد٬ پس صفر را خودمان اضافه کردیم.

هنوز یک مقدر پهن است.

XY = [0.7 0;0 1]*XY;

الان باید ماتریس هندسه را بنویسم

هر ستون این ماتریس مربوط به یک شکل است یعنی اگر بخواهیم معادله را همزمان برای دو شکل حل کنیم ماتریس دو ستون خواهد داشت.

در این مثال یک شکل داریم و ماتریس یک ستونی است.

اگر هندسه :

به شکل دایره باشد در سطر اول عدد ۱ ٬ در سطر دوم و سوم به ترتیب مختصات x,y مرکز دایره و در سطر ۴ شعاع دایره را می نویسیم.

به شکل چند ضلعی باشد در سطر اول عدد ۲ ٬ در سطر دوم تعداد اضلاع(n) ٬ در n سطر بعدی xها و در n سطر بعد از آن y ها می نویسیم.

به شکل چهار ضلعی باشد در سطر اول  عدد ۳ ٬ را می نویسیم. سطرهای بعدی را مانند حلت چند ظلعی پر می کنیم.

به شکل بیضی باشد در سطر اول عدد ۴ ٬ در سطر دوم و سوم به ترتیب مختصات x,y مرکز بیضی در سطر چهار و پنج اندازه قطر بزرگ و کوچک بیضی و در سطر ششم زاویه چرخش قرار می گیرد.

حالا می‌توانیم ماتریس هندسه  بنویسیم.

شکل ما از نوع چند ضلعی است پس در سطر اول عدد 2 را می‌نویسیم.توابع x,y دو تا صفر دارند، یکی در بی‌نهایت و یکی در منفی بی‌نهایت ولی ما فقط یکی را اضافه کردیم چون هر نقطه را فقط 1 بار باید بیاوریم.

x = XY(1,:)';

y = XY(2,:)';

n = length(x);

gd = [2;n;x;y];

dl = decsg(gd);

تا اینجا توانستیم که ماتریس هندسه و ماتریس هندسه تجزیه شده را بدست آوریم. اگر به ارلی در اول مقاله آوردم نگاه کنید می‌بیند که 3 مرحله اول حل الگوریتم معادلات پاره‌ای را انجام داده‌ایم. این 3 مرحله را می‌توان در یک ام-فایل خلاصه کرد. برای نوشتن این ام-فایل از دستور wgeom استفاده می‌کنیم.

fid=wgeom(dl,'cardiodgeom')

این دستور یک ام-فایل با اسم cardiodgeom ایجاد می‌کند. اگر به هر دلیلی این فایل ایجاد نشود مقدار fid برابر با -1  می‌شود.

ادامه دارد...

 

 

 

[ دوشنبه 1 مرداد 1386 - 07:07 ق.ظ ]
[ویرایش شده در : یکشنبه 23 خرداد 1389 - 07:04 ق.ظ]

[ پیام ()|| امین باشی ] [حل pde , ] [+]

مثال مطلب قبلی در نظر بگیرید٬ یک بار دیگر می خواهیم این مثال را حل کنیم.

اولین تغییری که در حل مثال می دهیم٬ نحوه شبکه بندی مسئله است.به نظر می آید که تعداد گره ها برای حل مسئله کافی نباشد با استفاده از دستور refinemesh تعداد گره ها بیشتر می کنیم.با زیاد  شدن گره ها تعتداد محاسبات و زمان آن بالا تر می رود ولی ممکن است که نتیجه کار چندان تغییر نکند و یا بهتر باشد که نوع شبکه را عوض کنیم. برای درک بهتر می توانید به کتاب هایی در مورد روش تفاضل محدود-متلب از این روش استفاده می کند- مراجعه کنید.

  

[p,e,t]=initmesh(g);

[p,e,t]=refinemesh(g,p,e,t);

x=p(1,';

y=p(2,';

u0=atan(cos(pi/2*x));

ut0=3*sin(pi*x).*exp(cos(pi*y));

tlist=linspace(0,5,50);

uu=hyperbolic(u0,ut0,tlist,b,p,e,t,1,0,0,1);

[ دوشنبه 9 بهمن 1385 - 04:01 ق.ظ ]
[ویرایش شده در : - - -]

[ پیام ()|| امین باشی ] [حل pde , ] [+]

[ پنجشنبه 5 بهمن 1385 - 09:01 ق.ظ ]
[ویرایش شده در : دوشنبه 9 بهمن 1385 - 05:01 ق.ظ]

[ پیام ()|| امین باشی ] [حل pde , ] [+]

در مطالب قبلی مقدمات استفاده از pdetoolbox را مرور کردیم و فرض می کنیم که بلد هستید جواب pde را در متغیری مثل u و تقسیم بندی مسئله  (mesh) را در متغیر هاب p,e,t ذخیره کنید.

دستور pdemesh برای رسم تقسیم بندی استفاده می شود.

pdemesh(p,e,t)
pdemesh(p,e,t,u)
h=pdemesh(p,e,t)
h=pdemesh(p,e,t,u)

حالات اول دستور٬ تقسیم بندی سیستم را به صورت دو بعدی رسم می کند.

حالت دوم این کار رابه صورت سه بعدی نمایش می دهد٬ ارتفاع هر گره به متغیر u (دما) بستگی دارد.

 اما اگر بخواهیم دمای هر گره را بر روی خود آن نشان دهیم به صورت عمل می کنیم:

pdemesh(p,e,t)

hold on

متغیر p دو سطر دارد که سطر اول مولفه x و سطر دوم مولفه y است.تعداد ستون های p به اندازه تعداد گره ها است.

متغیر u یک ستون و به تعداد گره ها سطر دارد. یعنی هر سطر آن متناظر با یک ستون p است.

حالا متغیر u  را بتبدیل به رشته گرده و بر روی شکل قرار می دهیم

uu=num2str(u);

nn=length(u);

for i=1:nn

text(p(1,i),p(2,i),uu(i,:))

drawnow

end

---

دوستی پرسیدن که در هنگام استفاده از plot  با پیغام خطای زیر روبرو می شوند:

Attempt to execute SCRIPT plot as a function 

احتمالا این دوستمون از یک متغیر یا یک ام فایل به اسم plot استفاده کرده است. برای اینکه با این خطا و خطاهایی شبیه به این روبرو نشوید اسم متغیرها را به دقت انتخاب کنید  

[ پنجشنبه 14 دی 1385 - 09:01 ق.ظ ]
[ویرایش شده در : - - -]

[ پیام ()|| امین باشی ] [حل pde , ] [+]

حل معادلات دیفراسیل جزیی قسمت مهمی از ریاضیات مهندسی است در متلب جعبه ابزار جداگانه ای برای حل اینگونه  معادلات قرار داده شده است .مرور کردن این جعبه ابزار به اندازه نوشتن یک کتاب وقت میگیرد.در فایل زیر توضیحات مختصری در مورد نحوه حل این گونه معادلات را نوشته ام.

http://matlabedu.googlepages.com/heat.pdf

[ چهارشنبه 29 آذر 1385 - 06:12 ق.ظ ]
[ویرایش شده در : - - -]

[ پیام ()|| امین باشی ] [حل pde , ] [+]

+ محاسبه سری فوریه به صورت عددی
+ تعریف هندسه pde
+ معادلات هذلولی ۲
+ معادلات هذلولی
+ رسم pde
+ تبدیل معادلات سمبولیک به عددی
+ حل pde
+ حل تحلیلی معادلات جبری
+ حل معادلات خطی
+ BVP
+ چاپ محتویات یک دایرکتوری
+ pdepe
+ معادلات با درجه بالاتر
+ حل عددی IVP
+ fourier

---

صفحات :

---

هرگونه استفاده تجاری از مــطالب این سایت بصورت كتاب٬ نشریه٬ وب و ... ممنوع می‌باشد

All right reserved©2005 Amin Bashi