matlab之一元多项式、数值微分、数值积分

一元多项式

1.绘图,求值

$$f(x)=9x^3-5x^2+3x+7(-2\le x\le 5),f(x)=?$$

a = [9 -5 3 7];%代表多项式各项系数
x = -2:0.01:5;
f = polyval(a,x);%求值
%美化部分
plot(x,f,'LineWidth', 2);
xlabel('x'); ylabel('f(x)');
set(gca, 'FontSize', 14)

2.求导

$$f(x)=5x^4-2x^2+1,f^{\prime }(x)=?$$

p=[5 0 -2 0 1];
polyder(p);

3.多项式相乘

$$f(x)=(5x^3-7x^2+5x+10)(4x^2+12x-3)=?$$

a=[5 -7 5 10];
b=[4 12 -3];
f=conv(a,b);

4.求不定积分

$$已知f(x)=5x^4-2x^2+1,求\int f(x)$$

p=[5 0 -2 0 1];
f=polyint(p, 3);%第二个参数代表反导后的常数项k

数值微分

1.基本原理

$$f^{\prime }(x_0)=\underset{h\to 0}{\lim }\frac{f(x_0+h)-f(x_0)}{h}$$

2.实现方法

$$计算不同h下，用数值微分求出的f(x)=sin(x)导数的函数$$

g = colormap(lines);%
hold on;
for i=1:4
    x = 0:power(10, -i):pi;
    y = sin(x);
    m = diff(y)./diff(x);%diff函数求出数组内相邻数组的差值
    plot(x(1:end-1), m, 'Color', g(i,:)); 
end
hold 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', {'0', 'p/4', 'p/2'});
h = legend('h=0.1','h=0.01','h=0.001','h=0.0001');
set(h,'FontName', 'Times New Roman'); 
box on;

3.扩展

$$知f(x)=x^3,求-2\le x\le 2范围内的f^{\prime },f^{\prime \prime }$$

x = -2:0.005:2;
y = x.^3;
m = diff(y)./diff(x);
m2 = diff(m)./diff(x(1:end-1));
plot(x,y,x(1:end-1),m,x(1:end-2),m2);
%美化部分
xlabel('x', 'FontSize', 18); 
ylabel('y', 'FontSize', 18);
legend('f(x) = x^3','f''(x)','f''''(x)', 4);
set(gca, 'FontSize', 18);

数值积分

根据积分的性质，可以在MATLAB上编写，但是已经有内置函数了，就只放图

1.一重积分

$${\int }_0^2\frac{1}{x^3-2x-5}dx$$

y = @(x) 1./(x.^3-2*x-5);
ans=integral(y,0,2);

2.二重积分

$$f(x,y)={\int }_0^{\pi }{\int }_{\pi }^{2\pi }(y\cdot sin(x)+x\cdot cos(y))dxdy$$

f = @(x,y) y.*sin(x)+x.*cos(y);
ans=integral2(f,pi,2*pi,0,pi);

3.三重积分

$${\int }_{-1}^1{\int }_0^1{\int }_0^{\pi }(y\cdot sin(x)+z\cdot cos(y))dxdydz$$

f = @(x,y,z) y.*sin(x)+z.*cos(y);
ans=integral3(f,0,pi,0,1,-1,1)