我们画如下函数的导数图
clc;
clear all
close all
syms x
eps = 0.1;
u = x-x*x+eps*((1.0/(4*pi))*sin(2*pi*x/eps)-(1.0/(2*pi))*x*sin(2*pi*x/eps)-(eps/(4*pi*pi))*cos(2*pi*x/eps)+(eps/(4*pi*pi)));
% 单个点赋值
U0 = double(subs(u,0)); % subs 为符号函数赋值,赋值结果不是double型数据做一下强制转换
U1 = double(subs(u,1.0)); % subs 为符号函数赋值,赋值结果不是double型数据做一下强制转换
diffU = diff(u,x);
X = linspace(0,1,500);
% 多个点赋值(数组赋值)
Udff = double(subs(diffU,X));
figure('name', 'diffU')
plot(X, Udff, 'm-', 'linewidth', 1.5)
xlabel('x', 'Fontsize',15)
hold on
eps = 0.1
eps =0.01