در اين تاپيك سورس ها و آموزش هاي مطلب قرار داده ميشه.
سورس معمولا" مربوط به M-File ها هستش و سعي شده كه اين كد ها در اكثر نسخه هاي Matlab قابل اجرا باشن.
در ادامه آموزش هاي كامل و كاربردي هم از Matlab قرار داده ميشه.
موفق باشيد:give_rose:
ravegoat
01-06-10, 10:34
با قرار دادن دستور زير دريك mfile و اجراي آن يه شش ضلعي در پنجره ي Figure ترسيم مي شود.
اين سورس مربوط به MFile ، مقدار h را از كاربر (توسط پنجره Commnad) گرفته و به تعداد h پله ي عرش را كه تعبيري شهودي از دنباله ي واگراي n^-1 است را ترسيم مي كند:
clc
h=input('Please insert a positive integer: ');
clc
if h>0
clf
hold on
x=[0 0 1];
y=[0 1 1];
plot(x,y)
for i=2:h
x=x+[(i-1)^-1 (i-1)^-1 i^-1];
y=y+[(i-1)^-1 i^-1 i^-1];
plot(x,y)
end
else
error(':-(')
end
ravegoat
01-06-10, 10:45
اين سورس مقدار ضلع مربع و يك مقدار حدي را از طريق پنجره ي Command از كاربر دريافت كرده و مربعات ارشميدس را تشكيل مي دهد تا مقدار ضلع كوچك ترين مربع به عدد حدي تعيين شده برسد:
clc
r=input('Plz determine R: ');
a=input('Now plz insert A: ');
b=a;
clc
clf
hold on
while abs(a)>=r
x=[-b a a];
y=[b b -a];
plot(x,y)
b=-a;
a=-.8*a;
end
ravegoat
01-06-10, 10:49
اين كد MFile حركت دو بعدي يه توپ را روي منحني |sint/t| شبيه سازي مي كند:
clear all
s=clock;
clc
x=linspace(0,50,500);
y=(abs(sin(x))./x);
for i=1:500
clf
hold on
plot(x(1:i),y(1:i),'-')
plot(x(i),y(i),'o r')
axis([-1,51,-0.1,1.1])
m(i)=getframe;
end
f=clock-s;
fprintf('Finished in %g second(s) \n',(f(4)*3600 + f(5)*60 + f(6)))
clear all
ravegoat
06-06-10, 14:03
با قراردادن سورس زير در يك MFile جديد و اجراي آن مي توانيد اين شبيه سازي را در پنجره ي Figure مشاهده كنيد:
clear all
clc
clf
dt=linspace(0,pi,120);
for i=1:120
surfl(peaks)
shading interp;
for j=1:round(abs( sin(dt(i))*20))
light('Position',[cos(dt(i))*50 0 sin(dt(i))*10],'Style','infinite');
end
if i<20
title('Hello World')
elseif i>20 && i<50
title('Morning')
elseif i>50 && i<70
title('Noon')
elseif i>70 && i<100
title('Afternoon')
else
title('Evening')
end
m(i)=getframe();
end
title('Thats a simple day!')
movie(m,1,12)
ravegoat
24-06-10, 16:04
سورس زير رو مي تونيد تو يه MFile بريزيد و اجرا كنيد. با اجراي كد برنامه مقدار n رو از طريق Command Windows از كاربر دريافت مي كنه و سپس n عدد از كاربر دريافت مي كند و سپس آن ها را به صورت صعودي به نزولي نمايش مي دهد. در اين قطعه كد از توابع Sorting مطلب استفاده نشده چون هدف نمايش الگوريتم و منطق اين مرتب سازي بوده است:great::
clear all
clc
n=input('How many numbers you want to insert: ');
for i=1:n
clc
m(i)=input('Please insert number: ');
end
clc
c=n;
for i=1:n
p=m(1);
inx=1;
for j=2:c
if p<m(j)
p=m(j);
inx=j;
end
end
m(inx)=m(c);
c=c-1;
fprintf('%g \n',p)
end
ravegoat
02-07-10, 18:12
براي تعريف يك تابع در MATLAB بايد يك MFile جديد ايجاد كنيد. سپس دستورات را اين گونه داخل تابع وارد كنيد:
function FuncName(a,b,c,...,n)
%دستورات
end
FuncName نام تابع است.
مقادير a تا n نيز نام آرگومان ها هستند. (توجه كنيم كه منظور از n اين است كه به مقدار دلخواه مي توان براي تابع پارامتر تعريف كرد.)
سپس بايد MFile را با نام تابع (در اين جا FuncName) در Work Path ذخيره كرد. در اين صورت همواره به تابع دسترسي داريم.
توجه كنيد كه متغير هاي تعريف شده در درون تابع پس اجراي دستورات درون function پاك مي شوند. پس اگر به مقادير اين متغير ها نياز داريد با نوع متغير را global انتخاب كنيد.
در آينده سعي مي شود كه بيش تر در اين مورد توضيح دهيم.
به عنوان يك مثال ساده، تابع اي با نام kore تعريف مي كنيم كه با گرفتن مقادير a , b و c به عنوان مركز كره و مقدار r براي شعاع، آن كره را ترسيم مي كند:
function kore(a,b,c,r)
[x,y,z]=sphere(30);
x=x*r+a;
y=y*r+b;
z=z*r+c;
surfl(x,y,z)
end
حال با دستور زير يك كره با شكل دلخواه رسم مي شود:
kore(1,2,3,0.5)
ravegoat
03-07-10, 23:18
با تو جه به تابع كره در پست 9 (Only the registered members can see the link)در اين قسمت مي تونيم حركت يه كره بر مسير دايره اي رو شبيه سازي كنيم:
clear all
k=linspace(0,2*pi,100);
for j=1:100
z(j)=5;
end
for t=1:100
clf
plot3(5*sin(k),5*cos(k),z,'-.')
hold on
kore(5*sin(t),5*cos(t),5,0.3)
axis([-6,6,-6,6,4,6])
m(t)=getframe();
end
ravegoat
03-07-10, 23:21
در اين مثال هم از تابع kore استفاده مي كنيم و حركت نوساني يك توپ رو كه در مسيري فنري شكل حركت مي كنه رو نمايش مي دهيم. كافيه كه سورس زير رو در يه mfile قرار بديم و با اجراي كد شبيه سازي رو تو پنجره ي figure مشاهده كنيم:
clear all
k=linspace(0,2*pi,100);
for t=1:100
clf
plot3(5*sin(k),5*cos(k),k,'-.')
hold on
kore(5*sin(k(t)),5*cos(k(t)),k(t)+(sin(t)/t)*3,0.5)
axis([-6,6,-6,6,-1,7])
m(t)=getframe();
end
ravegoat
06-07-10, 13:57
اين سورس Mfile بيش تر جنبه ي مكانيكي-فيزيكي داره و نوسان يك فنر رو نشون ميده. بعد با دستور subplot نمودار هاي مكان-زمان و سرعت-زمان رو ترسيم مي كنه.
clear all
clc
t=linspace(0,8*pi,400);
x=sin(t);
y=cos(t);
figure(1)
for i=1:400
plot3(x,y,t+t*(sin(t(i))/t(i))*2*pi)
axis([-1 1 -1 1 -10 50])
m(i)=getframe;
end
figure(2)
subplot(1,2,1)
plot(t,t+t.*(sin(t)./t)*2*pi)
subplot(1,2,2)
plot(t,(cos(t)-sin(t).*t).*t.^2)
title('Position/Time - Velocity/Time')
yellowYeild
14-09-10, 23:14
با سلام خدمت دوستان ..
چنتا سوال دارم:
1- این سورسها در چه ورژنی آزمایش میشن؟
2-معنی این چیه؟
در مربع ارشمیدس یه جا نوشتید b=a و در انتها نوشتید b=-a !! این یعنی چی ؟ چجوری ممکنه؟
3-ظاهرا kore در متلب ورژن R8 تغییر کرده چون ارور میده ؟ میدونین باید چه چیزی به جاش نوشت؟
............................
خیلی ممنونم
من سوالات زیادی دارم واسه همین عضو شدم امیدوارم پاسخ بدید.
ravegoat
17-09-10, 11:34
با سلام خدمت دوستان ..
چنتا سوال دارم:
1- این سورسها در چه ورژنی آزمایش میشن؟
2-معنی این چیه؟
در مربع ارشمیدس یه جا نوشتید b=a و در انتها نوشتید b=-a !! این یعنی چی ؟ چجوری ممکنه؟
3-ظاهرا kore در متلب ورژن R8 تغییر کرده چون ارور میده ؟ میدونین باید چه چیزی به جاش نوشت؟
............................
خیلی ممنونم
من سوالات زیادی دارم واسه همین عضو شدم امیدوارم پاسخ بدید.
با سلام!
دوست عزيز عضويت شما در PCN رو تبريك ميگم.
1-نسخه ي 2008
2-خب الگوريتم ترسيم به اين شكل هست. براي بار اول دستور b=a اجرا ميشه تا دو خط از مربع رسم شه. بعد برنامه وارد يه حلقه ي while ميشه و در انتهاي هر بار اجراي حلقه دستور b=-a اجرا ميشه. در واقع فقط در يك ترسيم نياز به b=a و در بقيه موارد b=-a نياز است. اجراش كرده بودم درست كار مي كرد. با اين حال اگر تو هر جاييش مشكل داريد بفرماييد تا توضيح بدم خدمتتون.
3-دستوري به اسم kore در هيچ يك از نسخه هاي متلب نيست. بنده خودم يه بار چنين دستوري رو در يك mfile تعريف كردم و بعد اون kore.m رو داخل شاخه ي work قرار دادم. بدين شكل دستور kore براي متلب تعريف ميشه. چون رسم كره يه چيز كاربرديه تصميم گرفتم يه بار تابع ترسيمش رو تعريف كنم تا كار راحت تر بشه.
موفق باشيد
آرمين
yellowYeild
23-09-10, 02:39
با اين حال اگر تو هر جاييش مشكل داريد بفرماييد تا توضيح بدم خدمتتون.
ممنونم..
بنده خودم يه بار چنين دستوري رو در يك mfile تعريف كردم و بعد اون kore.m رو داخل شاخه ي work قرار دادم
میشه این دستور رو بذارید تا من هم استفاده کنم..
------------------------------------------------------------
سوالات دیگه:
1-چرا مینویسیم m(i)=getframe چرا با g یا u اجرا نمیشه؟ از پیش فرضهای نرم افزاره؟ ( به مانند اینکه حروف باید کوچک باشند -در برخی موارد-)
2- s=clock; در مورد این کمی توضیح میدید؟و f=clock-s;
.
.
من باز سوال دارم ..
بسیار ممنونم
ravegoat
23-09-10, 07:12
ممنونم..
میشه این دستور رو بذارید تا من هم استفاده کنم..
------------------------------------------------------------
سوالات دیگه:
1-چرا مینویسیم m(i)=getframe چرا با g یا u اجرا نمیشه؟ از پیش فرضهای نرم افزاره؟ ( به مانند اینکه حروف باید کوچک باشند -در برخی موارد-)
2- s=clock; در مورد این کمی توضیح میدید؟و f=clock-s;
.
.
من باز سوال دارم ..
بسیار ممنونم
خواهش مي كنم!
تابع مربوط به رسم كره دراينجا (Only the registered members can see the link) بود.
------------------------------------------------------------
1-اين سوال رو متوجه نميشم. ميشه دقيق تر توضيج بديد. تابع getframe درواقع عكسي رو از پنجره ي figure ميگيره و در آرايه ي m ذخيره مي كنه. متغير m يك آرايه است يعني متغير شامل تعداد متناهي عنصر. با دستور m(i) به عنصر i ام اين آرايه ميشه دسترسي پيدا كرد. حالا سوالتون دقيقا" در كدام مورد است؟
2-دستور clock در واقع زمان سيستم رو به هنگام اجراي دستور بر مي گردونه. كلا" اون كد براي تعيين مقدار زمان اجراي يك كد در mfile به كار ميره. ابتدا زمان شروع به كار برنامه در متغير s ذخيره ميشه سپس با پايان كد، مقدار s از زمان فعلي سيستم (clock) كه زمان پايان كار محسوب ميشه، كم ميشه و حاصل در متغير f قرار ميگيره. البته ميشه به جاي f=clock-s دستور etime(clock,s) رو هم به كار برد.
.
.
سوالات ديگه تون رو هم بفرماييد؛ تا اون جا كه مي تونم پاسخ ميدم:wink:.
آرمين
ravegoat
30-09-10, 08:28
پيچ ارشميدوس همون طور كه مي دونيد يه خم باز و ساده به شكل دايره اي هست كه به تدريج شعاع اون كم ميشه. صرف نظر از ظاهر جالب اين پيچ، كاربرد هاي مفيدي رو در بخش مكانيك داره. بنده سعي كردم يه كد ساده براش بنويسم كه البته كمي دقتش پايينه. ترسيم اين پيچ در مختصات قطبي خيلي ساده تر ميشه براي همين از دستور polar استفاده كردم:
همين طور كه ملاحظه مي كنيد اول شعاع رو تعريف كردم و مقدار اون رو از 10 تا صفر تغيير دادم و بعد مقدار زاويه تتا رو از 0 تا ده پي تغيير دادم تا تعداد دور هاي پيچ به 5 برسه كه البته مي تونيد بر حسب نياز خودتون اين مقادير رو تغيير بديد.
ravegoat
12-02-11, 20:46
نمونه ي سورس سري فوريه در نرم افزار MATLAB پيوست شده. اين سورس تابع f رو دريافت مي كنه در در بازه از a تا b به اندازه ي k تعداد جملات سري فوريه تابع را محاسبه مي كنه. البته رسم نمودارش يكم نياز به بهينه سازي داره...
ravegoat
12-02-11, 20:59
اين كد براي تحليل ارتعاش يك صفحه ي مرتعش هستش. PDE مربوط به چنين ارتعاشي نياز به حل انتگرال هاي دوگانه داره كه در اين جا از حل عددي انتگرال خود نرم افزار بهره گرفته شده. در ضمن پس از حل معادله ي ديفرانسيل پاره مي تونيد نحوه ي ارتعاش پوسته رو به صورت پويانمايي مشاهده كنيم. سورس اين امكان رو براي شما فراهم مي كنه تا شرايط اوليه رو بر حسب نياز خودتون تغيير بديد.
ravegoat
12-02-11, 21:04
اين سورس اختصاص به حل عددي PDE به روش Mesh بندي داره. يك معادله ي گرما به روش عددي (Numerical) و دقيق (Exact) حل ميشه تا شما بتونيد دو روش رو با هم مقايسه كنيد.
ravegoat
27-06-11, 00:53
اين سورس متلب بر اساس جنس انتخابي شما، در يك بازه ي معين قطر مفتول هاي مناسب براي طراحي فنر تحت بار استاتيكي فشاري و پارامتر هاي وابسته به آن (مانند طول بحراني، دور هاي موثر و ...) را نمايش مي دهد. در نهايت طراح با توجه به معيار هايي كه در نظر دارد مناسب ترين قطر را انتخاب مي كند. به علاوه اين سورس خروجي اكسل نيز دارد.
ravegoat
09-07-11, 22:37
سورسي كه در اين پست قرار داده شده مربوط به تحليل ارتعاشات يك سيستم دو درجه آزادي است. مسئله ي مربوطه از كتاب ارتعاشات ستو انتخاب شده است. اين مسئله با روش انرژي لاگرانژ حل و سپس خطي سازي شده است درنتيجه جواب مسئله براي سرعت هاي نسبتا" زياد با خطا همراه خواهد بود.
5593
ravegoat
12-07-11, 12:48
سورس هاي پيوست شده برنامه هايي هستند كه با توجه به اجزايي كه روي شفت قرار داده مي شوند، طراح را در جهت تعيين قطر مناسب شفت ياري مي كنند.
بخشي از برنامه نمودار گشتاور خمشي شفت را با توجه به مكان عناصري چون چرخ دنده، چرخ زنجير، بيرينگ و پولي ترسيم مي كند. كاربر مي تواند با تعيين مكان اين اجزا و توان انتقالي به عنوان ورودي، نمودار گشتاور خمشي را ملاحظه كند.
بخش ديگر برنامه يافتن قطر متوسط شفت است. اين كار بر اساس ضرايب مارين، رابطه ي Heywood براي تمركز تنش و همچنين تركيب معيار هاي von Mises و Soderberg صورت مي گيرد.
ravegoat
24-07-11, 15:10
بسته ي پيوست شده شامل 5 سورس متلب در زمينه ي مقاومت مصالح است. اين سورس در واقع حل مسايل كامپيوتري برگزيده ي كتاب مقاومت بير-جانستون است كه بخش هاي زير را در بر مي گيرد:
طراحي ستون بر اساس فرمول اويلر
ذخيره ي انرژي كرنشي در ميله
تغيير مكان تير ها
معيار شكست von Mises و Tresca
كرنش پيچشي در چرخدنده ها
ravegoat
06-08-11, 15:27
سورس پيوست شده مربوط به شبيه سازي يك سيستم آزادي يك درجه آزادي است. سيستم جرم و فنر به ميله ي دواري متصل شده كه با سرعت ثابت مي چرخد. چرخش ميله سبب ايجاد شتاب كريوليس در سيستم شده كه اين امر معادله ي قانون دوم نيوتون را ناهمگن مي كند همانند آن چه كه در ارتعاشات اجباري روي مي دهد. اين معادله ي ديفرانسيل با روش عددي ode45 در متلب حل شده و پاسخ به كمك توابع گرافيكي رسم مي شود.
ravegoat
04-11-11, 20:13
Simulink از جمله ابزار هاي قدرتمند نرم افزار MATLAB جهت شبيه سازي سيستم هاي مختلف نظير سيستم هاي ديناميكي، كنترلي و الكتريكي است.
در اين پست نمونه اي از شبيه سازي يك سيستم كنترلي Velocity Feedback با Simulink نشان داده شده است. در اين سيستم از خروجي يك فيدبك گرفته شده و اين مقدار در يك عدد به عنوان Gain ضرب مي شود. كاربر مي تواند با تغيير Gain ها، سعي كند سيستم مرتبه ي II را به طور بهينه كنترل كند.
نسخه ي 2008 Simulink اين سيستم پيوست شده است.
pakniat.s
11-11-11, 13:50
function kore(a,b,c,r)
با سلام ،
آرایه رو می تونیم به عنوان آرگومان ورودی به تابع بدیم( تو matlab )
ravegoat
11-11-11, 14:24
function kore(a,b,c,r)
با سلام ،
آرایه رو می تونیم به عنوان آرگومان ورودی به تابع بدیم( تو matlab )
با سلام!
عضويت شما رو در SSA تبريك مي گم.
بله چنين چيزي امكان پذير است.
function kore(arr)
a=arr(1);
b=arr(2);
c=arr(3);
r=arr(4);
end
موفق باشيد
آرمين
pakniat.s
11-11-11, 21:05
ممنون از راهنمایی تون ، همون طور که مطلع هستین تو آرایه ما با اندیس کار می کنیم تو مطلب هم همینطوره با خیر؟ و اینکه کتاب خوبی در مورد matlab سراغ دارید
ممنون
ravegoat
11-11-11, 22:32
ممنون از راهنمایی تون ، همون طور که مطلع هستین تو آرایه ما با اندیس کار می کنیم تو مطلب هم همینطوره با خیر؟ و اینکه کتاب خوبی در مورد matlab سراغ دارید
ممنون
خواهش مي كنم.
بله دقيقا" انديس براي آرايه در متلب هم وجود داره (در كدي هم كه به عنوان مثال قرار داده شده بود، از انديس هاي 1 تا 4 براي گرفتن درايه هاي يه آرايه استفاده شد!).
به علاوه براي آرايه هاي چند بعدي هم ميشه انديس هاي هر بعد رو با علامت ويرگول ( , ) از هم جدا كرد.
در زمينه ي كتاب بايد عرض كنم كه نرم افزار متلب بخش هاي متنوعي داره و منطقي هستش كه براي آشنايي بهتر با هر كدوم از اين بخش ها كتاب مربوط به اون بخش تهيه بشه.
با اين حال در زمينه ي كليات كار با MATLAB به نظر بنده كتاب MATLAB Demystified از David McMahon (انتشارات McGraw-Hill) مي تونه مفيد باشه. يك كتاب فارسي هم در بازار موجوده كه متاسفانه اسم نويسنده اش رو فراموش كردم (سعي مي كنم اسمش رو بپرسم).
يه PDF فارسي و خلاصه از متلب هم ضميمه شده كه اميدوارم براي شروع كار كمكتون كنه.
موفق باشيد
آرمين
ravegoat
13-12-11, 06:10
Root Locus روشی در مهندسی کنترل جهت تحلیل رفتار یک سیستم به ازای تغییر مقدار یک Gain به شمار می آید. در MATLAB کافیست تنها با تعریف ضرایب چند جمله های صورت و مخرج تابع تبدیل حلقه باز، مکان هندسی ریشه ها را رسم نمود:
a=[0 1 2];
b=[1 2 3];
rlocus(a,b)
در مثال فوق 1s + 2 در صورت کسر و 1s^2 + 2s + 3 در مخرج کسر قرار دارد.
ravegoat
22-01-12, 15:45
نمودار بوده (Bode Diagram) نموداري مفيد براي تحليل سيستم ها در نظريه ي كنترل است. بر اساس اين نمودار مي توان پاسخ يه سيستم نسبت به فركانس هاي مختلف را پيش بيني كرد و دريافت كه سيستم چه فركانس هايي را فيلتر مي كند.
براي رسم نمودار بود همانند مثال مكان هندسي، صورت و مخرج تابع تبديل را تعيين كرده و سپس با دستور bode نمودار دامنه و فاز رسم مي گردد:
a=[0 1 2];
b=[1 2 3];
bode(a,b)
ravegoat
22-01-12, 15:52
با رسم تابع تبديل حلقه باز يك سيستم در فضاي فركانسي و كمك گرفتن از قضيه ي نگاشت (Mapping) مي توان پايداري آن سيستم را بررسي كرد. با دستور زير و مطابق مثال هاي گذشته مي توان نمودار نايكوئيست يك تابع تبديل را رسم كرد و طبق معيار پايداري نايكوئيست، پايداري سيستم را بررسي نمود:
a=[0 1 2];
b=[1 2 3];
nyquist(a,b)
ravegoat
22-01-12, 16:08
با دستور زير مي توانيد تابع تبديل (Transfer Function) يك سيستم را به نمايش فضاي حالت (State Space) تبديل كنيد:
[A,B,C,D]=tf2ss(num,den)
كه num صورت تابع تبديل و den مخرج آن است.
همچنين اگر يك نمايش فضاي حالت به صورت:
xdot=Ax+Bu
y=Cx+Du
داشته باشيم (اين شكل نمايش در نظريه ي كنترل نوين بسيار كاربرد دارد)، با دستور زير مي توان تابع تبديل اين فضاي حالت را به دست آورد:
[num,den]=ss2tf(A,B,C,D)
ravegoat
25-01-12, 22:35
سورس متلب پیوست شده نمونه ای از روش تشخیص حروف انگلیسی دست نوشته (OCR) با شبکه ی عصبی است.
منبع (Only the registered members can see the link)
ravegoat
30-01-12, 12:28
در چرخدنده هاي مارپيچ (Helical) توان به تدريج روي دنده ها اثر مي كند و همين امر باعث كاهش ضربه و كاهش سر و صدا مي مي شود. فرآيند طراحي يك جعبه دنده ي هليكال (Helical Gearbox) شايد كمي زمان بر باشد. به علاوه چنين گيربكسي با تغيير نسبت هاي تبديل (Transmission Ratio) و مدول (Module) چرخدنده ها نياز به بهينه سازي (Optimization) دارد. سورس پيوست شده ي متلب در اين پست، با دريافت پارامتر هاي ورودي نظير نسبت تبديل كلي، توان ورودي و جنس چرخدنده ها، پارامتر هاي طراحي يك گيربكس هليكال دو مرحله اي را در اختيار كاربر قرار مي دهد و براساس معيار حجم گيربكس، آن را بهينه مي كند. لازم به ذكر است كه اين سورس تحليل خمش و تحليل خستگي سطحي را پوشش مي دهد.
ravegoat
30-01-12, 23:05
روش سيمپسون يكي از روش هاي متداول براي انتگرال گيري عددي است. براي استفاده از اين روش در انتگرال يگانه ي معين در متلب مي توان از دستور quad استفاده كرد:
F = @(x)exp(x.^2).*(cos(x).^3);
Q = quad(F,-0.1,1.4)
مثال فوق انتگرال تابع exp(x^2)*cos(x)^3 را بازه ي -0.1 تا 1.4 محاسبه مي كند.
ravegoat
06-02-12, 09:29
سورس متلب پيوست شده با دريافت ماتريس A و B در دستگاه معادلات:
AX=B
ابتدا ماتريس ضرايب-A-را به صورت حاصل ضرب يك ماتريس پايين مثلثي (Lower Triangular Matrix) و يك ماتريس بالا مثلثي (Upper Triangular Matrix) نمايش مي دهد و سپس به كمك اين ماتريس هاي مثلثي جواب دستگاه را حساب مي كند. همچنين اين سورس دستگاه را بار ديگر به روش ماتريس معكوس و با استفاده از دستور inv حل مي كند:
A-1(AX=B) => X=A-1B
و نتيجه را با روش قبلي مقايسه خواهد كرد.
ravegoat
08-02-12, 17:08
دستور fzero در متلب در حل عددی معادله ی f(x)=0 کاربرد دارد. فرضا" اگر تابع f به صورت x2 - 1 باشد آن گاه تابع بر حسب متغیر x به شکل زیر تعریف می شود:
f=@(x)x.^2-1
سپس یک جواب تقریبی معادله را حدس می زنیم؛ فرضا" 2 . آن گاه جواب معادله به شکل زیر به دست می آید:
fzero(f,2)
ravegoat
16-03-12, 11:49
Simulink® 7
Getting Started Guide
از MathWorks
در 93 صفحه به زبان انگليسي
ravegoat
24-06-12, 11:08
سورس پيوست شده مربوط به بهينه سازي چند هدفه با NSGA II است. كاربر مي تواند با تعريف دو تابع هدف (Cost Function) در فايل cf.m و تعيين پارامتر هايي نظير بازه ي تغييرات متغير هاي طراحي (به تعداد دلخواه)، طول كروموزوم، تعداد جمعيت اوليه و تعداد نسل ها در فايل main.m ، فرآيند بهينه سازي را اجرا نمايد. در نهايت برنامه Pareto Front به علاوه نقطه ي مصالحه ي طراحي (Trade-Off Point) را به كاربر نمايش مي دهد.
لازم به ذكر است كه اين سورس شامل بهينه سازي مقاوم (Robust Optimization) با ديدگاه هاي RBDO و RDO نيز است.
ravegoat
25-06-12, 20:06
این سورس با استفاده از مدل Single Track یک خودرو و محاسبه ی ضریب انتقال بار (LTR) ، واژگونی (RollOver) خودرو را پیش بینی می کند.
ravegoat
15-02-13, 19:07
در این بسته می توانید سورس های مفیدی را پیرامون نحوه ی طراحی فیلتر های سیگنال با متلب بیابید.
منبع (Only the registered members can see the link)
ravegoat
20-02-13, 19:29
عنوان:
System Identification Toolbox 7
Reference
نویسنده:
Lennart Ljung
ناشر:
MathWorks
این کتاب مرجع کاملی برای جعبه ابزار System Identification در Simulink نرم افزار متلب محسوب می شود که در 365 صفحه به زبان انگلیسی تدوین شده است.
ravegoat
28-02-13, 13:33
این آموزش خلاصه شده نحوه ی به کار گیری کتابخانه ی Open CV را در نرم افزار MATLAB نشان می دهد.
ravegoat
02-03-13, 13:57
فایل پیوست شده مربوط به آموزش طراحی Filter دیجیتال با روش Butterworth در نرم افزار MATLAB است.
arnouri
02-03-13, 14:45
دوست عزیز با سلام و تشکر از m-file خوبتون
یک درخواست داشتم: امکان داره لطف کنین و توضیح بدین برای بیشتر از 2 هدف چطور میشه از برنامه استفاده نمود؟
ravegoat
02-03-13, 20:00
دوست عزیز با سلام و تشکر از m-file خوبتون
یک درخواست داشتم: امکان داره لطف کنین و توضیح بدین برای بیشتر از 2 هدف چطور میشه از برنامه استفاده نمود؟
با سلام!
خواهش می کنم...
دوست گرامی؛ NSGA-II فقط برای دو تا تابع هدف می تونه متغیر های طراحی رو بهینه کنه و اگر تعداد توابع هدف تون بیش تر از دو تا باشه باید از SPEA و یا MUGA استفاده کنید که متاسفانه فعلا" سورس شون رو ندارم.
موفق باشید
آرمین
ravegoat
28-03-13, 16:21
کد متلب زیر الگوریتمی ساده و موثر را با پیاده سازی روش Total Variation جهت کاهش نویز سیگنال دیجیتال یک بعدی معرفی می کند:
function [x,J] = denoiseTV(y,lambda,Nit)
% [x,J] = denoiseTV(y,lambda,a,Nit)
% Total variation filtering (denoising) using
% iterative clipping algorithm.
% INPUT
% y - noisy signal (row vector)
% lambda - regularization parameter
% Nit - number of iterations
% OUTPUT
% x - result of denoising
% J - objective function
J = zeros(1,Nit); % objective function
N = length(y);
z = zeros(1,N-1); % initialize z
alpha = 4;
T = lambda/2;
for k = 1:Nit
x = y - [-z(1) -diff(z) z(end)]; % y - D’ z
J(k) = sum(abs(x-y).^2) + lambda * sum(abs(diff(x)));
z = z + 1/alpha * diff(x); % z + 1/alpha D z
z = max(min(z,T),-T); % clip(z,T)
end
منبع (Only the registered members can see the link)
bmaze
19-05-13, 09:28
سلام
ببخشید من یه برنامه میخواستم بنویسم سوالش اینه
برنامه ای بنویسید که یک عدد صحیح و یک عدد طبیعی کوچکتر از 9 را به عنوان مبنا از کاربر دریافت کند و عدد خوانده شده ی اولیه را به مبنای داده شده ببرد و نتیجه را در خروجی نمایش دهد.
میشه کمکم کنید؟
bmaze
19-05-13, 13:28
جواب سوالما پیدا کردم میزارم که بقیه هم استفاده کنن
clc;
clear all;
a=input('adad=');
b=input('mabna=');
i=1;
while a>=b
c(i)=mod(a,b);
a=floor(a/b);
i=i+1;
end
z=fliplr(c);
fprintf('%g%g',a,z);
behzadk777
20-05-13, 23:21
سلام . خسته نباشید
یکی از سورس ها به درد من خورد در این تاپیک
میخواستم ببینم مهندس شما سورس مربوط به ارتعاشات آزاد سیستم های یک درجه آزادی میرا دارید؟؟؟
اگر کمک کنید خیلی لطف کردید
ravegoat
21-05-13, 15:41
سلام . خسته نباشید
یکی از سورس ها به درد من خورد در این تاپیک
میخواستم ببینم مهندس شما سورس مربوط به ارتعاشات آزاد سیستم های یک درجه آزادی میرا دارید؟؟؟
اگر کمک کنید خیلی لطف کردید
با سلام!
عضویت تون رو در شهر سخت افزار تبریک می گم.
سورس نمونه ی ارتعاشات یک درجه آزادی زیر میرا (Underdamped) پیوست شد.
موفق باشید
آرمین
ravegoat
31-05-13, 10:38
این مجموعه اسلاید پیوست شده تحت عنوان " چکیده ای از مبانی کامپیوتر ها، الگوریتم ها و برنامه نویسی به کمک MATLAB " که توسط دکتر محمودی مهر از گروه مهندسی مکانیک دانشگاه گیلان گردآوری شده می تونه یه منبع مناسب برای کسانی باشه که دوست دارن برنامه نویسی رو با زبان متلب آغاز کنن.
ravegoat
11-07-13, 08:38
کد زیر پاسخ ضربه ی واحد (Impulse Response) و پاسخ پله ی واحد (Step Response) مربوط به یک تابع تبدیل دلخواه (Transfer Function) را رسم می کند:
s=tf('s');
g = 10 /(s^2 + 2*s +10); %Transfer Function
subplot(2,1,1)
impulse(g) %Plot Impulse Response
subplot(2,1,2)
step(g) %Plot Step Response
ravegoat
11-07-13, 23:12
فرمان lsim برای شبیه سازی پاسخ هر نوع تابع تبدیلی نسبت به یک ورودی دلخواه در فضای زمان کاربرد دارد. کد زیر پاسخ تابع تبدیل مثال قبل را نسبت به دو ورودی شیب و سهموی از زمان صفر تا 10 ثانیه با گام زمانی 0.1 ثانیه رسم می کند:
s=tf('s');
g = 10 /(s^2 + 2*s +10); %Transfer Function
t = 0:0.1:10; %0:dt:Tfinal
subplot(1,2,1)
u = t; %Ramp Input
lsim(g,u,t);
subplot(1,2,2)
u = t .^ 2; %Parabolic Input
lsim(g,u,t);
ravegoat
13-07-13, 22:58
فرمان randperm(n) در متلب یک ماتریس سطری با n ستون ایجاد می کند که هر درآیه ی آن عدد تصادفی بین یک تا n است به طوری که مقدار هیچ دو درآیه با هم برابر نیست.
randperm(6)
خروجی:
3 2 6 4 1 5
ravegoat
19-07-13, 16:29
برای تصویر کردن یک رویه در صفحات مختصات می توان از تابع جالبی به اسم ShadowPlot استفاده کرد. البته این تابع به طور پیش فرض در متلب وجود ندارد. پس برای اجرای کد زیر ابتدا این تابع را از پیوست دانلود کنید و سپس MFile آن را در مسیر جاری نرم افزار متلب کپی کنید. کد زیر رویه ای مطابق شکل را ترسیم و سایه ی آن رویه را در صفحات X و Y ایجاد می کند:
8912
f = @(x,y) exp(-x.^2 -y.^(-2)).*(x.^2+y.^2);
[X,Y] = meshgrid(-10:0.5:10,-10:0.5:10);
surf(X,Y,f(X,Y))
xlim([-11,11])
ylim([-11,11])
shadowplot x
shadowplot y
منبع (Only the registered members can see the link)
ravegoat
27-07-13, 08:24
برای ایجاد تابع چگالی احتمالات (PDF) می توان از فرمان hist استفاده کرد که نمونه ای از کاربرد این فرمان جهت رسم PDF در زیر آورده شده است:
data = [1 2 3 3 4]; %# Sample data
xRange = 0:10; %# Range of integers to compute a probability for
N = hist(data,xRange); %# Bin the data
plot(xRange,N./numel(data)); %# Plot the probabilities for each integer
xlabel('Integer value');
ylabel('Probability');
منبع (Only the registered members can see the link)
ravegoat
20-08-13, 14:46
با بهره گیری از فرمان sprintf می توانید مقادیر متغیر های مختلف را به هنگام نمایش عبارت فرمان input نشان دهید:
for i =1:9
for j=1:9
tmp(i,j) = input(sprintf('Plz enter value of row %g and column %g: ', i, j));
end
end
ravegoat
22-08-13, 15:59
اگر تمایل دارید تا زمان اجرای بخشی از کد خود را در متلب بدانید، کافیست آن بخش از کد را مطابق الگوی زیر بین فرمان های tic و toc قرار دهید:
tic;
%Your Code Here
toc;
ravegoat
23-08-13, 19:44
برای آن که دریابید چقدر زمان صرف اجرای هر بخش از کد شما در متلب می شود، می توانید کد خود را مطابق الگوی زیر اجرا نمایید:
profile on
%Your Code Here
profile off
profile viewer
ravegoat
24-08-13, 08:12
سورس پیوست شده نحوه ی پیاده سازی الگوریتم های مختلف Edge Detection نظیر Prewitt، Sobel، Roberts و Canny را روی تصاویر نشان می دهد.
9171
ravegoat
02-09-13, 16:36
کد زیر یک ترفند ساده را برای تبدیل یک ماتریس به یک بردار نشان می دهد:
x=rand(4);
x=x(:)
ravegoat
11-10-13, 20:11
فرمان زیر مقادیر ویژه و بردار های ویژه ی (Eigenvalues and Eigenvectors) ماتریس A را محاسبه کرده و حاصل را به ترتیب در D و V قرار می دهد:
[V, D] = eig(A)
ravegoat
23-10-13, 16:51
کد زیر جهت تحلیل یک سیستم دینامیکی بیان شده در فضای حالت به کار می رود که بدون ورودی بوده و تنها تحت تاثیر شرایط اولیه است:
سلام دوستان ، من یه سوال دارم از اونایی که با برنامه نویسی Mathlab آشنا هستند ،
سوالم اینه چطور می تونم دو تا عدد 100 رقمی رو بر هم تقسیم کنم ؟
البته سوال اصلی من یه چیز دیگه است ولی تو قسمت بالا موندم ، فقط چون برنامه من خیلی طولانیه و مورد بالا هم خیلی زیاد تو برنامم استفاده میشه خواهشن جوابها کوتاه باشه .
فقط خواهشی که از دوستان دارم اینه که اگه به طور عملی جواب گرفتند برنامه رو اعلام کنند .
ravegoat
24-10-13, 22:05
سلام دوستان ، من یه سوال دارم از اونایی که با برنامه نویسی Mathlab آشنا هستند ،
سوالم اینه چطور می تونم دو تا عدد 100 رقمی رو بر هم تقسیم کنم ؟
البته سوال اصلی من یه چیز دیگه است ولی تو قسمت بالا موندم ، فقط چون برنامه من خیلی طولانیه و مورد بالا هم خیلی زیاد تو برنامم استفاده میشه خواهشن جوابها کوتاه باشه .
فقط خواهشی که از دوستان دارم اینه که اگه به طور عملی جواب گرفتند برنامه رو اعلام کنند .
با سلام!
کلا" برای اجرای عملیات روی اعداد بزرگ یه جعبه ابزار به نام VPI وجود داره که البته به طور پیش فرض با متلب نصب نمیشه و باید جداگونه اون رو نصب کرد. پایه ی کار این جعبه ابزار این طوره که اعداد بزرگ رو به شکل رشته در میاره و اعمال رو روی اون ها پیاده می کنه.
به شکل مشابه خودتون می تونید اعداد صد رقمی رو به صورت رشته (string) تعریف کنید و مطابق همون الگوی تقسیم که در دوران ابتدایی بهمون آموزش دادن رقم به رقم اعداد رو به هم تقسیم کنید.
اگر فرصت کنم کدش رو در همین جا قرار میدم.
موفق باشید
آرمین
علی حسنی
26-10-13, 18:06
سلام مجدد و تشکر فراوان ،
من خدمتتون عرض کردم چون برنامم خیلی طولانیه نمی تونم برای عمل تقسیم برنامه بنویسم (این برنامه حتی اگه چند خط باشه جواب برنامه اصلی که باید در کمتر از یک ثانیه جواب بده ساعتها طول خواهد کشید ! ) فقط می خوام یه برنامه یا نرم افزاری باشه که خودش دو متغیر 100 رقمی رو بخونه و جواب تقسیم یا ضربشون رو بده، تازه اعداد در پروژه من حداقل صد رقمی هستند و ممکنه صد ها رقم داشته باشه !
برنامه ها و نرم افزارهای مختلفی رو هم تست کردم ولی به جواب نرسیدم اگه باز هم کسی از دوستان نرم افزار یا راهکاری سراغ داره ممنون میشم راهنمایی کنه تا من این برنامه به ظاهر ساده رو حل کنم .
ravegoat
27-10-13, 20:40
اگر یک معادله ی چند جمله ای از مرتبه ی n به شکل زیر داشته باشیم:
anxn + an-1xn-1 + ... a1x + a0 = 0
با فرمان roots می توان ریشه های آن را پیدا کرد.
ابتدا یک ماتریس سطری از ضرایب چند جمله ای به شکل زیر شکل می دهیم:
p = [an an-1 ... a1 a0]
سپس با فرمان زیر ریشه ها محاسبه می شوند:
roots(p)
مثال: حل معادله
x3 + 2x2 + 2x + 1
p = [1 2 2 1];
roots(p)
ravegoat
03-11-13, 16:13
برای یافتن فرم Jordan Canonical یک ماتریس می توان از دستور jordan به شکل زیر استفاده کرد:
a=[1 2 3;3 0 1;0 1 2];
jordan(a)
ravegoat
03-11-13, 16:42
برای یافتن ماتریس انتقال حالت (State Transition Matrix) در متلب به شکل زیر عمل می کنیم:
a=[1 0;-1 -1];
syms t
expm(a*t)
ravegoat
04-11-13, 20:28
برای یافتن rank ماتریس دلخواه a به شکل زیر عمل می کنیم:
a=[1 1 0 0;-2 1 0 1;3 -3 -1 2]
rank(a)
ravegoat
05-11-13, 19:30
برای انتقال تابع f از فضای زمان به فضای لاپلاس به شکل زیر عمل می کنیم:
syms t
f = 1 - exp(t)
laplace(f)
ravegoat
05-11-13, 19:33
برای بردن g از حوزه ی لاپلاس به حوزه ی زمان به شکل زیر عمل می کنیم:
syms s
g = 1/s - 1/(s-1)
ilaplace(g)
ravegoat
09-11-13, 13:35
برای یافتن Sine و Cosine ماتریس دلخواه A به شکل زیر عمل می کنیم:
A = magic(3);
funm(A, sin)
funm(A, cos)
ravegoat
20-11-13, 09:11
تجزیه ی مقدار تکین یا به اختصار SVD یک روش تجزیه ی ماتریس های حقیقی و مختلط است که کاربرد های جالبی در مهندسی دارد. اگر M یه ماتریس دلخواه باشد، تجزیه ی مقدار تکین آن متشکل از سه ماتریس U، S و V خواهد بود که به شکل زیر تعریف می شود:
USV* = M
که در آن V* ترانهاده ی مزدوج ماتریس V است.
برای یافتن SVD ماتریس M در متلب به شکل زیر عمل می کنیم:
چگالی طیفی توانی یا به اختصار PSD ابزاری در پردازش سیگنال جهت یافتن نحوه ی توزیع قدرت سیگنال برحسب فرکانس است. نمونه ای از محاسبه ی PSD با MATLAB در زیر آمده است:
t = 0:.001:.25;
x = sin(2*pi*50*t) + sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(y(1:50))
title('Noisy time domain signal')
Y = fft(y,251);
Pyy = Y.*conj(Y)/251;
f = 1000/251*(0:127);
plot(f,Pyy(1:128))
title('Power spectral density')
xlabel('Frequency (Hz)')
ravegoat
24-11-13, 13:30
برای یافتن تبدیل Z یک تابع گسسته در MATLAB به شکل زیر عمل می کنیم:
syms n;
f = n^2;
ztrans(f)
ravegoat
24-11-13, 13:34
برای یافتن تبدیل معکوس Z در متلب به این شکل عمل می کنیم:
syms z
F = z/(z-1);
iztrans(F)
ravegoat
26-11-13, 09:13
برای یافتن تبدیل فوریه یک تابع در حوزه ی زمان می توان از دستور fourier استفاده کرد. مثال:
syms x;
f = exp(-x^2);
fourier(f)
ravegoat
28-11-13, 10:32
در این پست مدل Simulink کنترل یک سیستم جرم و فنر با روش State Vector Feedback قرار گرفته است. در SVFC با فیدبک گرفتن از تمام حالت های سیستم، نوع پاسخ آن کنترل می شود. اما در مدل جرم و فنر در عمل امکان فیدبک گرفتن از سرعت وجود ندارد. لذا جهت یافتن مقدار تقریبی سرعت باید از یک تخمین گر حالت (State Estimator) استفاده کرد. جهت اجرای این کار از یک مشاهده گر با رتبه ی کاهش یافته (Reduced Order Observer) بهره گرفته شده است.
ravegoat
02-12-13, 19:21
جهت تولید پالس مربعی در MATLAB از فرمان gensig به شکل زیر عمل می کنیم:
[u,t] = gensig('square',T,d,Ts)
که در آن:
T دوره تناوب پالس
d مدت زمان پالس
Ts زمان نمونه برداری
مثال:
[u,t] = gensig('square',5,30,0.1)
ravegoat
07-12-13, 15:00
اگر یک سیستم دینامیکی در فضای حالت به شکل زیر تعریف شده باشد:
x. = Ax + Bu
y = Cx + Du
آن گاه ماتریس کنترل پذیری و مشاهده پذیری به ترتیب به شکل زیر در به دست می آیند:
ctrb(A,B)
obsv(A,C)
مثال:
A = [1 2;-1 1];
B = [2;1];
C = [1 0];
ctrb(A,B)
obsv(A,C)
ravegoat
12-01-14, 15:00
برای محاسبه ی گرامیان کنترل پذیری و مشاهده پذیری در متلب از فرمان gram استفاده می شود.
جهت طراحی یک کنترلر بهینه ی مربعی-خطی در متلب می توان از فرمان LQR به شکل زیر بهره برد:
[K,S,e] = LQR(A,B,Q,R,N)
که در آن A و B ماتریس های فضای حالت سیستم و Q، R و N پارامتر های تابع هدف هستند. در آن صورت K ماتریس ضرایب بردار فیدبک حالت، S پاسخ معادله ی ریکاتی و e مقادیر ویژه ی سیستم مطلوب خواهد بود.
tarane66
26-01-14, 00:41
من به کد متلب روش ضمنی دو بعدی در مکعب مستطیل نیاز دارم. لطفا کمکم کنید.
ravegoat
26-01-14, 08:34
من به کد متلب روش ضمنی دو بعدی در مکعب مستطیل نیاز دارم. لطفا کمکم کنید.
دوست گرامی عضویت شما رو در شهر سخت افزار تبریک می گم.
اگر منظورتون روش ضمنی تفاضل محدود در حل معادلات دیفرانسیل کارتزین دو بعدی هستش می تونید به دو پیوند زیر مراجعه کنید:
Parallel Spectral Numerical Methods/Examples in Matlab - Wikibooks, open books for an open world (Only the registered members can see the link)
Math467 - Scientific Computation (Only the registered members can see the link)
ولی اگه منظورتون چیز دیگه ای بوده، لطفا سوالتون رو واضح تر مطرح بفرمایید.
با سپاس
آرمین
ravegoat
05-02-14, 23:29
الگوریتم Local Linear Model Tree یا به اختصار LOLIMOT یک الگوریتم فازی (Only the registered members can see the link)-عصبی (Only the registered members can see the link) براساس مدل TSK است که در یادگیری سیستم های غیر خطی و تشخیص الگو کاربرد دارد.
ravegoat
18-02-14, 15:33
برای یافتن تبدیل فوریه ی دو بعدی یک تصویر می توان به شکل زیر عمل کرد:
I = imread('test.jpg'); % Read image from the path
F = fft2(I); % 2D FFT
F = fftshift(F); % Center FFT
F = abs(F); % Get the magnitude
F = log(F+1); % Use log, for perceptual scaling, and +1 since log(0) is undefined
F = mat2gray(F); % Use mat2gray to scale the image between 0 and 1
imshow(F,[]); % Display the result
منبع (Only the registered members can see the link)
ravegoat
24-02-14, 14:17
با دستور diag به شکل زیر می توان عناصر روی قطر اصلی یک ماتریس مربعی را استخراج کرد:
a = rand(3);
diag(a)
ravegoat
06-04-14, 08:21
سورس پیوست شده یک آشوب را در سیستم دینامیکی را نشان می دهد که به ازای دو شرط اولیه ی نزدیک به هم، پاسخ های سیستم پس از مدتی تفاوت قابل ملاحظه ای از هم پیدا می کنند.
ravegoat
10-04-14, 14:33
دستور margin فرمانی مفید جهت یافتن حد بهره (Gain Margin) و حد فاز (Phase Margin) یک سیستم دینامیکی است. مثال:
hd = tf([0.04798 0.0464],[1 -1.81 0.9048],0.1)
[Gm, Pm] = margin(hd)
ravegoat
13-04-14, 07:18
جهت ایجاد یک عدد تصادفی با توزیع نرمال یا گاوسی از فرمان randn در متلب استفاده می شود. مثال زیر یک عدد تصادفی را براساس یک توزیع نرمال با میانگین 1 و انحراف معیار 2 ایجاد می کند:
r = 1 + 2.*randn(1)
M A H R A D
13-04-14, 12:54
جهت ایجاد یک عدد تصادفی با توزیع نرمال یا گاوسی از فرمان randn در متلب استفاده می شود. مثال زیر یک عدد تصادفی را براساس یک توزیع نرمال با میانگین 1 و انحراف معیار 2 ایجاد می کند:
r = 1 + 2.*randn(1)
آرمین جان برای توزیع های دیگه مثل نمایی یا برنولی هم دستوری وجود داره؟
ravegoat
13-04-14, 14:51
آرمین جان برای توزیع های دیگه مثل نمایی یا برنولی هم دستوری وجود داره؟
بله مهراد جان... فرمان random در آرگومان خودش اسم توزیع رو دریافت می کنه که شامل توزیع های متنوعی میشه. مثلا" برای تولید یک عدد تصادفی با توزیع ریلی فرمان زیر رو وارد می کنیم:
x1 = random('rayl',1)
برای برخی از توزیع های پرکاربرد هم از جمله همین توزیع نمایی و توزیع برنولی دستور اختصاصی وجود داره که به ترتیب exprnd و binornd هستش.
موفق باشی
ravegoat
24-04-14, 10:36
برای ذخیره کردن محتویات پنجره ی Figure در یک فایل از فرمان saveas به شکل زیر استفاده می کنیم:
h = figure(1);
saveas(h, 'filePath', 'jpg')
sitara
03-05-14, 13:36
سلام. ببخشید میخواستم بدونم کد آماده ای در متلب هست که بتونه یک ماتریس رو به عنوان ورودی بگیره و درایه های اون رو به طور اتوماتیک تغییر بده و این کار رو چندین بار انجام بده و در هر بار یک ماتریس جدید به ما بده (مثلا با 100 تکرار) . خیلی ممنون.
ravegoat
03-05-14, 19:26
سلام. ببخشید میخواستم بدونم کد آماده ای در متلب هست که بتونه یک ماتریس رو به عنوان ورودی بگیره و درایه های اون رو به طور اتوماتیک تغییر بده و این کار رو چندین بار انجام بده و در هر بار یک ماتریس جدید به ما بده (مثلا با 100 تکرار) . خیلی ممنون.
با سلام!
دوست گرامی عضویت شما رو در شهر سخت افزار تبریک می گم.
منظورتون از این که " درآیه های ماتریس رو به طور اتوماتیک تغییر بده " چیه؟
ravegoat
06-05-14, 16:47
برای پیچاندن (Convolve) دو سیگنال تک بعدی x و y در هم از فرمان conv به شکل زیر استفاده می کنیم:
conv(x,y)
ravegoat
12-05-14, 13:18
سورس پیوست شده مربوط به تشخیص چهره با الگوریتم Kanade–Lucas–Tomasi یا به اختصار KLT است که برخلاف الگوریتم Voila-Jones قابلیت تعیین زاویه ی چرخش چهره و یا تجانس آن را دارد.
ravegoat
18-05-14, 19:41
کد زیر مثالی از نوشتن مقدار سه متغیر در یک فایل متنی است:
a = 2;
b = 4;
c = 10;
fid = fopen('output.txt','w');
fprintf(fid,'a=%g\r\n',a);
fprintf(fid,'b=%g\r\n',b);
fprintf(fid,'c=%g',c);
fclose(fid);
ravegoat
24-05-14, 19:57
برای تبدیل یک رنگ از فضای RGB به HSV به شکل زیر عمل می کنیم:
rgb2hsv([255 128 0])
ravegoat
26-05-14, 21:25
کد MATLAB پیوست شده تابعی است که کاربر را قادر می سازد تا فریم های دلخواهی را از یک فایل ویدئویی حذف کند و ویدئوی جدید را در خروجی تولید نماید.
ravegoat
18-06-14, 19:56
جهت Crop کردن یک تصویر در متلب به شکل زیر عمل می کنیم:
که در آن x1، y1 و x2، y2 مختصات دو گوشه ی مستطیل برش اند.
ravegoat
24-06-14, 21:42
برای تولید اعداد تصادفی صحیح بین دو عدد a و b به شکل زیر عمل می کنیم:
randi([a,b])
ravegoat
07-07-14, 23:19
فایل پیوست شده مربوط به شبیه سازی یک Classic Washout Filter جهت ایجاد حرکات شتاب دار در شبیه ساز ها است. برای اجرای فایل نیاز به نرم افزار MATLAB 2013 است.
gomnam15010
19-07-14, 22:58
در صورتی که بخواهیم بهینه سازی مقاوم انجام دهیم چگونه باید از این کد استفاده کرد؟
سورس پيوست شده مربوط به بهينه سازي چند هدفه با NSGA II است. كاربر مي تواند با تعريف دو تابع هدف (Cost Function) در فايل cf.m و تعيين پارامتر هايي نظير بازه ي تغييرات متغير هاي طراحي (به تعداد دلخواه)، طول كروموزوم، تعداد جمعيت اوليه و تعداد نسل ها در فايل main.m ، فرآيند بهينه سازي را اجرا نمايد. در نهايت برنامه Pareto Front به علاوه نقطه ي مصالحه ي طراحي (Trade-Off Point) را به كاربر نمايش مي دهد.
لازم به ذكر است كه اين سورس شامل بهينه سازي مقاوم (Robust Optimization) با ديدگاه هاي RBDO و RDO نيز است.
ravegoat
20-07-14, 21:05
در صورتی که بخواهیم بهینه سازی مقاوم انجام دهیم چگونه باید از این کد استفاده کرد؟
با سلام!
در مجموعه ای که شما بهش اشاره کردید بهینه سازی مقاوم بر اساس دو رویکرد RBDO و RDO صورت می گیره. در RBDO توابع هدف شما احتمال های شکست خواهند بود که در فایلی که در انجمن قرار داده شده بود احتمال شکست در مورد مقاومت و جرم یک تیر تحت بار لحاظ شده بود. در RDO میانگین و انحراف مربوط به یک پارامتر مثلا مقاومت تیر به عنوان توابع هدف مطرح میشه. برای استفاده از این فایل ها شما کافیه تعداد نمونه ها و PDF و CDF پارامتر های عدم قطعیت رو تعیین کنید و خروجی رو به عنوان توابع هدف در NSGA استفاده کنید.
موفق باشید
آرمین
ravegoat
26-07-14, 19:19
برای یافتن عدد شرطی (Condition Number) ماتریس A به شکل زیر عمل می کنیم:
cond(A)
ravegoat
20-08-14, 19:16
کد زیر محتویات داخل حلقه را تنها تا 10 ثانیه تکرار می کند:
start_tm = clock();
loopTime = 10;
while etime(clock(),start_tm)<loopTime,
%... Put your code in here ...
end
منبع (Only the registered members can see the link)
ravegoat
04-09-14, 19:34
برای یافتن تابع تبدیل حلقه بسته ی (Closed Loop) سیستم مدار باز G که تابع تبدیل مسیر فیدبک آن H است، می توان از فرمان feedback به شکل زیر استفاده کرد:
feedback(G,H)
ravegoat
12-09-14, 11:26
برای استخراج هیستوگرام روشنایی یک تصویر در متلب به شکل زیر عمل می کنیم:
برای یافتن مقدار میانه در یک آرایه از فرمان median به شکل زیر استفاده می کنیم:
median([1 5 3 10 6])
ans = 5
JiMmY.j
21-10-14, 09:39
سلام دوست عزیز.
من خیلی فوری یه برنامه متلب می خوام که 1000 تا تابع سینوسی با فازهای مرتب ( با فاصله های یکسان) رو ایجاد کنه. بعد این توابع رو بر حسب زمان رسم کنه.
مرسی.:give_rose:
JiMmY.j
21-10-14, 10:28
من خودم این برنامه رو نوشتم شما اطلاع داردید درست هست یا نه؟
سلام مهندس
من میخوام تبدیل فوریه یک پالس مربعی رو رسم کنم.
کمک میکنید؟
ravegoat
22-10-14, 20:47
سلام دوست عزیز.
من خیلی فوری یه برنامه متلب می خوام که 1000 تا تابع سینوسی با فازهای مرتب ( با فاصله های یکسان) رو ایجاد کنه. بعد این توابع رو بر حسب زمان رسم کنه.
مرسی.:give_rose:
با سلام!
دوست گرامی منظورتون از فاز مرتب چیه؟ اگر منظورتون اینه که 1000 تا تابع سینوسی مد نظر تون است که هر کدوم با تابع قبلی یک اختلاف فاز ثابت داشته باشن کدی که قرار دادید رو باید به شکل زیر اصلاح بفرمایید:
clear all
clf
clc
t=linspace(0,10);
f=linspace (0,2*pi,1000);
for i=1:1000
y=sin((2*pi*t)+f(i));
hold on
plot(t,y);
end
title('Constant Phase & Constant Amp');
ylabel('Summation');
xlabel('Time');
موفق باشید
آرمین
ravegoat
22-10-14, 21:00
سلام مهندس
من میخوام تبدیل فوریه یک پالس مربعی رو رسم کنم.
کمک میکنید؟
با سلام!
برای یافتن بیش ترین مقدار در یک آرایه و همچنین اندیس آن مقدار در MATLAB می توان به شکل زیر عمل کرد:
A = [1 3 5 2 4];
[maxVal,maxInd] = max(A);
خروجی
maxVal = 5
maxInd = 3
JiMmY.j
24-10-14, 17:18
با سلام!
دوست گرامی منظورتون از فاز مرتب چیه؟ اگر منظورتون اینه که 1000 تا تابع سینوسی مد نظر تون است که هر کدوم با تابع قبلی یک اختلاف فاز ثابت داشته باشن کدی که قرار دادید رو باید به شکل زیر اصلاح بفرمایید:
clear all
clf
clc
t=linspace(0,10);
f=linspace (0,2*pi,1000);
for i=1:1000
y=sin((2*pi*t)+f(i));
hold on
plot(t,y);
end
title('Constant Phase & Constant Amp');
ylabel('Summation');
xlabel('Time');
موفق باشید
آرمین
مرسی آرمین جان. ولی نتیجه عجیب و غریب شد.:1. (27):
میشه بفرمایید هدف از این دستور Hold On چیه و چرا دستور Plot داخل حلقه قرار داده شده؟
ممنون.
ravegoat
24-10-14, 17:57
مرسی آرمین جان. ولی نتیجه عجیب و غریب شد.:1. (27):
میشه بفرمایید هدف از این دستور Hold On چیه و چرا دستور Plot داخل حلقه قرار داده شده؟
ممنون.
خواهش می کنم...
نتیجه عجیب نیست. چون اختلاف فاز توابع کوچک در نظر گرفته شده، وقتی همه ی این توابع رو در یک صفحه رسم کنیم به چنین شکلی می رسیم. اگر اختلاف فاز ها افزایش پیدا کنه، توابع به شکل تفکیک شده رسم میشن.
دستور hold on برای اینه که وقتی نمودار جدیدی در داخل پنجره ی Figure رسم شد، نمودار قبلی پاک نشه و نمودار جدید روی اون رسم بشه (همین دستور باعث شده بود که خروجی کد قبلی این طور بشه).
حلقه هم برای این بود که مقدار فاز رو در هر گام به میزان ثابتی افزایش بده.
باز هم اگه سوالی بود در خدمتم.
آرمین:11():
JiMmY.j
24-10-14, 18:02
خواهش می کنم...
نتیجه عجیب نیست. چون اختلاف فاز توابع کوچک در نظر گرفته شده، وقتی همه ی این توابع رو در یک صفحه رسم کنیم به چنین شکلی می رسیم. اگر اختلاف فاز ها افزایش پیدا کنه، توابع به شکل تفکیک شده رسم میشن.
دستور hold on برای اینه که وقتی نمودار جدیدی در داخل پنجره ی Figure رسم شد، نمودار قبلی پاک نشه و نمودار جدید روی اون رسم بشه (همین دستور باعث شده بود که خروجی کد قبلی این طور بشه).
حلقه هم برای این بود که مقدار فاز رو در هر گام به میزان ثابتی افزایش بده.
باز هم اگه سوالی بود در خدمتم.
آرمین:11():
بازهم متشکر.
می خواستم ببینم که فرق خروجی کدی که شما زحمت کشیدید و کدی که بنده نوشته بودم در چی هست. به طور مثال اگه من به جای f(i) نوشته باشم f و از حلقه for هم استفاده نکرده باشم چه اتفاقی تو محاسبات متلب میوفته؟ قبل از این فکر می کردم که این دو کد یکی هستند.
ممنون از پاسخ های زیباتون....
ravegoat
24-10-14, 19:28
بازهم متشکر.
می خواستم ببینم که فرق خروجی کدی که شما زحمت کشیدید و کدی که بنده نوشته بودم در چی هست. به طور مثال اگه من به جای f(i) نوشته باشم f و از حلقه for هم استفاده نکرده باشم چه اتفاقی تو محاسبات متلب میوفته؟ قبل از این فکر می کردم که این دو کد یکی هستند.
ممنون از پاسخ های زیباتون....
همچنان خواهش می کنم...
خیر این دو کد معادل هم نیستند:
شما فرموده بودید که باید 1000 تابع سینوسی داشته باشیم که با هم اختلاف فاز مرتب داشته باشند. بنده این گونه برداشت کردم که ابتدا یک تابع سینوسی با اختلاف فاز فرضا" phi بین زمان صفر تا ده ثانیه رسم می شود. در گام بعدی همین تابع با اختلاف فاز 2phi رسم می شود و این روند تا مرحله ی 1000 ادامه می یابد. در نهایت ما 1000 تابع سینوسی خواهیم داشت.
در کد شما ابتدا سینوس یک نقطه در زمان 0.1 ثانیه با اختلاف فاز phi رسم می شود. در زمان 0.2 ثانیه سینوس آن نقطه با اختلاف فاز 2phi رسم می شود و این روند تا زمان 10 ادامه می یابد. در نهایت ما 1000 نقطه خواهیم داشت که تنها یک تابع سینوسی را شکل می دهند.
اگر جایی ابهام داشت بفرمایید تا بیش تر توضیح بدم.
آرمین:11():
JiMmY.j
24-10-14, 19:53
همچنان خواهش می کنم...
خیر این دو کد معادل هم نیستند:
شما فرموده بودید که باید 1000 تابع سینوسی داشته باشیم که با هم اختلاف فاز مرتب داشته باشند. بنده این گونه برداشت کردم که ابتدا یک تابع سینوسی با اختلاف فاز فرضا" phi بین زمان صفر تا ده ثانیه رسم می شود. در گام بعدی همین تابع با اختلاف فاز 2phi رسم می شود و این روند تا مرحله ی 1000 ادامه می یابد. در نهایت ما 1000 تابع سینوسی خواهیم داشت.
در کد شما ابتدا سینوس یک نقطه در زمان 0.1 ثانیه با اختلاف فاز phi رسم می شود. در زمان 0.2 ثانیه سینوس آن نقطه با اختلاف فاز 2phi رسم می شود و این روند تا زمان 10 ادامه می یابد. در نهایت ما 1000 نقطه خواهیم داشت که تنها یک تابع سینوسی را شکل می دهند.
اگر جایی ابهام داشت بفرمایید تا بیش تر توضیح بدم.
آرمین:11():
بله درست می فرمایید. متشکر.
من به مباحث مربوط به Image Processing علاقه زیادی دارم. بخش های مقدماتی و پایه اون رو هم با کمک استادم کار کردم. ولی به تعدادی منابع برای بخش های پیشرفته نیاز دارم. شما می تونید منابع دیگه ای رو معرفی کنید. انگلیسی یا فارسی باشه فرقی نمی کنه.
مرسی...
:11()::11()::11():
ravegoat
25-10-14, 15:19
بله درست می فرمایید. متشکر.
من به مباحث مربوط به Image Processing علاقه زیادی دارم. بخش های مقدماتی و پایه اون رو هم با کمک استادم کار کردم. ولی به تعدادی منابع برای بخش های پیشرفته نیاز دارم. شما می تونید منابع دیگه ای رو معرفی کنید. انگلیسی یا فارسی باشه فرقی نمی کنه.
مرسی...
:11()::11()::11():
مباحث پردازش تصویر رو می تونید در این تاپیک (Only the registered members can see the link) دنبال کنید. پست نهم این تاپیک حاوی یک کتاب است که به طور پیشرفته به مباحث بینایی ماشین پرداخته. اگر سوالی هم داشتید می تونید در همون جا مطرح کنید.
موفق باشید
آرمین
ravegoat
20-11-14, 17:53
اگر بخواهیم در MATLAB ردیف های یک ماتریس را بر حسب مقادیر یک ستون آن مرتب کنیم از دستور sortrows به شکل زیر استفاده می کنیم:
A = magic(5)
sortrows(A, 4)
کد فوق ماتریس A را برحسب ستون 4 ام آن از کوچک به بزرگ مرتب می کند. اگر قصد داشتیم که ماتریس A بر حسب همان ستون از بزرگ به کوچک مرتب شود، به جای 4 از -4 استفاده می کردیم.
ravegoat
12-02-15, 15:53
کد زیر مثالی از خواندن اطلاعات از درگاه سریال توسط MATLAB است:
منبع (Only the registered members can see the link)
ravegoat
28-04-15, 23:22
متغیر alpha در کد زیر شامل میزان شفافیت پیکسل ها در موقعیت متناظر در تصویر Image.png خواهد بود:
[A, map, alpha] = imread('Image.png')
ravegoat
29-12-15, 08:02
فایل پیوست شده در بر گیرنده ی آموزش جبر خطی، جبر ماتریسی و عملگر های آماری در MATLAB همراه با مثال های کاربردی است که می تواند برای شروع راهنمای بسیار مناسبی باشد.
mehdi jafari h
12-02-16, 17:49
سلام دوست عزیز
کد مربوط به ارتعاش سیم را نیازدارم
ممنون
- - - Updated - - -
سلام دوست عزیز
کد مربوط به ارتعاش سیم را نیازدارم
ممنون
ravegoat
12-02-16, 19:26
سلام دوست عزیز
کد مربوط به ارتعاش سیم را نیازدارم
ممنون
- - - Updated - - -
سلام دوست عزیز
کد مربوط به ارتعاش سیم را نیازدارم
ممنون
با سلام،
دوست گرامی به شهر سخت افزار خوش آمدید.
به پیوند های زیر مراجعه فرمایید:
The vibrating string example (Only the registered members can see the link)
MATLAB GUIs - Vibrating String (Only the registered members can see the link)
Only the registered members can see the link
Slingshot Effect (Only the registered members can see the link)
موفق باشید
آرمین
فرین
22-11-16, 09:41
سلام دوستان ممنون میشم برای کد متلب این مساله که مریوط به درس یادگیری ماشین هوش مصنوعی هست کمکم کنید
یک صفحه که ده تا نقطه در اون هست که یک خط به این۱۰ تا نقطه فیت کنیم رد بشه و ماشین خطای فیتنس میده اگر نبود برمیداریم دو تا خط فیت میکنیم دوباره خطا میده اگر نداشت برمیداریم اینبار سه خط و بعد از اون یک منحنی رسم میکنیم و اینکه چه اتفافی افتاده :-(
ravegoat
26-11-16, 08:40
سلام دوستان ممنون میشم برای کد متلب این مساله که مریوط به درس یادگیری ماشین هوش مصنوعی هست کمکم کنید
یک صفحه که ده تا نقطه در اون هست که یک خط به این۱۰ تا نقطه فیت کنیم رد بشه و ماشین خطای فیتنس میده اگر نبود برمیداریم دو تا خط فیت میکنیم دوباره خطا میده اگر نداشت برمیداریم اینبار سه خط و بعد از اون یک منحنی رسم میکنیم و اینکه چه اتفافی افتاده :-(
با سلام،
دوست گرامی به شهر سخت افزار خوش آمدید :1. (35):.
متاسفانه بنده به درستی متوجه عبارات برجسته شده در سوالتون نشدم:
اگر چی نبود؟
آیا منظور دو تا خط در ادامه ی هم است یا دو خط مجزا؟
اگر چی نداشت؟
چه نوع منحنی ای؟ منحنی نمایی؟
منظور کدوم اتفاقه؟
اگر لطف بفرمایید ابهامات فوق را رفع کنید، بهتر میشه راه حل ارایه کرد.
mhm233
22-01-17, 17:41
کسی میتونه لطف کنه سورس مربوط به رویه های درجه دوم رو اینجا بزاره ؟ متلب رو سیستم من کرک نمیشه خودم نمیتونم امتحان کنم فکر نمیکنم کدش پیچیده باشه!؟!
x^2/a^2)+(y^2/b^2)+(z^2/c^2)=1)
x^2/a^2)+(y^2/b^2)-(z^2/c^2)=1)
x^2/a^2)-(y^2/b^2)+(z^2/c^2)=1)-
x^2/a^2)+(y^2/b^2)=z/c)
x^2/a^2)+(y^2/b^2)=z/c)-
2^x^2/a^2)+(y^2/b^2)=z^2/c)
ravegoat
23-01-17, 07:19
کسی میتونه لطف کنه سورس مربوط به رویه های درجه دوم رو اینجا بزاره ؟ متلب رو سیستم من کرک نمیشه خودم نمیتونم امتحان کنم فکر نمیکنم کدش پیچیده باشه!؟!
x^2/a^2)+(y^2/b^2)+(z^2/c^2)=1)
x^2/a^2)+(y^2/b^2)-(z^2/c^2)=1)
x^2/a^2)-(y^2/b^2)+(z^2/c^2)=1)-
x^2/a^2)+(y^2/b^2)=z/c)
x^2/a^2)+(y^2/b^2)=z/c)-
2^x^2/a^2)+(y^2/b^2)=z^2/c)
با سلام،
می تونید با تابع ezsurf (Only the registered members can see the link) این کار رو انجام بدید ولی باید دو نکته رو مد نظر داشته باشید:
تابع باید به شکل صریح تعریف بشه یعنی z=f(x,y)
ضابطه نباید دارای مقادیر پارامتری نظیر a، b و c باشه!
اگر اصرار به استفاده از مقادیر پارامتری دارید، مطابق آموزش زیر عمل کنید:
ezsurf (Symbolic Math Toolbox) (Only the registered members can see the link)
Sr._.sahar
15-12-17, 16:15
با سلام و خسته نباشید
لطفا اگه امکان دارد در خصوص برنامه زیر راهنمایی کنید . ممنون
برنامه‌ای بنویسید که نمرات چند دانشجو را به صورت یک بردار بگیرد و عملیات زیر را انجام دهد:
در صورتیکه ورودی کاربر بردار نباشد(ماتریس یا اسکالر باشد) پیام خطا دهد.(راهنمایی برای دادن پیام خطا می‌توانید از تابع error به جای disp استفاده کنید)
با استفاده از حلقه for و دستورات شرطی if-else-end تک تک نمرات را چک کند و به صورت زیر آنها را تغییر دهد:
نمرات کمتر از ۵ را به ۹ تغییر دهد
نمرات بین ۵ و ۸ را به ۹.۵ تغییر دهد.
نمرات بین ۸ و ۱۰ را به ۱۰ تغییر دهد.
نمرات بین ۱۰ و ۱۵ را ۱ نمره افزایش دهد
نمرات بیشتر از ۱۵ و کمتر از ۲۰ را ۰.۵ نمره افزایش دهد.
sepehrdad.moradi
16-07-18, 11:56
سلام.راستش من نتونستم هیچکدوم از فایهای این تاپیک رو دانلود کنم.لطفا یه بررسی بکنید.ممنون
ravegoat
17-07-18, 00:15
سلام.راستش من نتونستم هیچکدوم از فایهای این تاپیک رو دانلود کنم.لطفا یه بررسی بکنید.ممنون
سلام دوست عزیز،
به شهر سخت افزار خوش آمدید.
بنده با مدیریت این مشکل رو مطرح می کنم. امیدوارم به زودی مشکل رفع بشه.
البته بنده برخی از فایل های پیوست صفحه ی نهم به بعد رو در آرشیوم موجود دارم. اگر به فایل خاصی نیاز دارید، بفرمایید تا در صورت موجود بودن مجددا بارگذاری کنم.
شاد باشید
آرمین:give_rose:
sm6m10
01-10-19, 15:20
سلام آیا امکانش هست در مورد آشوب در متلب کدی بذارین ممنون
ravegoat
21-10-19, 17:46
سلام آیا امکانش هست در مورد آشوب در متلب کدی بذارین ممنون
سلام دوست عزیز،
به شهر سخت افزار خوش آمدید.
بابت تاخیر پوزش میخوام.
امیدوارم پیوند زیر بهتون کمک کنه:
Only the registered members can see the link