PDA

مشاهده نسخه کامل : سورس ها و آموزش هاي Maple



ravegoat
09-12-10, 07:56
با سلام خدمت دوستان!

اين تاپيك به سورس ها و آموزش هاي نرم افزار Maple اختصاص داره. Maple يه نرم افزار محاسباتي و مدلينگ بر پايه ي سمبل هاست يه چيزي شبيه Mathematica. اين ويژگي باعث ميشه Maple قابليت حل تحليلي داشته باشه. به عبارتي ديگه شما مي تونيد يه معادله رو به طور تحليلي و بر اساس يه سري پارامتر ها حل كنيد و پارامتر ها رو در هر زماني كه مايل باشيد مي تونيد مقدار دهي كنيد اما در نرم افزاري مثل MATLAB چنين چيزي به طور پيش فرض ممكن نيست. براي روشن تر شدن قضيه يه مثال مي زنيم:

انتگرال نا معين تابع cos(ax) بر حسب dx رو در نظر بگيريد. در Maple حاصل اين انتگرال sin(ax)/a خواهد بود اما در MATLAB با اين خطا مواجه مي شويد كه تابع يا متغيري به نام a تعريف نشده.

محيط كد نويسي Maple در سه شاخه قرار ميگره: Text براي بيان توضيحات، Math براي كد نويسي دستورات و Drawing براي نمايش عناصر بصري. به علاوه قابليت نماد گذاري و قرار دادن Component ها رو ميشه در اين محيط ها قرار داد. بنابراين محيط كدنويسي در Maple يكپارچه هست و تحت عنوان Maple Worksheet شناخته ميشه.

اگر چه Maple در محاسبات عددي كند تر از MATLAB عمل مي كنه و به اندازه ي اون گستردگي شاخه نداره اما ابزار هاي تحليل رياضي Maple بسيار قدرتمنده و اون رو به يه نرم افزار قدرتمند در رياضيات كاربردي بدل كرده.

ravegoat
09-12-10, 08:14
سورس مربوط به سري فوريه رو در اينجا پيوست كردم كه تو Maple 14 تست شده.

در اين سورس مي تونيد تابع و دوره ي تناوب رو مشخص كنيد و بعد تعيين كنيد كه سري تا چند جمله محاسبه بشه. متغير flo دقت اعشار محاسبات رو تعيين مي كنه. پس از محاسبه ي سري هم نمودار اون رسم ميشه تا بتونيد مشاهده كنيد كه سري فوريه چقدر به تابع اصلي نزديك شده.



به اميد خدا بعد از امتحانات دو تا آموزش خوب هم از Maple مي ذارم.:great:

ravegoat
24-01-11, 14:06
همون طور كه قول داده بودم دو آموزش از Maple به زبان انگليسي مي ذارم (اين دو تا رو تو IrExpert.ir پيدا كردم.) يكي از اين آموزش ها مقدماتي و اون يكي پيشرفته (چون يكم حجمش بالا بود به دو قطعه تقسيمش كردم). اميدوارم كه از اين آموزش ها در كنار Help نرم افزار بهره ببريد.:great:

ravegoat
12-02-11, 20:53
اين برنامه در Maple براي حل معادله ي تار مرتعش به روش دالامبر هستش. حل دالامبر روش حالت خاصي از روش مشخصه ها هستش كه براي تحليل ارتعاش تاري كه بدون سرعت اوليه رها مي شود كاربرد داره. اين برنامه نحوه ي ارتعاش موج رو هم نشون ميده. به علاوه شما مي تونيد معادله ي تغيير حالت اوليه رو هم به دلخواه تغيير بديد...:wink:

منبع سورسي كه حجمش كمتره توي خود سورس وجود داره.

ravegoat
12-02-11, 21:01
اين برنامه براي تحليل ارتعاش دو بعدي نوشته شده كه به شما ارتعاش ايجاد شده در صفحه براساس IC در نظر گرفته شده، نمايش ميده.

ravegoat
12-02-11, 21:09
اگه هندسه ي يك شكل دو بعدي صورتي شبيه به دايره داشته باشه، معادله ي لاپلاس رو ميشه در مختصات قطبي به كار برد تا حل مسئله ساده تر بشه. اين سورس هم يه معادله ي گرماي لاپلاس رو براي يك دايره به نمايش ميذاره. شما قادر خواهيد بود شرايط مسئله رو بر حسب نياز خودتان تغيير بديد.

رويا 1
29-05-11, 16:02
سلام خسته نباشيد من برنامه حل معادله موج رو به يكي از زبانهاي برنامه نويسي ميخواستم.كسي هست كمكم كنه؟

ravegoat
21-06-11, 17:39
سلام خسته نباشيد من برنامه حل معادله موج رو به يكي از زبانهاي برنامه نويسي ميخواستم.كسي هست كمكم كنه؟

با سلام!

دوست گرامي عضويت شما رو در PCN تبريك مي گم.

لطفا" به پست هاي 4 و 5 مراجعه كنيد. پياده سازي حل معادله موج دالامبر در هر زبان به سادگي امكان پذيره. به علاوه با تعريف روش هاي عددي انتگرال كه در بخش .Net انجمن موجوده، مي تونيد با محاسبه ي سري فوريه معادله موج رو حل كنيد.

موفق باشيد
آرمين

ravegoat
02-08-13, 23:42
در خیلی از مواقع ما در نرم افزار Maple نیاز داریم که خروجی به شکلی منظم در یک فایل Excel ذخیره کنیم و یا آن که دسته ای از اطلاعات را از فایل Excel وارد یک آرایه درون نرم افزار کینم. فایل پیوست شده نحوه ی اجرای این عملیات را در قالب یک مثال به ما نشان می دهد.

منبع (Only the registered members can see the link)

aseman1369_s
01-11-13, 12:13
سلام

من می خوام 2 معادله را باهم حل کنم (البته با شرایط مرزی)
برنامه زیر را نوشته ام ;
eq1 := diff(f(x), x, x, x)+(1/2)*f(x)*(diff(f(x), x, x)) = 0;
eq2 := diff(g(x), x, x)+.4*f(x)*(diff(g(x), x)) = 0;
eq3 := f(0) = 0, (D(f))(0) = 0, (D(f))(10) = 1;
;eq4 := g(0) = 1, g(10) = 0
dd := dsolve({eq1, eq2, eq3, eq4}, {f(x), g(x)}, numeric, output = listprocedure)
g := eval(g(x), dd);
f := eval(f(x), dd)
for k from 5.3 by 0.5e-1 to 5.4 do
print(dd, eval(g(x), x = k))
end do
حالا می خواهم مشتق تابع g را در x=0 بدست بیاورم اما نمی دونم چجوری

لطفا راهنمایی کنید

با تشکر

ravegoat
01-11-13, 20:00
سلام

من می خوام 2 معادله را باهم حل کنم (البته با شرایط مرزی)
برنامه زیر را نوشته ام ;
eq1 := diff(f(x), x, x, x)+(1/2)*f(x)*(diff(f(x), x, x)) = 0;
eq2 := diff(g(x), x, x)+.4*f(x)*(diff(g(x), x)) = 0;
eq3 := f(0) = 0, (D(f))(0) = 0, (D(f))(10) = 1;
;eq4 := g(0) = 1, g(10) = 0
dd := dsolve({eq1, eq2, eq3, eq4}, {f(x), g(x)}, numeric, output = listprocedure)
g := eval(g(x), dd);
f := eval(f(x), dd)
for k from 5.3 by 0.5e-1 to 5.4 do
print(dd, eval(g(x), x = k))
end do
حالا می خواهم مشتق تابع g را در x=0 بدست بیاورم اما نمی دونم چجوری

لطفا راهنمایی کنید

با تشکر
با سلام!
دوست گرامی عضویت شما در شهر سخت افزار تبریک می گم.:party:

اگر شما ضابطه ی تابع g رو داشته باشید، با دستور diff می تونید مشتق اون رو به دست بیارید و با دستور eval مقدار مشتق رو در نقطه ی دلخواه محاسبه کنید:


dg := diff(g, x);
eval(dg, x = 0)


همچنین می تونید با داشتن مقدار g در همسایگی نقطه ی صفر مشتق رو با روش Five-point stencil (Only the registered members can see the link) تخمین بزنید.

موفق باشید
آرمین

aseman1369_s
01-11-13, 22:06
با سلام!
دوست گرامی عضویت شما در شهر سخت افزار تبریک می گم.:party:

اگر شما ضابطه ی تابع g رو داشته باشید، با دستور diff می تونید مشتق اون رو به دست بیارید و با دستور eval مقدار مشتق رو در نقطه ی دلخواه محاسبه کنید:


dg := diff(g, x);
eval(dg, x = 0)


همچنین می تونید با داشتن مقدار g در همسایگی نقطه ی صفر مشتق رو با روش Five-point stencil (Only the registered members can see the link) تخمین بزنید.

موفق باشید
آرمین

این کار را کردم اما مقدار 0 را می دهد !!

آیا در برنامه ای نوشتم مشکلی هست؟

ravegoat
01-11-13, 22:42
این کار را کردم اما مقدار 0 را می دهد !!

آیا در برنامه ای نوشتم مشکلی هست؟
متاسفانه بنده میپل نصب ندارم که سورس شما رو بررسی کنم ولی در یه نگاه کلی که بهش انداختم مشکلی نداشت.

شما اول بررسی کنید که به ازای dd(0) در خروجی آیا مشتق g در x=0 ظاهر میشه.

در صورتی که این ترم وجود داشت سورس زیر رو امتحان کنید؛ شاید جواب داد:


eq1 := diff(f(x), x, x, x)+(1/2)*f(x)*(diff(f(x), x, x)) = 0;
eq2 := diff(g(x), x, x)+.4*f(x)*(diff(g(x), x)) = 0;
eq3 := f(0) = 0, (D(f))(0) = 0, (D(f))(10) = 1;
eq4 := g(0) = 1, g(10) = 0;
dd := dsolve({eq1, eq2, eq3, eq4}, {f(x), g(x)}, numeric, output = operator)
dg := eval(D(g)(x), dd);
dg(0)


در اولین فرصت هم اگر میپل پیدا کردم سورس ها رو مجددا بررسی می کنم.

موفق باشید
آرمین

aseman1369_s
01-11-13, 23:01
متاسفانه بنده میپل نصب ندارم که سورس شما رو بررسی کنم ولی در یه نگاه کلی که بهش انداختم مشکلی نداشت.

شما اول بررسی کنید که به ازای dd(0) در خروجی آیا مشتق g در x=0 ظاهر میشه.

در صورتی که این ترم وجود داشت سورس زیر رو امتحان کنید؛ شاید جواب داد:


eq1 := diff(f(x), x, x, x)+(1/2)*f(x)*(diff(f(x), x, x)) = 0;
eq2 := diff(g(x), x, x)+.4*f(x)*(diff(g(x), x)) = 0;
eq3 := f(0) = 0, (D(f))(0) = 0, (D(f))(10) = 1;
eq4 := g(0) = 1, g(10) = 0;
dd := dsolve({eq1, eq2, eq3, eq4}, {f(x), g(x)}, numeric, output = operator)
dg := eval(D(g)(x), dd);
dg(0)


در اولین فرصت هم اگر میپل پیدا کردم سورس ها رو مجددا بررسی می کنم.

موفق باشید
آرمین

بالاخره تونستم خروجی بگیرم
اینم دستورش
dg := eval(diff(g(x), x), dd)
print(dd, eval(dg(x), x = 0))
با تشکر از پیگیری شما

sara_731
11-12-13, 17:08
سلام
اگه میشه این تمرین و برای من حل کنید
خیلی ممنونم
من هیچ آشنایی با با نرم افزار mapleندارم که از تاپیک های قبلی استفاده کنم و همچنین آخرین محلت تحویل شنبه است که در این مدت کوتاه این نرم افزار پیچیده یادگیریش خیلی برام سخته.

متشکرم

ravegoat
15-12-13, 07:41
سلام
اگه میشه این تمرین و برای من حل کنید
خیلی ممنونم
من هیچ آشنایی با با نرم افزار mapleندارم که از تاپیک های قبلی استفاده کنم و همچنین آخرین محلت تحویل شنبه است که در این مدت کوتاه این نرم افزار پیچیده یادگیریش خیلی برام سخته.

متشکرم
با سلام!

دوست گرامی عضویت شما رو در شهر سخت افزار تبریک می گم.

نحوه ی حل معادله ی موج یک بعدی و دو بعدی در پیوند های زیر تشریح شده:
pde.html (Only the registered members can see the link)
Math 450/Maple Demos/preview/2d-Wave.html (Only the registered members can see the link)

به علاوه جهت کسب اطلاعات بیش تر پیرامون حل معادلات با مشتقات جزئی در Maple می تونید آموزش زیر رو مطالعه بفرمایید:
Only the registered members can see the link

موفق باشید
آرمین

sara_731
15-12-13, 21:25
وقتی که من میخوام اینو حل کنم فقط یه نمودار خالی میکشه! اشکال از کجاست؟

u:=sum((c[n]*cos(lambda[n]*t)+d[n]*sin(lambda[n]*t))*sin(n*Pi*x/L),n=1..40);
(1.117372821 + 1.220847109 sinx cosx) sin(x)

+ (0.03002180081 + 2.540100125 sinx cosx) sin(3. x)

+ (-0.2040095464 + 1.523996800 sinx cosx) sin(5. x)

+ (-0.1703333837 - 0.8932621544 sinx cosx) sin(7. x)

+ (-0.02982150529 - 2.489260003 sinx cosx) sin(9. x)

+ (0.08202785276 - 1.796643686 sinx cosx) sin(11. x)

+ (0.09564846028 + 0.5477985513 sinx cosx) sin(13. x)

+ (0.02942331960 + 2.388597326 sinx cosx) sin(15. x)

+ (-0.04508855063 + 2.033330735 sinx cosx) sin(17. x)

+ (-0.06682310642 - 0.1913707567 sinx cosx) sin(19. x)

+ (-0.02883202008 - 2.240126858 sinx cosx) sin(21. x)

+ (0.02675490134 - 2.229320656 sinx cosx) sin(23. x)

+ (0.05081744906 - 0.1688873247 sinx cosx) sin(25. x)

+ (0.02805468571 + 2.046820234 sinx cosx) sin(27. x)

+ (-0.01558288316 + 2.380690709 sinx cosx) sin(29. x)

+ (-0.04018727876 + 0.5257651252 sinx cosx) sin(31. x)

+ (-0.02710059726 - 1.812546490 sinx cosx) sin(33. x)

+ (0.007982846823 - 2.484411222 sinx cosx) sin(35. x)

+ (0.03233081750 - 0.8721197336 sinx cosx) sin(37. x)

+ (0.02598110487 + 1.541994616 sinx cosx) sin(39. x)

plot(u,x=0..L);
Warning, expecting only range variable x in expression (1.117372821+1.220847109*sinx*cosx)*sin(x)+(.30021 80081e-1+2.540100125*sinx*cosx)*sin(3.*x)+(-.2040095464+1.523996800*sinx*cosx)*sin(5.*x)+(-.1703333837-.8932621544*sinx*cosx)*sin(7.*x)+(-.2982150529e-1-2.489260003*sinx*cosx)*sin(9.*x)+(.8202785276e-1-1.796643686*sinx*cosx)*sin(11.*x)+(.9564846028e-1+.5477985513*sinx*cosx)*sin(13.*x)+(.2942331960e-1+2.388597326*sinx*cosx)*sin(15.*x)+(-.4508855063e-1+2.033330735*sinx*cosx)*sin(17.*x)+(-.6682310642e-1-.1913707567*sinx*cosx)*sin(19.*x)+(-.2883202008e-1-2.240126858*sinx*cosx)*sin(21.*x)+(.2675490134e-1-2.229320656*sinx*cosx)*sin(23.*x)+(.5081744906e-1-.1688873247*sinx*cosx)*sin(25.*x)+(.2805468571e-1+2.046820234*sinx*cosx)*sin(27.*x)+(-.1558288316e-1+2.380690709*sinx*cosx)*sin(29.*x)+(-.4018727876e-1+.5257651252*sinx*cosx)*sin(31.*x)+(-.2710059726e-1-1.812546490*sinx*cosx)*sin(33.*x)+(.7982846823e-2-2.484411222*sinx*cosx)*sin(35.*x)+(.3233081750e-1-.8721197336*sinx*cosx)*sin(37.*x)+(.2598110487e-1+1.541994616*sinx*cosx)*sin(39.*x) to be plotted but found names [cosx, sinx]

ravegoat
17-12-13, 08:07
وقتی که من میخوام اینو حل کنم فقط یه نمودار خالی میکشه! اشکال از کجاست؟

u:=sum((c[n]*cos(lambda[n]*t)+d[n]*sin(lambda[n]*t))*sin(n*Pi*x/L),n=1..40);
(1.117372821 + 1.220847109 sinx cosx) sin(x)

+ (0.03002180081 + 2.540100125 sinx cosx) sin(3. x)

+ (-0.2040095464 + 1.523996800 sinx cosx) sin(5. x)

+ (-0.1703333837 - 0.8932621544 sinx cosx) sin(7. x)

+ (-0.02982150529 - 2.489260003 sinx cosx) sin(9. x)

+ (0.08202785276 - 1.796643686 sinx cosx) sin(11. x)

+ (0.09564846028 + 0.5477985513 sinx cosx) sin(13. x)

+ (0.02942331960 + 2.388597326 sinx cosx) sin(15. x)

+ (-0.04508855063 + 2.033330735 sinx cosx) sin(17. x)

+ (-0.06682310642 - 0.1913707567 sinx cosx) sin(19. x)

+ (-0.02883202008 - 2.240126858 sinx cosx) sin(21. x)

+ (0.02675490134 - 2.229320656 sinx cosx) sin(23. x)

+ (0.05081744906 - 0.1688873247 sinx cosx) sin(25. x)

+ (0.02805468571 + 2.046820234 sinx cosx) sin(27. x)

+ (-0.01558288316 + 2.380690709 sinx cosx) sin(29. x)

+ (-0.04018727876 + 0.5257651252 sinx cosx) sin(31. x)

+ (-0.02710059726 - 1.812546490 sinx cosx) sin(33. x)

+ (0.007982846823 - 2.484411222 sinx cosx) sin(35. x)

+ (0.03233081750 - 0.8721197336 sinx cosx) sin(37. x)

+ (0.02598110487 + 1.541994616 sinx cosx) sin(39. x)

plot(u,x=0..L);
Warning, expecting only range variable x in expression (1.117372821+1.220847109*sinx*cosx)*sin(x)+(.30021 80081e-1+2.540100125*sinx*cosx)*sin(3.*x)+(-.2040095464+1.523996800*sinx*cosx)*sin(5.*x)+(-.1703333837-.8932621544*sinx*cosx)*sin(7.*x)+(-.2982150529e-1-2.489260003*sinx*cosx)*sin(9.*x)+(.8202785276e-1-1.796643686*sinx*cosx)*sin(11.*x)+(.9564846028e-1+.5477985513*sinx*cosx)*sin(13.*x)+(.2942331960e-1+2.388597326*sinx*cosx)*sin(15.*x)+(-.4508855063e-1+2.033330735*sinx*cosx)*sin(17.*x)+(-.6682310642e-1-.1913707567*sinx*cosx)*sin(19.*x)+(-.2883202008e-1-2.240126858*sinx*cosx)*sin(21.*x)+(.2675490134e-1-2.229320656*sinx*cosx)*sin(23.*x)+(.5081744906e-1-.1688873247*sinx*cosx)*sin(25.*x)+(.2805468571e-1+2.046820234*sinx*cosx)*sin(27.*x)+(-.1558288316e-1+2.380690709*sinx*cosx)*sin(29.*x)+(-.4018727876e-1+.5257651252*sinx*cosx)*sin(31.*x)+(-.2710059726e-1-1.812546490*sinx*cosx)*sin(33.*x)+(.7982846823e-2-2.484411222*sinx*cosx)*sin(35.*x)+(.3233081750e-1-.8721197336*sinx*cosx)*sin(37.*x)+(.2598110487e-1+1.541994616*sinx*cosx)*sin(39.*x) to be plotted but found names [cosx, sinx]

اون طور که به نظر میاد معادله ی شما هم با مکان x تغییر می کنه و هم زمان t ! آیا شما از قبل t رو مقدار دهی کردید؟ چون شکل جواب صحیح به نظر نمی رسه. اگر شما t رو از قبل مقدار دهی کنید باید عبارت زیر:


(c[n]*cos(lambda[n]*t)+d[n]*sin(lambda[n]*t)


یک عدد بشه در حالی که حاصل به شکل تابعی از x شده:


p+q*sinx*cosx


خطای شما هم به این دلیل که sinx و cosx عبارت های معینی نیستند و شکل صحیح اون ها باید sin(x) و cos(x) باشه.

sara_731
17-12-13, 20:48
سلام خیلی از شما ممنونم که به سوالات من جواب می دین.
وقتی که من g:=sin2x به شکل g:=sin(2x) تعریف می کنم ارور میده نمی دونم چرا ؟؟؟؟
کل برنامه این شکلیه که می خوام تو زمانهای 0و1.1.5و2و2.5و3 شکل تارمرتعش به ازای x از 0 تا پی بکشم که نمی شه نمودار خالی می کشه
f:=k((Pi/2)^4-(x-Pi/2)^4);
/ 4\
|1 4 / 1 \ |
k|-- Pi - |x - - Pi| |
\16 \ 2 / /
g:=sin2x;
sin2x
lambda[n] := c*n*Pi/L;
c n Pi
------
L
c[n]:=(2/L)*int(f*sin(n*Pi*x/L),x=0..L);
2 f (cos(n Pi) - 1)
- -------------------
n Pi
d[n]:=(2/L*lambda[n])*int(g*sin(n*Pi*x/L),x=0..L);
2 c sin2x (cos(n Pi) - 1)
- -------------------------
L
L:=Pi;
Pi
k:=1;
1
c:=1;
1
g:=sin2x;
sin2x
lambda[n] := c*n*Pi/L;
n
c[n]:=(2/L)*int(f*sin(n*Pi*x/L),x=0..L);
2 f (cos(n Pi) - 1)
- -------------------
n Pi
d[n]:=(2/L*lambda[n])*int(g*sin(n*Pi*x/L),x=0..L);
2 sin2x (cos(n Pi) - 1)
- -----------------------
Pi
u:=sum((c[n]*cos(lambda[n]*t)+d[n]*sin(lambda[n]*t))*sin(n*Pi*x/L),n=1..40);
/4 f cos(t) 4 sin2x sin(t)\
|---------- + --------------| sin(x)
\ Pi Pi /

/4 f cos(3 t) 4 sin2x sin(3 t)\
+ |------------ + ----------------| sin(3 x)
\ 3 Pi Pi /

/4 f cos(5 t) 4 sin2x sin(5 t)\
+ |------------ + ----------------| sin(5 x)
\ 5 Pi Pi /

/4 f cos(7 t) 4 sin2x sin(7 t)\
+ |------------ + ----------------| sin(7 x)
\ 7 Pi Pi /

/4 f cos(9 t) 4 sin2x sin(9 t)\
+ |------------ + ----------------| sin(9 x)
\ 9 Pi Pi /

/4 f cos(11 t) 4 sin2x sin(11 t)\
+ |------------- + -----------------| sin(11 x)
\ 11 Pi Pi /

/4 f cos(13 t) 4 sin2x sin(13 t)\
+ |------------- + -----------------| sin(13 x)
\ 13 Pi Pi /

/4 f cos(15 t) 4 sin2x sin(15 t)\
+ |------------- + -----------------| sin(15 x)
\ 15 Pi Pi /

/4 f cos(17 t) 4 sin2x sin(17 t)\
+ |------------- + -----------------| sin(17 x)
\ 17 Pi Pi /

/4 f cos(19 t) 4 sin2x sin(19 t)\
+ |------------- + -----------------| sin(19 x)
\ 19 Pi Pi /

/4 f cos(21 t) 4 sin2x sin(21 t)\
+ |------------- + -----------------| sin(21 x)
\ 21 Pi Pi /

/4 f cos(23 t) 4 sin2x sin(23 t)\
+ |------------- + -----------------| sin(23 x)
\ 23 Pi Pi /

/4 f cos(25 t) 4 sin2x sin(25 t)\
+ |------------- + -----------------| sin(25 x)
\ 25 Pi Pi /

/4 f cos(27 t) 4 sin2x sin(27 t)\
+ |------------- + -----------------| sin(27 x)
\ 27 Pi Pi /

/4 f cos(29 t) 4 sin2x sin(29 t)\
+ |------------- + -----------------| sin(29 x)
\ 29 Pi Pi /

/4 f cos(31 t) 4 sin2x sin(31 t)\
+ |------------- + -----------------| sin(31 x)
\ 31 Pi Pi /

/4 f cos(33 t) 4 sin2x sin(33 t)\
+ |------------- + -----------------| sin(33 x)
\ 33 Pi Pi /

/4 f cos(35 t) 4 sin2x sin(35 t)\
+ |------------- + -----------------| sin(35 x)
\ 35 Pi Pi /

/4 f cos(37 t) 4 sin2x sin(37 t)\
+ |------------- + -----------------| sin(37 x)
\ 37 Pi Pi /

/4 f cos(39 t) 4 sin2x sin(39 t)\
+ |------------- + -----------------| sin(39 x)
\ 39 Pi Pi /
t:=0;
0
u:=sum((c[n]*cos(lambda[n]*t)+d[n]*sin(lambda[n]*t))*sin(n*Pi*x/L),n=1..40);
4 f sin(x) 4 f sin(3 x) 4 f sin(5 x) 4 f sin(7 x)
---------- + ------------ + ------------ + ------------
Pi 3 Pi 5 Pi 7 Pi

4 f sin(9 x) 4 f sin(11 x) 4 f sin(13 x) 4 f sin(15 x)
+ ------------ + ------------- + ------------- + -------------
9 Pi 11 Pi 13 Pi 15 Pi

4 f sin(17 x) 4 f sin(19 x) 4 f sin(21 x)
+ ------------- + ------------- + -------------
17 Pi 19 Pi 21 Pi

4 f sin(23 x) 4 f sin(25 x) 4 f sin(27 x)
+ ------------- + ------------- + -------------
23 Pi 25 Pi 27 Pi

4 f sin(29 x) 4 f sin(31 x) 4 f sin(33 x)
+ ------------- + ------------- + -------------
29 Pi 31 Pi 33 Pi

4 f sin(35 x) 4 f sin(37 x) 4 f sin(39 x)
+ ------------- + ------------- + -------------
35 Pi 37 Pi 39 Pi
with(plots):
plot(u,x=0..L);
Warning, expecting only range variable x in expression 4*f/Pi*sin(x)+4/3*f/Pi*sin(3*x)+4/5*f/Pi*sin(5*x)+4/7*f/Pi*sin(7*x)+4/9*f/Pi*sin(9*x)+4/11*f/Pi*sin(11*x)+4/13*f/Pi*sin(13*x)+4/15*f/Pi*sin(15*x)+4/17*f/Pi*sin(17*x)+4/19*f/Pi*sin(19*x)+4/21*f/Pi*sin(21*x)+4/23*f/Pi*sin(23*x)+4/25*f/Pi*sin(25*x)+4/27*f/Pi*sin(27*x)+4/29*f/Pi*sin(29*x)+4/31*f/Pi*sin(31*x)+4/33*f/Pi*sin(33*x)+4/35*f/Pi*sin(35*x)+4/37*f/Pi*sin(37*x)+4/39*f/Pi*sin(39*x) to be plotted but found name f

ravegoat
18-12-13, 09:33
سلام خیلی از شما ممنونم که به سوالات من جواب می دین.
وقتی که من g:=sin2x به شکل g:=sin(2x) تعریف می کنم ارور میده نمی دونم چرا ؟؟؟؟
....

خواهش می کنم...

مشکل به وجود اومده به خاطر اینه که sin(2x) به شکل sin(2*x) تعریف نشده!

کد به شکل زیر اصلاح شد:


restart;
f:=k((Pi/2)^4-(x-Pi/2)^4);
g:=sin(2*x);
lambda[n] := c*n*Pi/L;
c[n]:=(2/L)*int(f*sin(n*Pi*x/L),x=0..L);
d[n]:=(2/L*lambda[n])*int(g*sin(n*Pi*x/L),x=0..L);
L:=Pi;
k:=1;
c:=1;
lambda[n] := c*n*Pi/L;
c[n]:=(2/L)*int(f*sin(n*Pi*x/L),x=0..L);
d[n]:=(2/L*lambda[n])*int(g*sin(n*Pi*x/L),x=0..L);
u:=sum((c[n]*cos(lambda[n]*t)+d[n]*sin(lambda[n]*t))*sin(n*Pi*x/L),n=1..40);
with(plots):
animate(plot,[u,x=0..L],t=0..3,frames=7);



موفق باشید
آرمین

sara_731
19-12-13, 20:06
سلام استاد گرامی
خیلی گرفتار بودم با راهنمایی هایت در مورد مشکلاتم در پروژه ام مرا آزاد کردی و با بزرگواری بی دریغتان مرا یاری نموده تا نتیجه ی کامل بگیرم ممنونم.
با تبریک پیشاپیش شب یلدا, با تشکر شاگرد شما سارا

majidmahjoor
27-12-13, 22:45
سلام
کد های کیپل VIM,HPM,PM,ADM,EBM رو می خواستم.سیمیتونه کمکم بکنه
variation iteration method
homotopy perturbation method
perturbation mthod
Adomian decomposition method
energy balance method

ravegoat
28-12-13, 18:03
سلام
کد های کیپل VIM,HPM,PM,ADM,EBM رو می خواستم.سیمیتونه کمکم بکنه
variation iteration method
homotopy perturbation method
perturbation mthod
Adomian decomposition method
energy balance method
با سلام!

VIM:
Only the registered members can see the link

HPM:
Only the registered members can see the link

PM:
Only the registered members can see the link

ADM:
Only the registered members can see the link

موازنه انرژی هم یافت نشد.

موفق باشید
آرمین

ata1
24-02-14, 15:51
سلام خسته نباشید
تو میپل میشه روش حساب تغییرات رو انجام داد؟
در واقع برای حل معادله لاگرانژ از اصل همیلتون، می تونیم تو میپل وریشن (variation) بگیریم ؟
اگر تو این زمینه فایل آموزشی یا مثالی هم امکان داشته باشه آپلود کنین ممنون میشم
با تشکر

ravegoat
24-02-14, 17:50
سلام خسته نباشید
تو میپل میشه روش حساب تغییرات رو انجام داد؟
در واقع برای حل معادله لاگرانژ از اصل همیلتون، می تونیم تو میپل وریشن (variation) بگیریم ؟
اگر تو این زمینه فایل آموزشی یا مثالی هم امکان داشته باشه آپلود کنین ممنون میشم
با تشکر
با سلام!

دوست عزیز عضویت شما رو در شهر سخت افزار تبریک می گم.

بله، بسته ی VariationalCalculus نرم افزار Maple قادر به چنین کار هستش که برای اطلاعات بیش تر می تونید به پیوند زیر مراجعه بفرمایید:
Overview of the VariationalCalculus Package - Maple Help (Only the registered members can see the link)

اگر به مثال هایی هم در این زمینه علاقه مند هستید، مشاهده ی پیوند های زیر توصیه میشه:
Only the registered members can see the link
Variational Methods—Wolfram Mathematica 9 Documentation (Only the registered members can see the link)

موفق باشید
آرمین

mamadrst123
13-08-14, 22:15
سلام
اگه میشه این تمرین و برای من حل کنید
خیلی ممنون می شم.
من هیچ آشنایی با با نرم افزار mapleندارم از تاپیک ها هم استفاده کردم خوندم ولی چیزی نفهمیدم.

m(d^2x/dt^2)=-kx-c(dx/dt)
11135

ravegoat
14-08-14, 19:44
سلام
اگه میشه این تمرین و برای من حل کنید
خیلی ممنون می شم.
من هیچ آشنایی با با نرم افزار mapleندارم از تاپیک ها هم استفاده کردم خوندم ولی چیزی نفهمیدم.

m(d^2x/dt^2)=-kx-c(dx/dt)
11135
با سلام!

دوست گرامی عضویت شما رو در شهر سخت افزار تبریک می گم.

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

موفق باشید
آرمین

mamadrst123
21-08-14, 12:28
با سلام!

دوست گرامی عضویت شما رو در شهر سخت افزار تبریک می گم.

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

موفق باشید
آرمین

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

ravegoat
21-08-14, 16:08
سلام آقا آرمین ممنون از تبریکتون..
و ممنون از پاسخی که دادید.
امکانش هست لطف کنید معادله حل کنید؟
فایل پیوست عکس از روی معادله گرفته بودم که باز نمیشه..
خواهش می کنم...

همان طور که عرض کردم با دستور dsolve می تونید که البته حلش براتون پ.خ. شد.

ولی برای رسم مواردی که در عکس بهش اشاره شده می تونید از دستور plot (Only the registered members can see the link) استفاده فرمایید.

موفق باشید
آرمین

- - - Updated - - -


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

همان طور که عرض کردم با دستور dsolve می تونید که البته حلش براتون پ.خ. شد.

ولی برای رسم مواردی که در عکس بهش اشاره شده می تونید از دستور plot (Only the registered members can see the link) استفاده فرمایید.

موفق باشید
آرمین

mahmass
08-12-14, 14:10
با سلام
برنامه اسمبل کردن ماتریس ها در میپل رو نیاز دارم ممنون میشم کمکم کنید
برای حل مسائل المانمحدود نیاز دارم

ravegoat
08-12-14, 20:02
با سلام
برنامه اسمبل کردن ماتریس ها در میپل رو نیاز دارم ممنون میشم کمکم کنید
برای حل مسائل المانمحدود نیاز دارم
با سلام!

دوست گرامی به شهر سخت افزار خوش آمدید.

لطفا" به پیوند زیر مراجعه کنید:
Only the registered members can see the link
اینم کدش:
Only the registered members can see the link

برای اطلاعات بیش تر:
Finite Elements Using Maple: A Symbolic Programming Approach - Artur Portela, A. Charafi - Google Books (Only the registered members can see the link=onepage&q=fem%20matrix%20assembling%20in%20maple&f=false)

موفق باشید
آرمین

ho3ein6900
18-01-15, 23:25
سلام.من میخوام این مقاله رو با میپل بنویسم.کسی میتونه کمکم کنه؟
پرشین گیگ | دانلود Sciencedirect paper received - libdl.ir.pdf (Only the registered members can see the link)

ravegoat
20-01-15, 08:24
سلام.من میخوام این مقاله رو با میپل بنویسم.کسی میتونه کمکم کنه؟
پرشین گیگ | دانلود Sciencedirect paper received - libdl.ir.pdf (Only the registered members can see the link)
با سلام!

دوست گرامی به شهر سخت افزار خوش آمدید.

متاسفانه امکان دانلود PDF مورد نظر وجود ندارد. لطفا" PDF رو در همین جا آپلود کنید و بفرمایید در مورد پیاده سازی کدوم قسمتش مشکل دارید.

موفق باشید
آرمین

ho3ein6900
20-01-15, 22:37
Only the registered members can see the linkمیشه ی لطفی کنین بگین چرا ئستور اخری تو عکس ارو داده؟Only the registered members can see the link


با سلام!

دوست گرامی به شهر سخت افزار خوش آمدید.

متاسفانه امکان دانلود PDF مورد نظر وجود ندارد. لطفا" PDF رو در همین جا آپلود کنید و بفرمایید در مورد پیاده سازی کدوم قسمتش مشکل دارید.

موفق باشید
آرمین
مهندس ی نگاه ب این عکس بنداز ببین مشکل از کجاس؟

ravegoat
09-04-15, 08:18
مهندس ی نگاه ب این عکس بنداز ببین مشکل از کجاس؟
دوست عزیز لطفا" دستور coeff(equ2,p) رو به تنهایی اجرا کنید ببینید نتیجه چی میشه؟

ali10 sa91
01-11-15, 22:08
سلام ببخشید کسی کدهای درون یابی خطی از جمله روش لاگرانژ و نیوتن رو داره

ravegoat
02-11-15, 08:15
سلام ببخشید کسی کدهای درون یابی خطی از جمله روش لاگرانژ و نیوتن رو داره
با سلام!

دوست گرامی به شهر سخت افزار خوش آمدید...

می تونید از دستور PolynomialInterpolation (Only the registered members can see the link omialInterpolation) استفاده بفرمایید.

برای اطلاعات بیش تر:
Only the registered members can see the link

موفق باشید
آرمین

naser38
15-02-16, 20:48
سلام دوست عزیز بنده نیاز شدید به حل معادله موج دو بعدی با میپل دارم ولی هر کاری میکنم فایلتون دانلود نمیشه اگه امکانش هست واسم ایمیل کنید از همکاری شما بسیار سپاس گذارم
ایمیل:naser_shahbazi38yahoo.com

naser38
18-02-16, 19:49
سلام دوستان کسی می تونه معادله موج دو بعدی و با شرایط اولیه و مرزی دلخواه با maple حل کنه:
diff(u(x, y, t), t, t) = c^2*(diff(u(x, y, t), y, y)+diff(u(x, y, t), x, x))

ravegoat
27-02-16, 09:25
سلام دوست عزیز بنده نیاز شدید به حل معادله موج دو بعدی با میپل دارم ولی هر کاری میکنم فایلتون دانلود نمیشه اگه امکانش هست واسم ایمیل کنید از همکاری شما بسیار سپاس گذارم
ایمیل:naser_shahbazi38@yahoo.com (Only the registered members can see the link)
با سلام!

دوست گرامی به شهر سخت افزار خوش آمدید.

این سورس رو می تونید از پست 5 دریافت کنید:
Only the registered members can see the link

پی نوشت:
دانلود پیوست بررسی شد که خدا رو شکر مشکلی نداشت.

شاد و موفق باشید
آرمین

eenvid
29-11-16, 14:19
اقا اين سورس سري فوريه كه گذاشتيد واسه maple 16 هم كار ميكنه؟؟ دو روزِ اينترنت رو ميگردم، هرچي سورس ميگيرم تو maple 16 كار نميكنه

eenvid
29-11-16, 18:49
آقا نمیدونم تا حالا چند نفر بهت گفتن عاشقتن ولی من یکی که عاشقتم.:love::love::love:
اون موقع سر کار بودم،اومدم خونه سورس رو چک کردم کار کرد.
فقط من این دوتا دستور رو هم به دستورات شما اضافه کردم که k رو عدد ثابت در نظر بگیره و جواب معادله ساده تر بشه
k:='k'; assume(k,integer)
به جای دو نقطه ها هم باید سمیکالن باشه.
اگه بخوایم نتیجه عبارت به صورت کسری باقی بمونه باید evalf[flo] رو حذف کنیم.

میدونم اینارو شما میدونید،اینارو گفتم برای کسایی که مثل من خواستن نتیجه عبارت به صورت کسری باقی بمونه

Seresty
24-01-17, 13:19
با سلام
نحوه کد نویسی برای حل معادله دیفرانسیل با سه متغیر به روش hpm ,vim, numerid رو در میپل میخواستم...مشتقات بالا...
با تشکر

mehrshadtiger
31-08-18, 11:22
سلام دوستان میتونید کمکم کنید تا معادله هایی که دراونها hypergeom و bessel دارنو چجوری باید حل کنم و ضرایب اونو بدست بیارم

padideh_hs
06-04-19, 11:25
سلام روز بخیر
ممنون از سایت بسیار خوبتون
این فایل متاسفانه بنده نمیتونم دانلود کنم لطف میکنید مشکل را بر طرف نمایید
با تشکر