MATLAB复杂隐函数求解!求大神指导!

发布网友发布时间：2022-05-16 00:59

我来回答

共1个回答

热心网友时间：2023-10-06 03:46

题目写的比较乱。

明确几点：

1、是不是说对于D1=1.3:0.1:1.8，D2=1.3:0.1:1.8的每一种组合，按照P=P1+P2=3*10^4来求解a1？

2、对于这样的复杂方程，求解析解是不太现实的，应该求数值解。

3、方程可能存在多个解，怎么处理？例如，取D1=D2=1.5，画出P-3E4的曲线如下：

显然，在0-30范围内有两个解。

追问你好，谢谢你的回答。以上函数中的al在我的设定中是有一定物理意义的，应该是个小于1大于0的正数。是否能够给出：在已知P的值时，得出D1，D2，al的三维关系图。麻烦大神把程序传上来！我快被这事弄崩溃了！

追答

参考代码：

syms D1 D2 G1 G2 al E1 E2 Km w6
G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);
[DD1,DD2] = meshgrid(1.3:0.1:1.8);
A1 = DD1*NaN;
for ii=1:size(DD1,1)
    for jj=1:size(DD2,2)
        D1 = DD1(ii,jj);
        D2 = DD2(ii,jj);
        ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));

        （此处有部分代码无改动，因字数*省略）
        u1=2/(4-D2-D1)*(ac2^((4-D2-D1)/2));
        P1=w1*r1+t1*u1;
        P=P1+P2;
        try
            A1(ii,jj)=fzero(@(a1)subs(P-3E4,a1),[eps 1]);
        catch
            fprintf('Error: D1=%g, D2=%g, %s\n',D1,D2,lasterr)
        end
    end
end
mesh(DD1,DD2,A1)
view(3)

注意，有些情况下（例如D1=1.3, D2=1.4）方程无解。