بایگانی‌های دانلود کد آماده - ايران متلب

اموزش متلب_pdepe(حل عددی معادله حرارت پاره ای وابسته به زمان)

ایران متلب — Wed, 11 Mar 2015 08:20:55 +0000

امروز قراره معادلات گرمایی در فضای تک بعد رو با چند مثال بررسی کنیم.

برای حل معادله گرما در زمان t و در فضای تک بعدی از دستور pdepe استفاده می کنیم.

فرم کلی معادله سھمی وار متلب بصورت زیر می باشد .

و ھمین طور فرم کلی شرایط مرزی بصورت زیر می باشد.

که x_l نمایانگر نقطه ابتدا و x_r نمایانگر نقطه انتهایی شرایط مرزی می باشد.توجه کنید که b دارای مقدار ثابت در هر دو معادله است.هر کدام از معادله های بالا را در m.file جداگانه قرار می دهیم.مقادیر s,b,c را از مقایسه معادله حاکم در مثال با فرم کلی معادله سهمیگون در متلب بدست اورده و در m.file قرار می دهیم.مقادیر p,q را نیز از مقایسه شرایط مرزی مسئله با فرم کلی شرایط مرزی تعریف شده برای متلب بدست می اوریم.و در اخر نیز شرایط کرانه ای را در m.file سوم قرار می دهیم.سپس به کمک دستور pdepe سه m.file را ترکیب کرده معادله را حل می کنیم.

:Example

معادله اصلی بصورت:

شرایط مرزی بصورت:

و شرایط اولیه بصورت:

کاملا مشخص است که مقدار c=1 می باشد.همچنین ما در طرف راست معادله حاکم بر مسئله مشتق دوم بر حسب x را داریم.پس مقدار b باید برابر u_x باشد.مقدار s نیز که به وضوح مشخص است که برابر 0 می باشد. ملاحضه کردید که مقادیر c,b,s طوری در معادله سهمیگون مقدار دهی شدند که به معادله حاکم بر مسئله برسیم.مقدار m نیز باید برابر 0 باشد تا عملا x از معادله حذف شود. m را در اخر مقدار دهی میکنیم.

کد زیر نمایانگر معادله اصلی و در eqn1.m ذخیره می کنیم.

function [c,b,s] = eqn1(x,t,u,DuDx)

%EQN1: MATLAB function M-file that specifies

%a PDE in time and one space dimension.

c = 1;

b = DuDx;

s = 0;

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

به شرایط ابتدایی دقت کنید.مقدار q برابر 0 است تا b حذف شود.در متیجه p برابر با u می شود.

و برای شرایط انتهایی نیز به همین صورت, q برابر صفر و مقدار p برابر با u-1 است.حالا نتایج رو در فایل bc1 ذخیره می کنیم.

function [pl,ql,pr,qr] = bc1(xl,ul,xr,ur,t)

%BC1: MATLAB function M-file that specifies boundary conditions

%for a PDE in time and one space dimension.

pl = ul;

ql = 0;

pr = ur-1;

qr = 0;

و همینطور شرط اولیه را در initial1.m ذخیره می کنیم.

function value = initial1(x)

%INITIAL1: MATLAB function M-
le that speci
es the initial condition

%for a PDE in time and one space dimension.

value = 2*x/(1+x^2);

و حالا نوبت به حل مساله به کمک دستور pdepe رسیده است.مساله را در زمان 10ثانیه نخست و طول واحد(مثلا یک لوله به طول واحد) حل می کنیم.

%PDE1: MATLAB script M-file that solves and plots

%solutions to the PDE stored in eqn1.m

m = 0;

% NOTE: m=0 specifies no symmetry in the problem. Taking

% m=1 specifies cylindrical symmetry, while m=2 specifies

% spherical symmetry.

% Define the solution mesh

x = linspace(0,1,20);

t = linspace(0,10,10);

%Solve the PDE

u = pdepe(m,@eqn1,@initial1,@bc1,x,t);

% Plot solution

figure(1)

surf(x,t,u);

title(‘Surface plot of solution’);

xlabel(‘Distance x’);

ylabel(‘Time t’);

figure(2)

for i=1:length(t)

plot(x,u(i,:))

hold on

end

جواب u بصورت یک ماتریس به اندازه t و x ذخیره شده است.برای مثال (u(1,5 مقدار دما در نقطه((t(1),x(5)

مشخص می کند.

ما می تونیم به کمک دستور((:,plot(x,u(1 بردار دما را در اغاز (t=0) رسم می کنیم.

مثال بعدی :

وکد زیر رو برای مساله به این صورت می نویسیم:

function [c,f,s] = eqn2(x,t,u,DuDx)

c = pi^2;

f = DuDx;

s = 0;

function [pl,ql,pr,qr] = bc2(xl,ul,xr,ur,t)

pl = ul;

ql = 0;

pr = pi * exp(-t);

qr = 1;

function u0 = initial2(x)

u0 = sin(pi*x);

و از دستور pdepe استفاده می کنیم:

m = 0;

x = linspace(0,1,20);

t = linspace(0,5,20);

u = pdepe(m,@eqn3,@initial3,@bc3,x,t);

figure(1)

surf(x,t,u);

title(‘Surface plot of solution’);

xlabel(‘Distance x’);

ylabel(‘Time t’);

figure(2)

plot(x,u(2,:))

title(‘Solution at t = 2’)

xlabel(‘Distance x’)

ylabel(‘u(x,2)’)

% Creat animation

figure(3)

fig = plot(x,u(1,:),’erase’,’xor’);

for k=1:length(t)

set(fig,’xdata’,x,’ydata’,u(k,:),’color’,’r’)

pause(0.1)

end

دریافت کتابچه حل معادلات ديفرانسيل با مشتقات پاره اي با MATLAB

نوشته اموزش متلب_pdepe(حل عددی معادله حرارت پاره ای وابسته به زمان) اولین بار در ايران متلب. پدیدار شد.

حل معادلات PDE با pdetool متلب

ایران متلب — Mon, 09 Mar 2015 20:03:25 +0000

در این قسمت می خواهیم شما را با جعبه ابزار pdetool متلب که برای حل معادلات دبفرانسیل جزیی (Partial Differential Equation) آشنا کنیم.

همان طور که می دونید برای حل معادلات PDE ما نیاز داریم که ناحیه ای را مشخص کنیم تا معادله ما در این ناحیه حل شود. به عنوان مثال می خواهیم تغییرات گرما در یک ناحیه مستطیلی را بدست بیاوریم. می بایست ابتدا یک مستطیل تعریف کنیم یا بکشیم و سپس آنرا مش بندی کنیم.

در مرحله اول دستور pdetool را در پنجره فرمان بزنید :

با این کار پنجره محیط pdetool باز می شود و شما می توانید ناحیه های دلخواه خود را رسم کنید.

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

pderect([-0.5 1 -0.5 1],’R1′)

یک مستطیل می کشد که گوشه چپ پایین آن در مختصات -0.5 و -0.5 قرار دارد و گوشه بالای سمت راشت در 1 و 1 قرار دارد. بر چسب آن را R1 می گذارد.

pderect([-0.2 0.8 -0.2 0.8],’R2′)

یک مستطیل می کشد که گوشه چپ پایین آن در مختصات -0.2 و -0.2 قرار دارد و گوشه بالای سمت راست در 0.8 و 0.8 قرار دارد. بر چسب آن را R2 می گذارد.

pdecirc(0.3,0.3,0.4,’C1′)

یک دایره به مرکز 0.3 و 0.3 می کشد که اندازه شعاع آن 0.4 است.

بعد از مشخص کردن ناحیه ها، می توانیم یک ناحیه را از دیگری کم کنیم و بین آنها را بگیریم.

برای انجام این رابطه ها از نوار ابزار set formula استفاده کنید و فرمول مورد نظرتون را وارد کنید.

مثلا در اینجا ما وارد می کنیم :

R1-R2+C1-C2

در مرحله بعد باید مش بندی کنیم برای این کار از منو Mesh

Initialize mesh را انتخاب می کنیم :

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

شما می توانید این مش بندی را تغییر دهید. مش بندی مثلثی پیش فرض نوع مش بندی در pdetool در متلب می باشد.

همانطور که در شکل بالا می بینید مثلث های مش تقریبا بزرگ هستند. برای کوچک تر کردن مش و دقت حل معادله، می توانید از منو mesh

refine mesh را انتخاب کنید. که مثلث ها را بر چهار تقسیم می کند.

پست قبلی ما در مورد حل معادلات مشتق جزیی در متلب را نیز ببینید.

https://matlab1.ir/?p=1369

نوشته حل معادلات PDE با pdetool متلب اولین بار در ايران متلب. پدیدار شد.

دانلود کتاب آشنایی با نرم افزار های مفید ریاضی

ایران متلب — Thu, 31 Jul 2014 12:33:41 +0000

همانطور که می دانید نرم افزارهای زیادی جهت انجام محاسبات ریاضی پیچیده و علمی وجود دارد که هر یک مزایا و معایب خاص خود را دارند .
بدون شک ۳ نرم افزار پرکاربرد در این زمینه MATLAB ، Mathematica و Maple هستند. اما شاید بیشتر افراد با تفاوت ها و مزایای هر یک از این نرم افزارها نسبت به سایرین آشنا نباشند.
کتابی که در این پست برای کاربران قرار داده ایم، توسط دکتر عقیله حیدری، از اساتید دانشگاه پیام نور، نوشته شده است. در این کتاب به تفصیل به بررسی ویژگی ها و مزایا و معایب هر یک از این نرم افزارها پرداخته شده است.

تعداد صفحات کتاب : ۱۸۰ صفحه

اندازه فایل : 3MB

دانلود فايل

نوشته دانلود کتاب آشنایی با نرم افزار های مفید ریاضی اولین بار در ايران متلب. پدیدار شد.

Particle filter (فیلتر ذره)

ایران متلب — Thu, 31 Jul 2014 10:48:54 +0000

یکی از ابزارهای خیلی پر کاربرد در پردازش ویدئو particle filter می باشد . که خیلی زیاد در کاربرد دنبال کردن (tracking) در ویدئو مورد استفاده قرار می گیرد.

ما در این پست می خواهیم برای شما کدهای تست شده MATLAB فیلتر ذره (particle filter) را قرار دهیم و یک توضیح خلاصه از روند الگوریتم آن برای شما بیاوریم.

دو PDF بسیار بسیار عالی که با زبانی ساده و روان مفهوم particle filter را آموزش می دهد :

PDF 1

Algorithm development using MATLAB – Object tracking by particle filter –

PDF 2

Object Tracking

رديابي اشياء(Object Tracking)

سيستمهاي هوشمند، سيستمهايي هستند كه خود به تنهايي و بطور خودكار اعمال خاصي انجام مي‌دهند . در برخي از اين سيستمها ، رديابي از بخشهاي عمده سيستم محسوب مي‌شود . سيستمهاي نظامي هوشمند و روبوت‌ها نمونه‌هايي از اين سيستمها هستند.

رديابي روشهاي خاصي دارد كه براي موارد گوناگون قابل استفاده است. روشهاي تخمين Kalman filter ، استفاده از Particle filter ، روشهاي Wavelet و Segmentation از جمله روشهاي قابل استفاده براي رديابي هستند .

موارد استفاده از آنها در سيستمهاي تعيين موقعيت جهاني(GPS) ، رادار ، رديابي هدف ، ماشين بينايي و كاربردهاي مربوط به رباتيك است.

چند کد نوشته شده در m-file متلب برای فیلتر ذره ای

دانلود پنچ كد آماده فيلتر ذره

نوشته Particle filter (فیلتر ذره) اولین بار در ايران متلب. پدیدار شد.