我是一名MATLAB初学者,希望高手的帮助和指导,谢谢!下面是一段关于分析滑动轴承承载能力的程序:
function[Omax,absvis,ex,temp]=JournalLoad(p)
global Tin SAE n d cd L
temp=Tin;Told=0;
while abs(temp-Told)>0.1
Told=temp;
absvis=AbsViscosity(temp);
ex=Ocvirk(absvis,p);
deltaT=TempRise(temp,ex,p);
temp=Tin+deltaT;
end
absvis=AbsViscosity(temp);
[ex,ocv]=Ocvirk(absvis,p);
Omax=25-ocv;
function absvis=AbsViscosity(tempF)
global SAE
AandB=[9.1209 3.5605;9.1067 3.5386;8.9939 3.4777;8.9133 3.4292;...
8.5194 3.2621;8.3666 3.1884];
SAE10=SAE/10;
absvis=10^(10^(AandB(SAE10,1)-AandB(SAE10,2)*log10(255.2...
+5/9*tempF))-7);
function[ex,ocv]=Ocvirk(vis,p)
global Tin SAR n d cd L
QH=pi*d*(n)*L*ex*cd/2;
phi=atan(pi*sqrt(1-ex^2)/4/ex);
vis=AbsViscosity(temp);
Tr=pi^2*d^3*L*(n)*vis/cd/sqrt(1-ex^2)+0.5*p*ex*cd*sin(phi);
deltaT(0.0640)*Tr*(n)/QH;
global Tin SAE n d cd L
Tin=120;SAE=30;n=3000/60;
d=2;cd=.0026;L=1.1;
options=optimset('display','off');
p=fzero(@JournalLoad,1000,options);
[Omax,absvis,ex,temp]=JournalLoad(p);
disp(['Maximum load='num2str(P)'lb.'])
disp(['Temperature rise='num2str(temp-Tin)'degrees F'])
disp(['Absolute viscosity='num2str(absvis)'reyn'])
disp(['Minimum oil film thickness='num2str(0.5*cd*(1-ex)'inches'])
上面程序在运行的时候怎么提示如下的错误:??? Error: File: E:\段辉-基于Matlab的凸轮设计\JournalLoad.m Line: 37 Column: 22
Missing MATLAB operator.
[本帖最后由 mooni 于 2009-4-16 11:33 编辑]