function s=Ct(v)
if v<11.4
    error
    return
end
eq1='2*1.29*pi*63^2*v^3*x*(1-x)^2-5e6';
x=solve(eq1);
x=eval(x);
x=real(x);
a=x(x>0&x<0.5);
s=4*a*(1-a);
return

15楼: Originally posted by sunshine819 at 2013-10-19 17:04:47
大概的意思就是这样，不过你在simulink模型里最好加一个初始化文件，这个文件里面包含你计算出来的数据，这样就不用每次打开的时候都重新生成一遍这两个数据了...

ZF=320-70;
w=698.094;
L=902.2;
EA=384243000;
CB=1;   %静摩擦系数
x=[4.9144e5,5.4374e5];   %fsolve所用初值 (XF=828.67   ZF=250)  
XF=853.87-5.2-20:0.1:853.87-5.2+20;
n=length(XF);
for i=1:n  %for循环 改变方程中XF的值
    if L-x(2)/w>0    %有躺底
    f1=@(x)([L-x(2)/w+(x(1)/w)*log(x(2)/x(1)+sqrt(1+(x(2)/x(1))^2))+x(1)*L/EA+CB*w*(-(L-x(2)/w)^2+(L-x(2)/w-x(1)/(CB*w))*max(L-x(2)/w-x(1)/(CB*w),0))/(2*EA)-XF(i);
        (x(1)/w)*(sqrt(1+(x(2)/x(1))^2)-1)+x(2)^2/(2*w*EA)-ZF]);
    z=fsolve(f1,x,optimset('Display','off'));
    HF(i)=z(1);
    VF(i)=z(2);
    XF(i)
    x=z     %将方程的解作为下一步fsolve初值
else   %无躺底
    f2=@(x)([(x(1)/w)*(log(x(2)/x(1)+sqrt(1+(x(2)/x(1))^2))-log((x(2)-w*L)/x(1)+sqrt(1+((x(2)-w*L)/x(1))^2)))+x(1)*L/EA-XF(i);
        (x(1)/w)*(sqrt(1+(x(2)/x(1))^2)-sqrt(1+((x(2)-w*L)/x(1))^2))+(x(2)*L-w*L^2/2)/EA-ZF]);
    z=fsolve(f2,x,optimset('Display','off'));
    HF(i)=z(1);
    VF(i)=z(2);
    XF(i)
    x=z     %将方程的解作为下一步fsolve初值
    end
end

4楼: Originally posted by sp_tried at 2013-10-19 13:04:16
不知我理解的是否正确：你似乎是想通过那段程序，用v值求出s，把计算得到的结果用于仿真，根据不同的v值查表得到s值。试试这样做：
1. 把你的v和算出的s保存为work space中的行向量；
2. 在lookup table参数设置页 ...