返回信息流
这是一条镜像帖。来源：北邮人论坛 / matlab / #12869同步于 2019/6/3
该镜像源已超过 30 天没有更新，可能在源站已被删除。
Matlab机器人发帖
求救 红线内的怎么用matlab实现啊
zhangxiaomei
2019/6/3镜像同步3 回复
订阅后，新回复会通过你的通知中心匿名送达。
3 条回复
只是复合辛普森的话, 可以参考这个
syms x
func1 = x^3; % f(x)
tmp = matlabFunction(diff(func1,x)); % f(x)求导并转化为function类型
L_func = @(x)sqrt(1+(tmp(x)).^2); % 弧长公式
% 输出结果
disp(compositeSimpsonRule(L_func, 800, 0, 1))
%%
% func 函数名称
% M 网格划分大小
% 积分下限
% 积分上限
function I = compositeSimpsonRule(func, M, left, right)
x = linspace(left, right, M);
I = 0;
h = (right-left)/M;
for i = 1:M-1
I = I+h/6*(func(x(i))+4*func((x(i)+x(i+1))/2)+func(x(i+1)));
end
end
【 在 dhown 的大作中提到: 】
: 只是复合辛普森的话, 可以参考这个
: [code=py]
: syms x
: ...................
谢谢啊 我去试一试