很多的工程問題中會有微積分的問題,我們需要對微分方程求解。
微分方程的求解分為符號微積分和數值微積分兩類。
本文介紹matlab中對微分方程進行符號求解
在matalb中對微分方程進行符號求解的時候使用的函式為int函式
呼叫格式:
f = int(expr) expr為函式f = int(expr,var) var為指定變數f = int(expr,a,b) a,b分別為下限和上限f = int(expr,var,a,b)f = int(___,name,value) 使用乙個或多個「名稱/值」對引數指定其他選項。例如,「 ignoreanalyticconstraints」,true指定int對被積數應用其他簡化
呼叫示例:
執行結果:%% f = int(expr)
syms x
expr = sin(x);
f1 = int(expr)
%% f = int(expr,var)
syms x t
expr = t*sin(x);
f2 = int(expr,x)
f3 = int(expr,t)
%% f = int(expr,a,b)
syms x
expr = sin(x);
f4 = int(expr,-0.5,0.3)
%% f = int(expr,var,a,b)
syms x t
expr = t*sin(x);
f5 = int(expr,x,[-0.5,0.3])
%% f = int(___,name,value)
syms f(x)
f(x) = acos(cos(x));
f = int(f,x)
f = int(f,x,'ignoreanalyticconstraints',true) %詳細參考下面解釋
此處對上面呼叫格式5做一解釋說明:預設情況下,int使用嚴格的數學規則。這些規則不允許int將acos(cos(x))重寫為x。 如果您想要乙個簡單實用的解決方案,請將「 ignoreanalyticconstraints」設定為true。f1 = -cos(x)
f2 = -t*cos(x),f3 = (t^2*sin(x))/2
f4 = cos(1/2) - cos(3/10)
f5 = t*(cos(1/2) - cos(3/10))
f1 = 1/2,f2 = 1/2
在執行上述程式的時候有疑問?程式執行發生了錯誤?
大家即解決了自己的問題,也會方便和你有同樣問題的同學學習交流。
matlab、python程式設計問題總結zhuanlan.zhihu.com
期待大家將自己的問題部落格投稿至上面專欄!
本人能力有限,錯誤之處望大家批評指正!
微積分筆記 無窮級數
研究無窮級數關心的問題 到底能不能收斂成乙個數?本質是研究數列的收斂性 begin sum frac end,lnx的導數就是1 x,所以這個級數是類似於ln n 的,而 begin lim ln x infty end 所以可以知道此時該級數不是乙個數。然後p 1的時候它更大所以肯定也不是乙個數 ...
Matlab 數值微積分與方程求解
專題六 數值微積分與方程求解 目錄 一 數值微分與數值積分 1.數值微分 2.數值積分 二 線性方程組求解 先說差分。差分可以分為向前差分,向後差分和中心差分。三者分別如下表述 matlab提供了求向前差分的函式diff,其呼叫格式有3種 舉個例子 求向量 1,34,54,32,56,78 的一階向...
MATLAB 數值微積分
學習筆記 郭彥甫 yan fu kuo 台大生機系 matlab教學 10數值微積分 representing polynomials in matlab values of polynomials polyval a 9,5,3,7 x 2 0.01 5 f polyval a,x plot x,...