matlab分段 程序,在Matlab中构建分段符号函数

一个选择是使用

heaviside函数使每个方程在给定范围之外等于零,然后将它们全部合并成一个方程：

syms x;

f = (heaviside(x)-heaviside(x-1))*x^3/6 + ...

(heaviside(x-1)-heaviside(x-2))*(1/6)*(-3*x^3+12*x^2-12*x+4) + ...

(heaviside(x-2)-heaviside(x-3))*(1/6)*(3*x^3-24*x^2+60*x-44) + ...

(heaviside(x-3)-heaviside(x-4))*(1/6)*(4-x)^3;

double(int(diff(f, 1)^2, x, 0, 4))

ans =

0.6667

另一个选择是对每个子范围执行每个功能的集成,然后添加结果：

syms x;

eq1 = x^3/6;

eq2 = (1/6)*(-3*x^3+12*x^2-12*x+4);

eq3 = (1/6)*(3*x^3-24*x^2+60*x-44);

eq4 = (1/6)*(4-x)^3;

total = int(diff(eq1, 1)^2, x, 0, 1) + ...

int(diff(eq2, 1)^2, x, 1, 2) + ...

int(diff(eq3, 1)^2, x, 2, 3) + ...

int(diff(eq4, 1)^2, x, 3, 4)

total =

2/3

更新：

尽管在这个问题上提到piecewise功能没有起作用,但Karan’s answer表示,至少在较新的版本中.目前的文献分析表明,它是在R2016b中引入的,但它显然比较早.我在Symbolic Math Toolbox的文档中发现它远远超过R2012b,但是调用语法与现在不同.我在以前的Symbolic Math Toolbox文档中找不到它,但它确实显示为其他工具箱(如统计和样条工具箱)中的一个功能,这解释了它在问题中的提及(为什么它没有在当时为符号方程式工作).