Matlab程式設計計算任意階精度的離散格式係數

2021-09-27 13:09:55 字數 886 閱讀 7155

離散格式的任意階導數,只要精度和點數確定後。對應各點的係數可以通過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 實現 第一步 取得座標點資料,讀...