یکی کمکم کنه تا یه این پروژه رو تموم کنم.

[amirr13]

خوب این سوال اول

من یه سری برنامه نوشتم تا اینکه یه سه معادله سه مجهول رو برام حل کنه

حالا برای اینکه از این سه جواب تو ادامه محاسبات استفاده کنم چی باید بکنم؟

الان یه k دارم که میخوام یه جوری حالیش کنم خروجیم به صورت k1 , k2 , k3 باشه

فعلا همین

[666mhd]

پاسخ : یکی کمکم کنه تا یه این پروژه رو تموم کنم.

منظورتون از k1 k2 k3 را متوجه نمیشم
برنامه را بگذارید ببینم
شاید بتونم مشکلتون را حل کنم

[amirr13]

پاسخ : یکی کمکم کنه تا یه این پروژه رو تموم کنم.

>> A=[a1,a2,a3;b1,b2,b3;c1,c2,c3];
>> B=[a4;b4;c4];
>> k=A\B;

متغیر ها میریزن تو دو تا ماتریس که حلش میشه این

حالا ایم متلب چجوری نشون میده؟
این جوری

k =

1.0131
-84.7058
-84.6894

اما من میخوام به صورت K1=1.0131 , k2=-84.7058 ,k3=-84.6894

باشه چون تو ادامه باید روی این k ها باز عملیات انجام بشه و این اعداد جواب نهایی برنامه نیستند

ممنون میشم اگه بتونی کمک کنی

[amirr13]

پاسخ : یکی کمکم کنه تا یه این پروژه رو تموم کنم.

خوب بی خیال اونو مجبور شدم با کرامر حل کنم

اینو کسی میتونه جواب بده؟

با حلقه فور 5 تا جواب بدست آوردم

چجوری باید اونارو رو یه نمودار نشون بدم؟

[666mhd]

پاسخ : یکی کمکم کنه تا یه این پروژه رو تموم کنم.

خوب بی خیال اونو مجبور شدم با کرامر حل کنم

اینو کسی میتونه جواب بده؟

با حلقه فور 5 تا جواب بدست آوردم

چجوری باید اونارو رو یه نمودار نشون بدم؟

توی حلقه for قبل از این که حلقه ی بعدی شروع بشه اون جواب ها را نشون بده

برنامه را اگه بگذارید بهتره

[amirr13]

پاسخ : یکی کمکم کنه تا یه این پروژه رو تموم کنم.

for z=(0:0.001:pi/4)
f=f1+((df/dx)*(z-x1))
j=sin(f)
k=a1/a2-cos(f)
l=((a1^2+a2^2+a4^2-a3^2)/(2*a1*a4)-((a1/a4)*cos(f)))
S1=2*atan((j^2+sqrt(j^2+k^2-l^2))/k+l)
S2=2*atan((j^2-sqrt(j^2+k^2-l^2))/k+l)
if S1<=2.094395,S>=0
S=S1
else
S=S2
end
Y=((s1-S)*(dy/ds))+y1
T=sin(z)
disp(Y)
disp(T)
plot(Y,z,'pg&#039
plot(T,z)
end

اینم از حلقه for برنامه

ابنم کلش

>> xs=0;
xf=pi/4;
a=(xs+xf)/2;
h=(xf-xs)/2;
x1=a-h*cos(pi/6);
x2=a;
x3=a+h*cos(pi/6);
y1=sin(x1);
y2=sin(x2);
y3=sin(x3);
dx=xf-xs;
dy=sin(xf)-sin(xs);
df=pi;
ds=pi/3;
f12=(df/dx)*(x2-x1);
f13=(df/dx)*(x3-x1);
s12=(ds/dy)*(y2-y1);
s13=(ds/dy)*(y3-y1);
f1=pi/6;
s1=pi/6;
f2=f1+f12;
f3=f1+f13;
s2=s1+s12;
s3=s1+s13;
a1=cos(f1);
a2=cos(s1);
a3=-1;
a4=cos((f1-s1));
b1=cos(f2);
b2=cos(s2);
b3=-1;
b4=cos((f2-s2));
c1=cos(f3);
c2=cos(s3);
c3=-1;
c4=cos((f3-s3));
A=[a1,a2,a3;b1,b2,b3;c1,c2,c3];
B=[a4;b4;c4];
k=A\B;

M1=[a4,a2,a3;b4,b2,b3;c4,c2,c3];
M2=[a1,a4,a3;b1,b4,b3;c1,c4,c3];
M3=[a1,a2,a4;b1,b2,b4;c1,c2,c4];
k1=(det(M1)/det(A));
k2=det(M2)/det(A);
k3=det(M3)/det(A);
a1=10;
a2=a1/k2;
a4=a1/k1;
a3=sqrt((a1^2+a2^2+a4^2)-(2*a2*a4*k3));
b=[0,pi/10,pi/8,pi/6,pi/4];
for z=(0:0.001:pi/4)
f=f1+((df/dx)*(z-x1))
j=sin(f)
k=a1/a2-cos(f)
l=((a1^2+a2^2+a4^2-a3^2)/(2*a1*a4)-((a1/a4)*cos(f)))
S1=2*atan((j^2+sqrt(j^2+k^2-l^2))/k+l)
S2=2*atan((j^2-sqrt(j^2+k^2-l^2))/k+l)
if S1<=2.094395,S>=0
S=S1
else
S=S2
end
Y=((s1-S)*(dy/ds))+y1
T=sin(z)
disp(Y)
disp(T)
plot(Y,z,'pg&#039
plot(T,z)
end

اینو اگه گفتین ، اینم بگین که میشه فایل اجرایی ساخت با متلب؟

اگه نمیشه چی کنم که بتونم خارج از محیط متلب اجرا بشه

اینم نمیشه؟

چه خاکب تو سرش بریزم؟ :angry:

[amirr13]

پاسخ : یکی کمکم کنه تا یه این پروژه رو تموم کنم.

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

مجبور شدم نصفه تحویل بدم

شاد باشید