Printable View
اگه می شه این سوالات رو سریع جواب بدید خیلی لازم دارم:
اول اینکه چرا z که یه مجهوله تو طرف دوم تساوی وجود داره این طوری که نمی شه z رو بدست آورد
دوم اینکه این چند خط پشت سرهم هستن یعنی وقتی می خواد مثلا x رو بنویسه چند خط می نویسه اینا پشت سر هم هستن؟چون مثل اینکه خط ها همشون یه چیز نوشته
سوم اینکه می شه این کارا رو یه جوری به vb منتقل کنیم.چون من می خوام این دستگاه رو به صورت یه تابع تو visual basic بنویسم.
با سلام!
در دستور solve باید اینگونه تغییر داد :
با سپاس !کد:[x y z]=solve(equ,'x','y','z')
در جواب خاطرات باید بگم میشه ولی خیلی سخت. همونطور که میدونی متلب یه نرم افزار کد بازه پس هم میتونی توش تغییرات ایجاد کنی و هم ازش کپی بگیری اینو تو دایرکتوری متلب سرچ کن 'solve' اون وقت یه فایل میاد که کد حل معادله تو اونه ولی تو اول باید معدله رو به صوت ماتریس در بیاری و بعد از اون تو vb ایتفاده کنی در ضمن برنامه نویسی تو متلب خیلی خیلی شبه جاواس اگه بازهم سئوالی بود بپرس.
یا حق
برنامه ای که نوشته شده خیلی خیلی وقت گیره .چون از 1 تا p-1 رو دونه به دونه چک میکنه
میتونید به جای p-1 از جذر p استفاده کنید
:20: موافقمنقل قول:
اگه يه عدد خيلي بزرگ به برنامه تون بدين كه ....اصلا پاسكال نمي تونه اعداد خيلي بزرگ رو قبول كنه چه برسه شما باهاش برنامه بنويسي!يه عددي با 100 رقم بده ببين چطور ميشه؟---- با جذر گرفتن نميشه اول بودن عددي رو تشخيص داد.فقط: اگر جذر يك عددي گويا بود اول نيست.----اين سايت ها اول يه بار گذاشتن كامپيوتر هاي بزرگ اعداد اول رو محاسبه كردن ؛بعد اون اعداد رو ذخيره كردن با دادن عدد به سرعت تعيين ميشه اوله يا نه!موفق باشيد./
الان كه ديدم اين سايت با اسكريپت اين كار رو مي كنه
بروبچ سلام
pp8khat خوشحالم که میبینم هنوز داری داری پست میدی منم بخاطر همین دوباره اومدم دارم پست میدم.
بچه ها این کد هم واسه متلبه و کاربرد بیشتری هم داره:
اگه کد های بیشتری در باره متلب میخواین یه سر به وبلاگ من بزنید.
این کد 0و1 تولید میکنه وکد:function isp = isprime(X)
if isempty(X), isp = false(size(X)); return, end
if any(X(:) < 0) || any(floor(X(:)) ~= X(:))
error('MATLAB:isprime:InputNotPosInt',...
'All entries of X must be positive integers.');
end
isp = false(size(X));
n = max(X(:));
if n > 2^32
error('MATLAB:isprime:InputOutOfRange',...
'The maximum value of X allowed is 2^32.');
end
p = primes(ceil(sqrt(n)));
for k = 1:numel(isp)
isp(k) = all(rem(X(k), p(p<X(k))));
end
% p(p<1) would give an empty matrix and all([]) returns true.
% we need to correct isp for this case.
isp(X==1 | X==0)=0;
کد زیر هم فاکتور های عدد رو نشون میده:
یا حقکد:if numel(n)~=1, error('MATLAB:factor:NonScalarInput','N must be a scalar.'); end
if (n < 0) || (floor(n) ~= n)
error('MATLAB:factor:InputNotPosInt', 'N must be a positive integer.');
end
if n > 2^32
error('MATLAB:factor:InputOutOfRange','The maximum value of n allowed is 2^32.');
end
if n < 4
f=n;
return
else
f = [];
end
p = primes(sqrt(n));
while n>1,
d = find(rem(n,p)==0);
if isempty(d)
f = [f n];
break;
end
p = p(d);
f = [f p];
n = n/prod(p);
end
f = sort(f);
من نگفتم میشه با جذر گرفتن اول بودن عددي رو تشخيص داد من گفتم میشه حلقه ی برنامه رو تا جذر p کوچیکتر کرد.نقل قول:
با جذر گرفتن نميشه اول بودن عددي رو تشخيص داد.
سلام
من مي خوام اين 6 معادله رو تو متلب حل كنم
>> equ='(-((cos(q1)*cos(q2)*cos(q3)-cos(q1)*sin(q2)*sin(q3))*cos(q4)+sin(q1)*sin(q4))* sin(q5)+(cos(q1)*cos(q2)*sin(q3)+cos(q1)*sin(q2)*c os(q3))*cos(q5))*d6+(-cos(q1)*cos(q2)*sin(q3)-cos(q1)*sin(q2)*cos(q3))*d4+cos(q1)*cos(q2)*a3*cos (q3)-cos(q1)*sin(q2)*a3*sin(q3)+cos(q1)*a2*cos(q2)+a1*c os(q1)-f1,(-((sin(q1)*cos(q2)*cos(q3)-sin(q1)*sin(q2)*sin(q3))*cos(q4)-cos(q1)*sin(q4))*sin(q5)+(sin(q1)*cos(q2)*sin(q3)+ sin(q1)*sin(q2)*cos(q3))*cos(q5))*d6+(-sin(q1)*cos(q2)*sin(q3)-sin(q1)*sin(q2)*cos(q3))*d4+sin(q1)*cos(q2)*a3*cos (q3)-sin(q1)*sin(q2)*a3*sin(q3)+sin(q1)*a2*cos(q2)+a1*s in(q1)-f2,(-(-sin(q2)*cos(q3)-cos(q2)*sin(q3))*cos(q4)*sin(q5)+(-sin(q2)*sin(q3)+cos(q2)*cos(q3))*cos(q5))*d6+(sin( q2)*sin(q3)-cos(q2)*cos(q3))*d4-sin(q2)*a3*cos(q3)-cos(q2)*a3*sin(q3)-a2*sin(q2)+d1-f3,exp(q6/pi)*(-((cos(q1)*cos(q2)*cos(q3)-cos(q1)*sin(q2)*sin(q3))*cos(q4)+sin(q1)*sin(q4))* sin(q5)+(cos(q1)*cos(q2)*sin(q3)+cos(q1)*sin(q2)*c os(q3))*cos(q5))-f4,exp(q6/pi)*(-((sin(q1)*cos(q2)*cos(q3)-sin(q1)*sin(q2)*sin(q3))*cos(q4)-cos(q1)*sin(q4))*sin(q5)+(sin(q1)*cos(q2)*sin(q3)+ sin(q1)*sin(q2)*cos(q3))*cos(q5))-f5,-(-sin(q2)*cos(q3)-cos(q2)*sin(q3))*cos(q4)*sin(q5)+(-sin(q2)*sin(q3)+cos(q2)*cos(q3))*cos(q5)-f6'
equ =
(-((cos(q1)*cos(q2)*cos(q3)-cos(q1)*sin(q2)*sin(q3))*cos(q4)+sin(q1)*sin(q4))* sin(q5)+(cos(q1)*cos(q2)*sin(q3)+cos(q1)*sin(q2)*c os(q3))*cos(q5))*d6+(-cos(q1)*cos(q2)*sin(q3)-cos(q1)*sin(q2)*cos(q3))*d4+cos(q1)*cos(q2)*a3*cos (q3)-cos(q1)*sin(q2)*a3*sin(q3)+cos(q1)*a2*cos(q2)+a1*c os(q1)-f1,(-((sin(q1)*cos(q2)*cos(q3)-sin(q1)*sin(q2)*sin(q3))*cos(q4)-cos(q1)*sin(q4))*sin(q5)+(sin(q1)*cos(q2)*sin(q3)+ sin(q1)*sin(q2)*cos(q3))*cos(q5))*d6+(-sin(q1)*cos(q2)*sin(q3)-sin(q1)*sin(q2)*cos(q3))*d4+sin(q1)*cos(q2)*a3*cos (q3)-sin(q1)*sin(q2)*a3*sin(q3)+sin(q1)*a2*cos(q2)+a1*s in(q1)-f2,(-(-sin(q2)*cos(q3)-cos(q2)*sin(q3))*cos(q4)*sin(q5)+(-sin(q2)*sin(q3)+cos(q2)*cos(q3))*cos(q5))*d6+(sin( q2)*sin(q3)-cos(q2)*cos(q3))*d4-sin(q2)*a3*cos(q3)-cos(q2)*a3*sin(q3)-a2*sin(q2)+d1-f3,exp(q6/pi)*(-((cos(q1)*cos(q2)*cos(q3)-cos(q1)*sin(q2)*sin(q3))*cos(q4)+sin(q1)*sin(q4))* sin(q5)+(cos(q1)*cos(q2)*sin(q3)+cos(q1)*sin(q2)*c os(q3))*cos(q5))-f4,exp(q6/pi)*(-((sin(q1)*cos(q2)*cos(q3)-sin(q1)*sin(q2)*sin(q3))*cos(q4)-cos(q1)*sin(q4))*sin(q5)+(sin(q1)*cos(q2)*sin(q3)+ sin(q1)*sin(q2)*cos(q3))*cos(q5))-f5,-(-sin(q2)*cos(q3)-cos(q2)*sin(q3))*cos(q4)*sin(q5)+(-sin(q2)*sin(q3)+cos(q2)*cos(q3))*cos(q5)-f6
>> [q1 q2 q3 q4 q5 q6]=solve(equ)
و هر بار با اين پيغام مواجه ميشم.لطفا كمكم كنيد.
Error in ==> C:\MATLAB6p5\toolbox\symbolic\maplemex.dll
Error in ==> C:\MATLAB6p5\toolbox\symbolic\maple.m
On line 116 ==> [result,status] = maplemex(statement);
Error in ==> C:\MATLAB6p5\toolbox\symbolic\solve.m
On line 115 ==> [R,stat] = maple('solve',eqns,vars);:41:
ما رو سر کار گذاشتی؟
اخه این چیه نوشتی؟
معادله است؟
طرف اول و دومش کجان؟
یه ذره بهتر بنویس تا بتونیم کمکت کنیم.
خداییش اگه خودت یه همچین پستی می دیدی حالت گرفته نمی شد؟
در ضمن تو خط 6 ام پاراگراف دوم به جای sin نوشتی s in
تازه من اينو از متلب كپي كردم اينجا گذاشتم اگه تو متلب paste مي كردي سر كار نمي رفتي.
من ميخوام 6 معادله 6 مجهول در متلب حل كنم كه مجهولاتم شامل sin,cos,exp هستش ولي هربار متلب بهم error ميده و ميگه اشكال از maplemex.dll و maple وmaplemex.m است كسي ميتونه بهم كمك كنه در ضمن براي حل معادلاتم از تابع solve استفاده كردم.
لطفا هركي ميتونه فورا كمكم كنه.:45:
من ميخوام براي تحليل رباتم، 6 معادله 6 مجهول در متلب حل كنم كه مجهولاتم شامل sin,cos,exp هستش ولي هربار متلب بهم error ميده و ميگه اشكال از maplemex.dll و maple وmaplemex.m است كسي ميتونه بهم كمك كنه در ضمن براي حل معادلاتم از تابع solve استفاده كردم.
لطفا هركي ميتونه فورا كمكم كنه.
با سلام خدمت شما دوست عزیز ..
عنوانی که شما انتخاب کردید خلاف قوانین هست
[ برای مشاهده لینک ، با نام کاربری خود وارد شوید یا ثبت نام کنید ]
4- حتما براي طرح سوال خود در تاپيك جديد عنواني گويا انتخاب نماييد و از انتخاب عباراتي نظير «كمك!!، يك سوال داشتم و ...» خودداري نماييد.
لطفا عنوان مناسبی برای مطلب خود انتخاب کنید
معادلش رو بذار
اگه تو notepad كپي كنيد درست مي شود.
از اينكه معادلات طولانيه تعجب نكنيد ، مشكل خطاش حل شد ولي وقتي اجرا ميشه متلب براي مدت طولاني مشغوله تا 14 ساعت هم شد ولي جوابي نداد.من اينطوري پيش رفتم اگه ميشه تو متلب كپي كن و اجرا كني مي بيني كه چقدر طول مي كشه.
syms q1 q2 q3 q4 q5 q6 f1 f2 f3 f4 f5 f6
eq1=-13/40*((cos(q1)*cos(q2)*cos(q3)-cos(q1)*sin(q2)*sin(q3))*cos(q4)+sin(q1)*sin(q4))* sin(q5)+13/40*(cos(q1)*cos(q2)*sin(q3)+cos(q1)*sin(q2)*cos(q3 ))*cos(q5)-41/50*cos(q1)*cos(q2)*sin(q3)-41/50*cos(q1)*sin(q2)*cos(q3)+29/200*cos(q1)*cos(q2)*cos(q3)-29/200*cos(q1)*sin(q2)*sin(q3)+17/20*cos(q1)*cos(q2)+7/20*cos(q1)-f1;
eq2=-13/40*((sin(q1)*cos(q2)*cos(q3)-sin(q1)*sin(q2)*sin(q3))*cos(q4)-cos(q1)*sin(q4))*sin(q5)+13/40*(sin(q1)*cos(q2)*sin(q3)+sin(q1)*sin(q2)*cos(q3 ))*cos(q5)-41/50*sin(q1)*cos(q2)*sin(q3)-41/50*sin(q1)*sin(q2)*cos(q3)+29/200*sin(q1)*cos(q2)*cos(q3)-29/200*sin(q1)*sin(q2)*sin(q3)+17/20*sin(q1)*cos(q2)+7/20*sin(q1)-f2;
eq3=163/200-13/40*(-sin(q2)*cos(q3)-cos(q2)*sin(q3))*cos(q4)*sin(q5)+13/40*(-sin(q2)*sin(q3)+cos(q2)*cos(q3))*cos(q5)+41/50*sin(q2)*sin(q3)-41/50*cos(q2)*cos(q3)-29/200*sin(q2)*cos(q3)-29/200*cos(q2)*sin(q3)-17/20*sin(q2)-f3;
eq4=exp(q6/pi)*(-((cos(q1)*cos(q2)*cos(q3)-cos(q1)*sin(q2)*sin(q3))*cos(q4)+sin(q1)*sin(q4))* sin(q5)+(cos(q1)*cos(q2)*sin(q3)+cos(q1)*sin(q2)*c os(q3))*cos(q5))-f4;
eq5=exp(q6/pi)*(-((sin(q1)*cos(q2)*cos(q3)-sin(q1)*sin(q2)*sin(q3))*cos(q4)-cos(q1)*sin(q4))*sin(q5)+(sin(q1)*cos(q2)*sin(q3)+ sin(q1)*sin(q2)*cos(q3))*cos(q5))-f5;
eq6=exp(q6/pi)*(-(-sin(q2)*cos(q3)-cos(q2)*sin(q3))*cos(q4)*sin(q5)+(-sin(q2)*sin(q3)+cos(q2)*cos(q3))*cos(q5))-f6;
[q1,q2,q3,q4,q5,q6]=solve(eq1,eq2,q3,eq4,eq5,eq6)
می خواستم بدونم توی نرم افزارهای رسم نمودار چه طور میشه از جز صحیح و قدر مطلق و لگاریتم استفاده کرد
معادلتون همونیه که گفتید ؟
نمی دونم منظور شما کدوم نرم افزاره ولی تو
maple
برای جزء صحیح از
floor
برای قدر مطلق از
abs
و برای لوگاریتم از
ln
استفاده میشه. این هم مثال
[html]
plot(floor(x),x=-10..10);
plot(abs(x),x=-10..10);
plot(ln(x)/ln(3),x= 0.1..5);
[/html]
این آخری لگاریتم به پایه ی 3.
سلام
هر کسی در این مورد میتونه کمکی کنه خواهشا بگه ...
خیلی خیلی عجله دارم ....
ممنون از همه دوستان ................
يه كم توضيح بيشتر مي دي؟
من تا حالا 20 تا تابع چگالي خوندم. حجمي ، خطي ، شارشي و ....
كتابش رو هم بگي خوبه
منحنی چگالی توزیع (t(n را برای nهای مختلف ؛
n = 5, 10 , 20 , 30 , 50
رسم کنید و رفتار منحنی را برای nهای بزرگ تحلیل کنید ؛
با استفاده از قضیه 12.8 صفحه 303 از کتاب آمار ریاضی جان فروند ،
.............
کسی جوابمو نداد تا اینکه خودم به جواب رسیدم ...
برای اونایی میگم که میان و چیزی نمیگن ...
برای رسم منحنی کافیست اینطور بنویسیم ؛
[{Plot [f,{x,a,bf : تابع شما
x : متغییر تابع
a : ابتدای بازه
b : انتهای بازه
سري فوريه در matlab لطفا
:20: سری فوریه ،معادله یک موج پیچیده رو تبدیل به هارمونیکهای اون میکنه . موضوع زیر رو از help نرم افزار متلب گرفتم
fourier
--- help for sym/fourier.m ---
FOURIER Fourier integral transform.
F = FOURIER(f) is the Fourier transform of the sym scalar f
with default independent variable x. The default return is
a function of w.
If f = f(w), then FOURIER returns a function of t: F = F(t).
By definition, F(w) = int(f(x)*exp(-i*w*x),x,-inf,inf), where
the integration above proceeds with respect to x (the symbolic
variable in f as determined by FINDSYM).
F = FOURIER(f,v) makes F a function of the sym v instead of
the default w:
FOURIER(f,v) <=> F(v) = int(f(x)*exp(-i*v*x),x,-inf,inf).
FOURIER(f,u,v) makes f a function of u instead of the
default x. The integration is then with respect to u.
FOURIER(f,u,v) <=> F(v) = int(f(u)*exp(-i*v*u),u,-inf,inf).
Examples:
syms t v w x
fourier(1/t) returns i*pi*(Heaviside(-w)-Heaviside(w))
fourier(exp(-x^2),x,t) returns pi^(1/2)*exp(-1/4*t^2)
fourier(exp(-t)*sym('Heaviside(t)'),v) returns 1/(1+i*v)
fourier(diff(sym('F(x)')),x,w) returns i*w*fourier(F(x),x,w)
See also IFOURIER, LAPLACE, ZTRANS.
Reference page in Help browser
اگه میخواهید که جزییات سری فوریه رو ببینید . از تئوری تا حل چند مثال ، در محیط متلب تایپ کنید:
doc fourier بطور دقیق جزییات رو بخونید
سری فوریه بیان می کنه که توابع متناوب را می توان به صورت مجموعی از توابع متناوب نوشت (شایدم برعکس)نقل قول:
[ برای مشاهده لینک ، با نام کاربری خود وارد شوید یا ثبت نام کنید ] سری فوریه ،معادله یک موج پیچیده رو تبدیل به هارمونیکهای اون میکنه . موضوع زیر رو از help نرم افزار متلب گرفتم
که بعدش نتیجه میشه مثلا موج صدا مجموعی از بینهایت موج تخت سینوسی می باشد
دوستان کسی میتونه منو تو رسم تابع زیر کمک کنه؟
[ برای مشاهده لینک ، با نام کاربری خود وارد شوید یا ثبت نام کنید ]
ببخشید
این cte یعنی چی؟
اگه a عددی ثابت هست ، اول تابع یا رابطه رو بر حسب x مرتب کن بعد با نرم افزار مطلب رسمش کن
X=y+( a^2/19.612) y^-2
y مخالف صفر هست . از x نسبت به y مشتق میگیریم :
dx/dy=1+( 2a^2/19.612) y^-3=0
( 3/1 ) ^(y^3=(2a^2/19.612)--------y=(2a^2/19.612
این از اکسترمم حالا یک بار دیگه از 'y مشتق بگیر برابر صفر قرار بده 0(="y )
نقطه عطف پیدا میشه .....................
البته در تابع x بر حسب y رسم میشه
در نرم افزار حتما مقدار a رو در اول برنامه مشخص کن
راستش من کار با مطلب رو خیلی بلد نیستم
میشه یه آموزش مقدماتی برای رسم این تابع برام بزاری:11:
راستی من میخوام y بر حسب x رسم بشه.
cte هم به معنی ثابت هست.نقل قول:
ok
سعی خودمو میکنم.عجله ای که ندارید؟؟؟
خب تابع رو بصورت زیر مرتب کنید بصورت زیر :
xy^2-y^3=(a^2/19.612)=C
سمت راست معادله رو ضریب ثابت گرفتم . تا اینجا رو داشته باشید
یه پیشنهاد بهتر
همون تابع x بر حسب y رو که داری اگه جای x, y رو عوض کنیتابع معکوس بدست میاد . حالا تابع معکوس رو که یک تابع یکدست هست (x یک طرف و y طرف دیگه ) تا آخر رسم کن . از اونجایی که هر تابع با معکوسش نسبت به نیمساز ربع فکر کنم اول و سوم قرینه هستند ، قرینه تابع رسم شده نسبت به نیمساز ربع اول و سوم (y=x) جواب شما رو میده
فکر کنم ابهامی نباشه
ببین a=1 گرفتم
شده این
[ برای مشاهده لینک ، با نام کاربری خود وارد شوید یا ثبت نام کنید ]
البته با maple
به ازای aهای مختلف:
کد:http://ufile.info/up/files/ae9yev9kv9fooylsba21.png
کد:http://ufile.info/up/files/txs58vha6nphka3bye2p.png
راستی یادم رفت بگم که این یک تابع هموگرافیک هست هر چه a رو بیشتر به سمت صفر میل بدیم ، شکل منحنی به خط y=x نزدیکتر میشه . در حالت حدی به جای a=0 قرار بدبد y=x حاصل میشه . منحنی زیر به ازای a=1,0.5,0.25,0.05 رسم شده. ملاحظه کنید که با کاهش a شکل منحنی به نیمساز ربع اول وسوم y=x نزدیکتر میشه .
رسم تابع معکوس با استفاده از matlab
این برنامه فوق العاده ساده رو به مطلب دادم ، منحنی ها رو رسم کرد
clear
a=1;
c=(a^2)/19.612;
k=0;
for i=-3:0.01:3
k=k+1;
y(k)=i;
x(k)=y(k)+(c/((y(k))^2));
end
plot(x,y)
برای a=1 و y در بازه [3,3-]:
************************************************** ******
a=0.5 کمتر از 1:
************************************************** *******************************
a=0.25 :
*********************************
a=0.05 و y در بازه [1,1-]
فقط داده های برنامه رو تغییر دادم . به همین سادگی:
clear
a=0.05;
c=(a^2)/19.612;
k=0;
for i=-1:0.01:1
k=k+1;
y(k)=i;
x(k)=y(k)+(c/((y(k))^2));
end
plot(x,y)
کار کردن با matlab خیلی ساده هست
دوستان من الان بخوام ریشه ی این معادله رو با مطلب بدست بیارم همچنین نتیجه ای میده:
[ برای مشاهده لینک ، با نام کاربری خود وارد شوید یا ثبت نام کنید ]
این به چه معنیه؟
فکر کنم
البته همه اش را باید خودتون بدونید
معادله به درجه اش ریشه داره. این یک ریشه اش حقیقی و 2 ریشه اش (نمی دونم چرا اسمش یادم نمیاد - آهها موهومی)
مثل x^2+1=0 که دو ریشه اش موهومی است
اون i هم که رادیکال منفی یک است
در ضمن اعداد موهومی به شکل a+bi نوشته می شوند