Matlab denklem çözme

feneresevdalı · 4 Haziran 2020

yazicizadee dedi:
Öncelikle script olarak yazmanızı öneriyorum.Denklemde çok bilinmeyen var diyor sanki

Hocam sizi alıkoyuyorum ama. Script olarak nasıl yazabilirim?
Hocam
dsolve('m*D2x(t)+c*Dx(t)+k*x(t)=f(t)','f(0)=0') yazınca

ans =

piecewise([f(0) == 0, {C6*exp(-(t*(c + (c^2 - 4*k*m)^(1/2)))/(2*m)) + exp(-(t*(c - (c^2 - 4*k*m)^(1/2)))/(2*m))*int((exp((x*(c - (c^2 - 4*k*m)^(1/2)))/(2*m))*f(x))/(c^2 - 4*k*m)^(1/2), x, 0, t, 'IgnoreAnalyticConstraints', true) + C5*exp(-(t*(c - (c^2 - 4*k*m)^(1/2)))/(2*m)) + exp(-(t*(c + (c^2 - 4*k*m)^(1/2)))/(2*m))*int(-(exp((x*(c + (c^2 - 4*k*m)^(1/2)))/(2*m))*f(x))/(c^2 - 4*k*m)^(1/2), x, 0, t, 'IgnoreAnalyticConstraints', true)}], [f(0) ~= 0, {}])

çıktısını aldım. Doğru mu?

yazicizadee · 4 Haziran 2020

feneresevdalı dedi:
Hocam sizi alıkoyuyorum ama. Script olarak nasıl yazabilirim?
Hocam.
Dsolve('m*D2X(t)+C*DX(t)+k*X(t)=F(t)','F(0)=0') yazınca.

Ans =

Piecewise([F(0) == 0, {C6*exp(-(t*(C + (C^2 - 4*k*m)^(1/2)))/(2*m)) + exp(-(t*(C - (C^2 - 4*k*m)^(1/2)))/(2*m))*int((exp((X*(C - (C^2 - 4*k*m)^(1/2)))/(2*m))*F(X))/(C^2 - 4*k*m)^(1/2), X, 0, T, 'ıgnoreanalyticconstraints', true) + C5*exp(-(t*(C - (C^2 - 4*k*m)^(1/2)))/(2*m)) + exp(-(t*(C + (C^2 - 4*k*m)^(1/2)))/(2*m))*int(-(exp((X*(C + (C^2 - 4*k*m)^(1/2)))/(2*m))*F(X))/(C^2 - 4*k*m)^(1/2), X, 0, T, 'ıgnoreanalyticconstraints', true)}], [f(0) ~= 0, {}])

Çıktısını aldım. Doğru mu?

Buna basarak yapın her türlü Matlab hesaplamaları diğer türlü hesap makinesinden farksız olur.

feneresevdalı · 4 Haziran 2020

yazicizadee dedi:
Eki Görüntüle 582020

Buna basarak yapın her türlü Matlab hesaplamaları diğer türlü hesap makinesinden farksız olur.

Hocam hatasız olduğunu biliyorum. Demek istediğim 'Df(0)=0' şartını sildim. Çıktı verdi. Güvenebilir miyim? Bir de Ezplot ifadesini nasıl kullanabilirim?

yazicizadee · 4 Haziran 2020

Plot expression or function - MATLAB fplot

This MATLAB function plots the curve defined by the function y = f(x) over the default interval [-5 5] for x.

www.mathworks.com

.
Ezplot çok eski kaç sürümünü kullanıyorsunuz. Yani Matlab'in kendisi aşırı yardım ediyor. Açıkcası tüm işimi oraya bakarak hallediyorum. Ezbere bilmem.

feneresevdalı dedi:
Hocam hatasız olduğunu biliyorum. Demek istediğim 'df(0)=0' şartını sildim. Çıktı verdi. Güvenebilir miyim? Bir de ezplot ifadesini nasıl kullanabilirim?

Cevabı buldun mu?

feneresevdalı · 4 Haziran 2020

Hocam bana çok güvenilir gelmedi cevap. Mühendislik uygulamalarında kullanılan 3 tane diferansiyel denklemin çözümüne dair matlab kodu isteniyor bizden. Matlab'i şimdiye kadar hep lineer cebir için kullanıyordum. Bana 1 tane denklem+kod örneği verirseniz onun üstünde çalışır diğerlerini yaparım.

Şu şartlar altında yazabilir misiniz?

Sanki bu sefer oldu gibi. Ne dersiniz hocam doğru mudur?

Hocam bir de bunun plotunu yazabilir misiniz?

Ahmet Hüdai Kaya · 4 Haziran 2020

Çıktıyı bu şekilde aldım;

Kod:

syms x(t) f(t) t;
m=0.25;
c=0.5;
k=1;
solution = dsolve('10*D2x(t)+0.2*Dx(t)+1*x(t)=2','x(0)=0','Dx(0)=0');
t = 0:10;
fplot(x(t),solution)

Verdiği çıktı ise bu.

feneresevdalı · 4 Haziran 2020

Ahmet Hüdai Kaya dedi:
Çıktıyı bu şekilde aldım;

Kod:

syms x(t) f(t) t; m=0.25; c=0.5; k=1; solution = dsolve('10*D2x(t)+0.2*Dx(t)+1*x(t)=2','x(0)=0','Dx(0)=0'); t = 0:10; fplot(x(t),solution)

Eki Görüntüle 582079

Verdiği çıktı ise bu.

Hocam attığınız kod bende çalışmadı.

Ahmet Hüdai Kaya · 4 Haziran 2020

feneresevdalı dedi:
Hocam attığınız kod bende çalışmadı.

Hata verdi mi?

feneresevdalı · 4 Haziran 2020

Hocam bu olmuş mu?

Ahmet Hüdai Kaya · 4 Haziran 2020

Katsayıları aynı hale getirince hemen hemen aynı çıkışı elde ediyoruz. 'ezplot' ve 'fplot' sadece sürüm farkından sıkıntı çıkarıyor.

Matlab denklem çözme

feneresevdalı

Megapat

yazicizadee

Decapat

feneresevdalı

Megapat

yazicizadee

Decapat

Plot expression or function - MATLAB fplot

feneresevdalı

Megapat

Ahmet Hüdai Kaya

Kilopat

feneresevdalı

Megapat

Ahmet Hüdai Kaya

Kilopat

feneresevdalı

Megapat

Ahmet Hüdai Kaya

Kilopat