哪位大侠能帮我看看下面matlab的代码,非常感谢
发布网友
发布时间：2023-01-04 21:53
共1个回答
热心网友
时间：2023-10-14 05:23
time=0:timeg:n*timeg;%积分区间
Ssp(1)=0; %积分值的初值
for i=1:(n-1)
xx=[time(i),time(i+1)]; %依次选取两个点
ff=[integfunction(i),integfunction(i+1)];%选取点所对应的积分函数值
pp=spline(xx,ff);%三次样条差值
int_pp=fnint(pp);%得出样条函数的积分表达式
Ssp(i+1)=ppval(int_pp,[time(i),time(i+1)])*[-1;1];%得到此处积分(相当于无限细分的矩形)值并存入Ssp
end
for i=1:(n-1)
Ssp(i+1)=Ssp(i)+Ssp(i+1);%所有积分值加起来得到总的积分值（所有小矩形面积相加）
end
Strain=2*C/B*Ssp;%得到表达式的值
把第一句的A换成A^2即可
