分享变系数数值积分问题。

weiniuzhu

四  变系数数值积分问题。
Ti=Ap*R*(p1+p3).*(sin(o)+(R*sin(2*o))./(2*L*(1-(R^2*sin(o).^2)/L^2).^(1/2)))-Ap*R*(p2+p4).*(sin(o)-(R*sin(2*o))./(2*L*(1-(R^2*sin(o).^2)/L^2).^(1/2)));%这个实际上是四缸机的扭矩，o为弧度，此处为自变量
要求Ti在[0,4*pi],上的数值积分，而且每隔一固定角度间隔，要从外部输入p1,p2,p3,p4汽缸压力数据;
function y=calTi(o,p1,p2,p3,p4)
R =52.5e-3;%曲柄半径
L=184e-3;%连杆长度
Ap= 0.0071;%活塞面积
y=Ap*R*(p1+p3).*(sin(o)+(R*sin(2*o))./(2*L*(1-(R^2*sin(o).^2)/L^2).^(1/2)))-Ap*R*(p2+p4).*(sin(o)-(R*sin(2*o))./(2*L*(1-(R^2*sin(o).^2)/L^2).^(1/2)));
%matlab空间代码
clc;clear;
load ptwo.txt;
p1=ptwo(:,1); p2=ptwo(:,2); p3=ptwo(:,3); p4=ptwo(:,4);
tspan=0:4*pi/length(p1):4*pi;
Ti1=[];
for i=1:length(tspan)-1
z=quadl(@calTi,tspan(i),tspan(i+1),[],[],p1(i),p2(i), p3(i), p4(i));
Ti1(i, :)=z ;
end