3cib.txt

计算机图形图像设计课程理论

function f=yeyamada(M,q,c,Mf,Q)
% 复制如下数据到C.W.窗口：
% hydraulic_motor([0:20:80],5e-5,5e-12,4,5e-3/6)
% M、Mf：液压马达的负载转矩和摩擦转矩。(单位：N.m)
% q：液压马达排量。(单位：m^3/s)
% c：泄漏系数。(单位：m^3/pa.s)
% Q：输入流量。(单位：m^3/s)
Mt=M+Mf
p=2*pi*Mt/q
deltaQ=c*p
n=60*(Q-deltaQ)/q
etaMv=(Q-deltaQ)/Q
etaMf=M./Mt
etaM=etaMv.*etaMf
x1=linspace(0,84,40);
c1=polyfit(Mt,etaM,3)
c1=polyval(c1,x1);
x2=linspace(10e5,105.5e5,40);
c2=polyfit(p,etaM,3)
c2=polyval(c2,x2);
subplot(2,1,1)
plot(x1,c1)
title('马达效率-扭矩曲线','fontweight','bold','fontsize',10)
xlabel('扭矩MtN.m)','fontweight','bold','fontsize',10)
ylabel('效率\eta','fontweight','bold','fontsize',10)
grid on
hold on
plot(Mt,etaM,'k+')
subplot(2,1,2)
plot(x2,c2)
title('马达效率-工作压力曲线','fontweight','bold','fontsize',10)
xlabel('工作压力paN/m^{2})','fontweight','bold','fontsize',10)
ylabel('效率\eta','fontweight','bold','fontsize',10)
grid on
hold on
plot(p,etaM,'k+')
figure;
S1=spapi(5,Mt,etaM)
S2=spapi(5,p,etaM)
subplot(211)
fnplt(S1)
title('(B样条插值)马达效率-扭矩曲线','fontweight','bold','fontsize',10)
xlabel('扭矩MtN.m)','fontweight','bold','fontsize',10)
ylabel('效率\eta','fontweight','bold','fontsize',10)
hold on
plot(Mt,etaM,'k+')
grid on
subplot(212)
fnplt(S2)
title('(B样条插值)马达效率-工作压力曲线','fontweight','bold','fontsize',10)
xlabel('工作压力paN/m^{2})','fontweight','bold','fontsize',10)
ylabel('效率\eta','fontweight','bold','fontsize',10)
hold on
plot(p,etaM,'k+')
grid on
figure;
x1=4:90;
y1=interp1(Mt,etaM,x1,'spline');
x2=[5.03e5:1e4:106e5];
y2=interp1(p,etaM,x2,'spline')
subplot(211)
plot(Mt,etaM,'k+',x1,y1)
title('(一维样条插值)马达效率-扭矩曲线','fontweight','bold','fontsize',10)
xlabel('扭矩MtN.m)','fontweight','bold','fontsize',10)
ylabel('效率\eta','fontweight','bold','fontsize',10)
grid on
subplot(212)
plot(p,etaM,'k+',x2,y2)
title('(一维样条插值)马达效率-工作压力曲线','fontweight','bold','fontsize',10)
xlabel('工作压力paN/m^{2})','fontweight','bold','fontsize',10)
ylabel('效率\eta','fontweight','bold','fontsize',10)
grid on