Matlab denklem çözme

Katılım
23 Temmuz 2019
Mesajlar
2.785
En iyi cevaplar
11
Beğeniler
1.317
Arkadaşlar fiziksel diferansiyel sistemlerinin çözümüne dair bir ödevim var. Ben yazıyorum ama hata alıyorum. Neyi eksik tanımlıyorum ya da hatam nerede?
Dipnot: Değerler soruya ait değil. Daha doğrusu hoca değer vermemiş. Değerler bizim tasarrufumuzda.
Denklem:

1591147544048-png.581073


Çözüm:
>> M=10;
>> K=1;
>> C=0.2;
>> Dsolve('m*D2X(t)+C*DX(t)+k*X(t)=F(X)','F(0)=0','df(0)=0')
>> Ezplot (X(t),[0,10])

1591148036343-png.581079
 
Son düzenleyen: Moderatör:
Katılım
29 Nisan 2016
Mesajlar
89
En iyi cevaplar
1
Beğeniler
18
Yer
Yozgat
'dsolve' kullanımı için böyle bir yardım sayfası mevcut aynı zamanda aşağıdaki t değişkenini belli aralıkta tanımlaman gerekmiyor mu?
Mesela 't = 1:10' gibi.
 
KS
fenerlimühendis
Katılım
23 Temmuz 2019
Mesajlar
2.785
En iyi cevaplar
11
Beğeniler
1.317
Katılım
29 Nisan 2016
Mesajlar
89
En iyi cevaplar
1
Beğeniler
18
Yer
Yozgat
Kod:
m=10;
k=1;
c=0.2;
syms a x(t)
syms b x(t)

eqn1 = diff(x,t)==a*x;

eqn2 = diff(x,t,2)==b*x;

eqn = m*eqn2+c*eqn1+1*x;
a = dsolve(eqn);
Merhaba biraz uğraştım fakat çıktı alamadım.

a(1) yazdırıldığında şu şekilde çıktı aldım;

C1*exp(-t*((200*a + 10000*b + 1)^(1/2)/100 + 1/100)) + C2*exp(t*((200*a + 10000*b + 1)^(1/2)/100 - 1/100))

Merhaba eğer bir şey bulabildiysen buraya yazarsan sevinirim.
Kod:
syms x(t)
syms f(t)

dsolve('10*D2x(t)+0.2*Dx(t)+1*x(t)=f(t)','f(0)=0','Df(0)=0')
fplot(f(t),[0,10])
Son olarak böyle birşey yaptım çıkış vermedi.
 
Son düzenleme:
KS
fenerlimühendis
Katılım
23 Temmuz 2019
Mesajlar
2.785
En iyi cevaplar
11
Beğeniler
1.317
Kod:
m=10;
k=1;
c=0.2;
syms a x(t)
syms b x(t)

eqn1 = diff(x,t)==a*x;

eqn2 = diff(x,t,2)==b*x;

eqn = m*eqn2+c*eqn1+1*x;
a = dsolve(eqn);
Merhaba biraz uğraştım fakat çıktı alamadım.

a(1) yazdırıldığında şu şekilde çıktı aldım;

C1*exp(-t*((200*a + 10000*b + 1)^(1/2)/100 + 1/100)) + C2*exp(t*((200*a + 10000*b + 1)^(1/2)/100 - 1/100))

Merhaba eğer bir şey bulabildiysen buraya yazarsan sevinirim.
Kod:
syms x(t)
syms f(t)

dsolve('10*D2x(t)+0.2*Dx(t)+1*x(t)=f(t)','f(0)=0','Df(0)=0')
fplot(f(t),[0,10])
Son olarak böyle birşey yaptım çıkış vermedi.
Hocam olmuyor ya :/
 
KS
fenerlimühendis
Katılım
23 Temmuz 2019
Mesajlar
2.785
En iyi cevaplar
11
Beğeniler
1.317
t'yi nasıl belirtebilirim hocam?
 

yazicizadee

Nanopat
Katılım
3 Aralık 2019
Mesajlar
159
En iyi cevaplar
2
Beğeniler
58
T yi vermem lazımki grafik çıksın. Çizimde zaten F(t).O t yi vermeden denklem çözülmezki zaten. T ye bir değer ver 10 gibi grafik çıkıyor mu bak. Ayrıca Matlab bütün değerleri saklar ve gösterir. Garip gördüğünüz yerlere bir daha bakın.
T yi vermem lazımki grafik çıksın. Çizimde zaten F(t).O t yi vermeden denklem çözülmezki zaten. T ye bir değer ver 10 gibi grafik çıkıyor mu bak. Ayrıca Matlab bütün değerleri saklar ve gösterir. Garip gördüğünüz yerlere bir daha bakın.
Hiç matlab bilmiyorsan başa t=10; yaz
 
KS
fenerlimühendis
Katılım
23 Temmuz 2019
Mesajlar
2.785
En iyi cevaplar
11
Beğeniler
1.317
1591228019096-png.582018

Hocam doğru yazdığıma hemen hemen eminim. Ama çıkmıyor.
 

yazicizadee

Nanopat
Katılım
3 Aralık 2019
Mesajlar
159
En iyi cevaplar
2
Beğeniler
58
Öncelikle script olarak yazmanızı öneriyorum.Denklemde çok bilinmeyen var diyor sanki
 


Yeni konular

Yukarı