Matlab 多項式及數值微積分

differentiation微分

兩個多項式相乘時,兩個多項式用兩個矩陣接收,利用conv卷積函式進行合併

integration積分數值微分

g = colormap(lines); hold on;

for i=1:4
x = 0:power(10, -i):pi;
y = sin(x); m = diff(y)./diff(x);
plot(x(1:end-1), m, 'color', g(i,:));
endhold off;
set(gca, 'xlim', [0, pi/2]); set(gca, 'ylim', [0, 1.2]);
set(gca, 'fontsize', 18); set(gca, 'fontname', 'symbol');
set(gca, 'xtick', 0:pi/4:pi/2);
set(gca, 'xticklabel', );
h = legend('h=0.1','h=0.01','h=0.001','h=0.0001');
set(h,'fontname', 'times new roman'); box on;

數值積分

柱形積分

梯形積分

辛普森simpson規則積分

自建函式繪製函式影象——函式指標(控制代碼)

寫乙個函式，傳遞函式指標和定義域作為引數，可以直接繪製函式指標指向的函式的影象。

function [y] = xy_plot(input,x)

% xy_plot receives the handle of a function
% and plots that function of x
y = input(x); plot(x,y,'r--');
xlabel('x'); ylabel('function(x)');
end

y = @(x) 1./(x.^3-2*x-5);      %一重積分求面積

integral(y,0,2) %傳參:函式指標,下限，上限

f = @(x,y) y.*sin(x)+x.*cos(y);   %二重積分求體積

integral2(f,pi,2*pi,0,pi) %傳參: 函式指標,下限，上限，下限，上限

f = @(x,y,z) y.*sin(x)+z.*cos(y);

integral3(f,0,pi,0,1,-1,1)

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,...

MATLAB多項式及多項式擬合

多項式均表示為陣列形式，陣列元素為多項式降冪係數 1.polyval函式 求多項式在某一點或某幾個點的值.p 1,1,1 x 2 x 1 x 1,0,1 y polyval p,x 另外求函式在某一點或某幾個點的值可以用函式feval.x 1,0,1 y feval x exp x x 注意用的乘法...

MATLAB數值法與微積分

函式之微分為求函式對自變數之導數，或為其斜率 利用數值方法則可以解出其他相關之問題，其應用部份已在前章討論。數值微分有兩種應用，其一是在資料收集完備後，分析其變化速度 其二為即時估計或量測速率。後者需要快速演演算法才能有立即反應。計算斜率，依其定義即為dy dx，在數值分析上必須轉化為可量測之變化量...