matlab高手进，隐函数求数值解，问题很急，需要编程！

Question

syms b cp;
u=2.65e+10;
c1b=5.057e+10;
cv=2596.15;
cs=3128.89;
us01=10;
us12=20;
h=1e-9;
cp=2400;
f=tan(b*sqrt((cp/cv)^2-1))-(c1b*sqrt((cp/cv)^2-1)*(u*sqrt(1-(cp/cs)^2)+(b/h)*us12-(b/h)*us01))/((u*(b/h)*sqrt(1-(cp/cs)^2)+((b/h)^2)*us01)*us01+c1b^2*((cp/cv)^2-1));
我的未知数是b和cp, 现在要求b从0到20时，每一个b对应的cp的值（b可以取0,1,2……20）。
不要告诉我用solve或fsolve，我用了根本求不出解来，cp是多解的，而且是数值解。
我要求解程序，求出的解应该是矩阵形式，高分悬赏，越快越好，谢了各位

Answer 1

你的公式中 cp的取值范围是cv到cs之间

用fzero在这个区间之内找结果就好了

u=2.65e+10;
c1b=5.057e+10;
cv=2596.15;
cs=3128.89;
us01=10;
us12=20;
h=1e-9;
b=(0:20)';
A=cell(21,2);
f=@(b,cp) tan(b*sqrt((cp/cv)^2-1))-(c1b*sqrt((cp/cv)^2-1)*(u*sqrt(1-(cp/cs)^2)+(b/h)*us12-(b/h)*us01))/((u*(b/h)*sqrt(1-(cp/cs)^2)+((b/h)^2)*us01)*us01+c1b^2*((cp/cv)^2-1));
for n=1:length(b)
    A{n,1}=b(n);
    for p0=cv+0.01:cs-0.01
        [s v flag]=fzero(@(cp) f(b(n),cp),p0);
        if flag==1
            if isempty(A{n})
                A{n,2}=s;
            elseif all(abs(A{n,2}-s)>1e-3)
                A{n,2}=[A{n,2},s];
            end
        end 
    end        
end

由于每个b的取值得到cp的解的个数不同，所以这里用cell矩阵来储存结果

最后结果再A中

A的第一列是b的取值

第二列是cp的可能取值，如果结果是[]，也就是无解

结果如下，第一列是b的取值，第二列是对应的cp的可能取值

追问

先谢谢你的回答，我用matlab画过此方程的曲线(matlab不需要数据准备可直接绘图)，常理上cp应该在cv和cs之间，可是图上表明在0到cv时cp也有值，这也是我要求数据的原因，因此我需要cp在0到cv上的解这才是关键，还有我需要k密一些，完成这些我就立即采纳绝不拖欠，谢谢啦！

追答

u=2.65e+10;
c1b=5.057e+10;
cv=2596.15;
cs=3128.89;
us01=10;
us12=20;
h=1e-9;
b=(0:20)';
f=@(b,cp) tan(b.*sqrt((cp/cv).^2-1))-(c1b*sqrt((cp/cv).^2-1).*(u*sqrt(1-(cp/cs).^2)+(b/h)*us12-(b/h)*us01))./((u*(b/h).*sqrt(1-(cp/cs).^2)+((b/h).^2)*us01)*us01+c1b^2*((cp/cv).^2-1));
ezplot(f,[0 20 0 3000]);

没看到小于cv的解

追问

我和你画的不一样啊？这是为什么呢？我的程序是：

公式太长没截完，我是先把f代入命令窗口把已知数都消除，就剩b和cp之后的式子整体代入ezplot，画出的图是：

我觉得不应该错啊，你试试？我老师就是觉得这个现象奇怪，才让我求数值解的

