題目:考察具有粘性阻尼的單擺
首先編寫一段 maple 子程式 pendphaseplot 如下:
pendphaseplot := proc (c, k, ui, uf, n1, vi, vf, n2, ti, tf, n, vu, vv)
local s, i, a, j, b, sys, fcns, q, p; s := {};
定義變數
for i from 0 to n1 do a := ui+i*(uf-ui)/n1;
a是u變數變化的步長
for j from 0 to n2 do b := vi+j*(vf-vi)/n2;
b是v變數變化的步長
sys := diff(u(t), t) = v(t), diff(v(t), t) = -k*sin(u(t))-c*v(t);
定義系統
fcns := ;
定義方程(帶入dsolve中用的,就是個過程變數)
q := dsolve(, fcns, type = numeric, method = rkf45);
帶入不同的初值求得相軌跡
p := plots[odeplot](q, [u(t), v(t)], ti .. tf, numpoints = n, view = [vu, vv], scaling = constrained);
將上一步求得的相軌跡作出相應的影象(通過數值解作圖)
s := s union p
union是取並集的操作,其實就是用s儲存p
end do
end do;
plots[display]([op(s)])
op的含義就是提取整個s用於作圖,從這裡感覺maple程式設計是乙個比較笨重的語言,當然了,用op也可以提取固定的一些項,詳細看幫助即可
end proc
至此整個程式結束
命令列:
無阻尼情況:
有阻尼的情況:
線性代數Ax b原理及工程上的應用
前言 線性代數在工程應用上十分廣泛,在座標系轉換,深度學習,求解演算法的優化解方面有著大量應用。因此掌握線性代數的基本理論,並且具有解決實際工程問題的能力尤為重要。線性方程組解的情況 線性方程組的解的三種情況 1.適定方程組 存在唯一解 2.欠定方程組 存在多解。變數數 方程組數 3.超定方程組 無...
c MVC 在工程中應用的實際處理
1.把collection放到另外乙個專案中。新增一類別庫 classlibrary1,在此類別庫中,定義所有的collection 2.在webmvc的專案中引用此類別庫。為了正常使用,需要在此專案中的global.asax檔案中,做如下定義。public static void register...
線性光耦在電流取樣中的應用
1 引言 在現代電氣測量和控制中,常常需要用低壓器件去測量 控制高電壓 強電流等模擬量,如果模擬量與數字量之間沒有電氣隔離,那麼,高電壓 強電流很容易串入低壓器件,並將其燒毀。線性光耦 hcnr200可以較好地實現模擬量與數字量之間的隔離,隔離電壓峰值達8000v 輸出跟隨輸入變化,線性度達0.01...