離散格式的任意階導數,只要精度和點數確定後。對應各點的係數可以通過taylor展開待定係數法計算得到。個人構造一些格式時,常常手動展開,計算係數的程式也已經寫過很多次。這裡乾脆掛在網上,避免每次使用的時候都得找半天甚至重寫。
當前格式我們只實現了最一般情形下的離散係數計算,緊緻格式係數的計算和重構係數的計算這裡均不考慮。另外,也不考慮給定超過精度所需的點,然後進行譜性質優化的情形。
程式構造考慮任意階精度,構造任意階導數的離散逼近,以及可以使用任意的點。
%計算任意階導數的離散係數
clc;
clear;
format rat;
%給定計算使用的點,其中0代表我們所需要構造導數點的位置
point = [0 1 2 3];
n = length(point);
%要模擬幾階導數
order = 1;
%計算taylor展開的係數
a = zeros(n);
for i = 1:n
for j = 1:n
a(i,j) = (point(j))^(i-1);
end
end
b = zeros(n,1);
b(order+1) = 1;
%求解係數並輸出
coefficient = a\b;
%輸出[point; coefficient']
當前是使用i,i+1,i+2,i+3
一共四個點離散逼近i點處的一階導數,輸出結果如下:
ans =
0 1 2 3
-11/6 3 -3/2 1/3
Matlab中數值計算精度
1.matlab中有三種運算精度,它們分別為數值演算法 符號演算法和可控精度演算法,將它們分別介紹如下 1 數值演算法把每個數取為16位,計算按浮點運算進行,它是運算速度最快的一種演算法.2 符號演算法把每個數都變為符號量,運算按有理量計算進行,它的優點是能夠得到精確結果,缺點是占用空間大,並且運算...
標題用C去計算行列式(任意階)
c cpp include include include struct node 請先輸入幾位的行列式 交換兩個資料 void swap int a,int b 通過這種方式去進行交換兩個數是無法達到目的的 計算n階行列式中的每乙個行列式的值,sum就是計算行列式的最終結果 void print ...
matlab計算任意多邊形面積
已知多邊形的各個點座標,且點按瞬時間或逆時間順序排列,求多邊形面積。首先根據計算多邊形面積的計算公式,用向量法計算比較方便。兩向量叉乘為這倆向量組成的三角形面積的兩倍,所以多邊形面積公式為 s 0.5 oa 叉乘 ob ob 叉乘 oc o 叉乘 oa matlab 實現 第一步 取得座標點資料,讀...