Printable View
با سلام به دوستان گرامی،
من یه تابع درجه دوم دارم که باید ازش مشتق بگیرم و بعد برابر یک ضریبی قرار بدم و سپس ریشه ی اون تابع مشتق گرفته شده (درجه یک) رو بدست بیارم.
من مشتق میگیرم و برابر اون ضریب قرار میدم ولی نمیتونم ریشه رو بدست بیارم.
کسی هست کمک کنه؟
منظورم اینه که دستور این معادله رو چطوری بنویسم؟ (x متغیره)
2ax+b=c
باید از این معادله x بدست بیاد.
مرسی
سلام
اگه می خواهید به صورت سمبولیک حل کنید می تونید از تولباکس symbolic استفاده کنید
نه. فعلن میخوام کد بزنم!نقل قول:
نوشته شده توسط rahnema1
دستوری نداره؟
ببینید این لینک به کارتون میاد:نقل قول:
نوشته شده توسط coronaa
http://www2.math.umd.edu/~immortal/2...quations.shtml
حالا اگه می خواهید چرخ را دوباره اختراع کنید باید کاری شبیه کاری که مثلا یک کامپایلر انجام میده انجام بدین
ممنونم. دقیقاً همینو میخواستم!نقل قول:
نوشته شده توسط rahnema1
مهندس مشکلم تقریباً حل شد. حالا میخوام نتیجه رو چاپ کنم یه پیغام خطا میده!
clc; clear;
syms P
for i=1:5
a(i)=5;
b(i)=10;
c(i)=15;
ro(i)=30;
lambda=20;
% -----------------------------
f=a*P^2+b*P+c;
d=diff(f,'P');
% -----------------------------
s=d(1,i)-lambda;
p(1,i)=solve(s) ; % P(i) bozorg
Pr=sum(p);
% -----------------------------
end
fprintf('Pi(Big) = %f',p);
fprintf('Pr = %f',Pr);
این خطا رو میده. ظاهراً P که بعنوان سیمبول شناخته شد دیگه دستور fprintf کار نمیکنه :
??? Error using ==> fprintf
Function is not defined for 'sym' inputs.
Error in ==> ProftCode at
fprintf('Pi(Big) = %f',p);
من تولباک سمبولیک ندارم حالا شما اگه بتونید اول p را تبدیل به عدد اعشاری بکنید چون f% که در fprintf استفاده میشه برای چاپ عدد اعشاری هست. مثلا با استفاده از دستور double یا امثال اینها
p هم عدد اعشاریه هم یک ماتریس یک در 5 هست. نکنه بخاطر ماتریس بودنش قبول نمیکنه؟نقل قول:
نوشته شده توسط rahnema1
ممنونم دوست عزیز.
حل شد تا اینجا. :لبخندساده:
اگه یه مجموعه داشته باشیم بصورت {1,2,3} و بخوایم حالت های ممکنی که می تونن با هم یه زیر مجموعه تشکیل بدن چطوریه؟
مثلاً به این صورت :
{1},{1,2},{1,3},{2,3},{1,2,3}
ضمن اینکه {1,2} و {2,1} فرقی نمیکنن با هم!
a=[1 2 3];
for i=1:numel(a)
s{i}=nchoosek(a,i);
end
disp(s)
nchoosek، جایگشته؟نقل قول:
نوشته شده توسط rahnema1
ولی من تمام حالتهای زیر مجموعه رو میخوام. تعدادشون نه! خود حالت ها رو نشون بده.
خروجی متلب :
[3x1 double] [3x2 double] [1x3 double]
s سک سلوله که داخلش تمام حالتها قرار داره.
خروجی octave:
octave-3.6.2.exe:3> disp(s)
{
[1,1] =
1
2
3
[1,2] =
1 2
1 3
2 3
[1,3] =
1 2 3
}
دستتون درد نکنه ولی من هنوز اون جوابی که میخواستم رو نگرفتم.
یکمی باز مقالاتم رو مطالعه میکنم تا برام روشنتر بشه که چی میخوام. بعد دوباره سوالم رو از شما میپرسم.
ممنونم از راهنمائیت.
مرسی
آقای مهندس،
برنامه ای که دادی رو یکمی دستکاری کردم و به یه جوابی رسیدم. بصورت زیر :
;
a=[1 2 3];
for i=1:numel(a);
v=nchoosek(a,i)
end
خروجی متلب :
v =
1
2
3
v =
1 2
1 3
2 3
v =
1 2 3
ولی میخوام خروجی رو به این صورت نمایش بده :
http://www.uplooder.net/img/image/1/...89cc88/aaa.jpg
راه حلی وجود داره به نظر شما؟
این را امتحان کنیدنقل قول:
راه حلی وجود داره به نظر شما؟
a=[1 2 3];
k=1;
for i=1:numel(a)
s=nchoosek(a,i);
for j=1:size(s,1)
v(k)=sum(x(s(j,:)).*y(s(j,:)));
k=k+1;
end
end
disp(v)
خطا :
??? Undefined function or method 'x' for input arguments of type
'double'.
Error in ==> sett1 at 6
v(k)=sum(x(s(j,:)).*y(s(j,:)));
توی فایل ورد که x و y داشتین. من فکر کردم x و y عبارتند از دو تا بردار مثلا با طول 3 (چون a هم از 1 تا 3 بود)یه چیزی شبیه این
x=rand(1,3);
y=rand(1,3);
بله. میشه گفت برداره.نقل قول:
نوشته شده توسط rahnema1
در اصل توی رابطه ام x,y وجود نداره. رابطه های دیگه هست ولی گفتم پیچیده نشه فقط بعنوان مثال یه جمله ی اون رو گذاشتم.
خب حالا جوابی که گرفتم بصورت زیره :
Columns 1 through 6
0.7441 0.5728 0.0124 1.3169 0.7565 0.5852
Column 7
1.3293
حالا v برای حالتهای مختلف به دست اومده
چطوری باید بفهمم کدوم v هست؟نقل قول:
نوشته شده توسط rahnema1
به همون ترتیبی که s را چاپ می کردیم
ادامه ش اینطوریه :
{i} یعنی یکی از اعضای s که در اینجا یکی از اعداد 1 و 2 یا 3 میتونه باشه
به. یعنی باید عضو i از اونا کم بشه.نقل قول:
نوشته شده توسط rahnema1
اگه i برابر 1 هست میشه :
v(1,2)-v(2)
v(1,2,3)-v(2,3)
فکر می کنم این فی باید جزء یک فرمول بزرگتر باشه درسته؟
نه.نقل قول:
نوشته شده توسط rahnema1
مستقله. فقط همون Vs هایی که بدست آوردیم باید درش قرار بدیم تا فی i بدست بیاد.
آخه یه مشکل وجود داره. توی فرمول v که در فایل ورد اولی بود اون اندکس i منظورش چیه؟ منظورش اعضای مجموعه s هست یا نه صرفا اندیس هست؟
تعداد اعضای مجموعه A هست. یعنی 3 تانقل قول:
نوشته شده توسط rahnema1
حل مسئله بدین صورته :
این یکی
a=[1 2 3];
x=rand(1,3);
y=rand(1,3);
function ret= indexV(Smini,S,v)
ret=0;
if (numel(Smini)==0)
return;
end
for i=1:numel(S)
if(size(S{i},2)<numel(Smini))
ret=ret+size(S{i},1);
else
[val,dx]=ismember(Smini, S{i},'rows');
ret=ret+dx;
break;
end
end
ret=v(ret);
end
k=1;
S={};
v=[];
for i=1:numel(a)
s=nchoosek(a,i);
S{i}=s;
for j=1:size(s,1)
v(k)=sum(x(s(j,:)).*y(s(j,:)));
k=k+1;
end
end
k=1;
phi=zeros(1,numel(a));
for i=1:numel(S)
for j=size(S{i},1)
for z=1:numel(a)
if (ismember(a(z),S{i}(j,:)))
phi(z)=phi(z)+factorial(numel(S{i}(j,:))-1)* factorial(numel(a) -numel(S{i} (j,:)))* (v(k)-indexV(setdiff( S{i}(j,:) ,a(z)),S,v))/ factorial(numel(a));
k=k+1;
end
end
end
end
disp(phi)
دستت درد نکنه. این خطا رو میده فعلاً :
??? Undefined function or variable 's'.
Error in ==> rphii at 9
S{i}=s;
آقا فعلاً خطا نداد. دارم بررسی میکنم ببینم جواب درست میده یا نه!
مرسی مهندس :تشویق::تشویق:
مهندس، من رابطه های خودمو توی برنامه جایگذاری کردم.
لینکی که براتون فرستادم ، رابطه ی دوم باید توی متلب نوشته بشه. من بصورت زیر نوشتم. چند تا جمله ی دیگه هم هست. این جمله رو گیر میده میگه ماتریس باید مربعی باشه. یعنی سایز این ماتریس با ماتریسهای دیگه که با هم جمع و کم میشن یکی نیست.
مشکلش کجاست ؟
http://www.uplooder.net/img/image/47...a6999/zzzz.jpg
v(k)=sum(a(s(j,:)).*(pi(s(j,:)))^2+b(s(j,:)).*pi(s (j,:))+c(s(j,:)));
طرز نوشتنش درسته؟(این a با a که قبلن تعریف کردیم فرق میکنه.
a,b,c,pi هر کدام، یک بردار یک در سه هستن.
آقا حل شد. ظاهرن!نقل قول:
نوشته شده توسط coronaa
بردار pi که به توان 2 بود رو ضرب نقطه ای کردم. درسته؟
v(k)=sum(a(s(j,:)).*(pi(s(j,:)).*pi(s(j,:)))+b(s(j ,:)).*pi(s(j,:))+c(s(j,:)));
حالا فعلن مشکلش اینه که v و فی رو منفی میده (البته با داده های خودم!)
دارم بررسیش میکنم.
تا اینجا خیلی چیزها رو از شما یاد گرفتم. واقعن سپاسگزارم.
مرسی
ببخشید اون قسمت آخر کد یه اشتباه کوچولو شده لطفا تغییر بدید
k=1;
phi=zeros(1,numel(a));
for i=1:numel(S)
for j=size(S{i},1)
for z=1:numel(a)
if (ismember(a(z),S{i}(j,:)))
phi(z)=phi(z)+factorial(numel(S{i}(j,:))-1)* factorial(numel(a) -numel(S{i} (j,:)))* (v(k)-indexV(setdiff( S{i}(j,:) ,a(z)),S,v))/ factorial(numel(a));
end
end
k=k+1;
end
end
disp(phi)
مهندس میشه به طور مختصر از تابعی که تعریف کردی رو توضیح بدی؟ احتمالن واسه محاسبه یs-{i}هست!
همچنین خط 38 برنامه.
ممنونم
بله تابع واسه s-{i هست. ا.ن سلول خودش شامل 3 تا زیر مجموعه هست ما تعداد اعضای این زیر مجموعه ها را می شماریم تا برسیم به عنصر مورد نظر و شماره اندکس اون را در میاریم تا بتونیم v با اون اندکس را پیدا کنیمنقل قول:
نوشته شده توسط coronaa
توی اون حلقه های تودرتوی آخر هم تقریبا همین روند را دنبال می کنیم با این تفاوت که الان ما v ها را داریم
یه اشتباه دیگه کشف کردم فکر کنم دیگه مشکلی نباشه:)
k=1;
phi=zeros(1,numel(a));
for i=1:numel(S)
for j=1:size(S{i},1)
for z=1:numel(a)
if (ismember(a(z),S{i}(j,:)))
phi(z)=phi(z)+factorial(numel(S{i}(j,:))-1)* factorial(numel(a) -numel(S{i} (j,:)))* (v(k)-indexV(setdiff( S{i}(j,:) ,a(z)),S,v))/ factorial(numel(a));
end
end
k=k+1;
end
end
disp(phi)
چرا فی منفی در میاد؟
همه ی v ها که مثبت هستن!
این اصلاح آخری که کردم دیگه با این داده ها منفی نمیشهنقل قول:
نوشته شده توسط coronaa