Question

用MATLAB 求曲线 f（x）=sin（x）+x^2 的曲率。 求编程

Answer 1

计算模型
f(x)=sin(x)+x^2曲率计算
计算方法
利用符号运算
算法流程
（1）用符号表示表示f(x)
（2）求f(x)的一阶导数，二阶导数
（3）根据曲率半径定义r=|(1+y'^2)^(3/2)/y''|
计算代码。
clc;clear all
syms x
f=sin(x)+x^2;
f1=diff(f,x);
f2=diff(f1,x);
r=(1+f1^2)^(3/2)/f2;
xx=0:0.1:30;rr=zeros(size(xx));

for i=1:length(xx)
rr(i)=abs(subs(r,x,xx(i)));
end
plot(xx,rr)
xlabel('x')
ylabel('曲率半径r')