MATLAB复杂隐函数求解!求大神指导!syms D1 D2 G1 G2 al E1 E2 Km w6G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);for D1=1.3:0.1:1.8;for D2=1.3:0.1:1.8;ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)
来源:学生作业帮助网 编辑:作业帮 时间:2024/08/14 05:32:10
![MATLAB复杂隐函数求解!求大神指导!syms D1 D2 G1 G2 al E1 E2 Km w6G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);for D1=1.3:0.1:1.8;for D2=1.3:0.1:1.8;ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)](/uploads/image/z/7788096-0-6.jpg?t=MATLAB%E5%A4%8D%E6%9D%82%E9%9A%90%E5%87%BD%E6%95%B0%E6%B1%82%E8%A7%A3%21%E6%B1%82%E5%A4%A7%E7%A5%9E%E6%8C%87%E5%AF%BC%21syms+D1+D2+G1+G2+al+E1+E2+Km+w6G1%3D10%5E%28-5%29%3BG2%3D10%5E%28-5%29%3B+E1%3D2%2A%2810%5E11%29%3BE2%3D2%2A%2810%5E11%29%3BKm%3D1.2%3Bw6%3D2%2A%2810%5E9%29%3Bfor+D1%3D1.3%3A0.1%3A1.8%3Bfor+D2%3D1.3%3A0.1%3A1.8%3Bac1%3D%28G1%5E2%29%2F%28%28Km%2A0.01%2F2%29%5E%282%2F%28D1-1%29%29%29%3Bac2%3D%28G2%5E2%29%2F%28%28Km%2A0.01%2F2%29%5E%282%2F%28D2-1%29)
MATLAB复杂隐函数求解!求大神指导!syms D1 D2 G1 G2 al E1 E2 Km w6G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);for D1=1.3:0.1:1.8;for D2=1.3:0.1:1.8;ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)
MATLAB复杂隐函数求解!求大神指导!
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);for D1=1.3:0.1:1.8;
for D2=1.3:0.1:1.8;
ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));
ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)));
w=2*(pi^(1/2))*E1*(G1^(D1-1))*D1*(al^(D1/2))/3;
r=(2/(3-2D1))*(al^((3-2D1)/2)-ac1^((3-2D1)/2))-((al^((2-D2)/2))*2/(5-2D1-D2)*(al^((5-2D1-D2)/2)-ac1^((5-2D1-D2)/2)));
t=Km*w6*D1*(al^(D1/2))/2;
u=(2/(2-D1)*(ac1^((2-D1)/2))-((al^((2-D2)/2)))*2/(4-D2-D1)*(ac1^((4-D2-D1)/2)));
P2=w*r+t*u;
w1=2*(pi^(1/2))*E2*(G2^(D2-1))*D1*(al^((D1+D2-2)/2))/3;
r1=2/(5-2D1-D2)*(al^((5-2D1-D2)/2)-ac2^((5-2D1-D2)/2));
t1=Km*w6*D1*(al^((D1+D2-2)/2))/2;
u1=2/(4-D2-D1)*(ac2^((4-D2-D1)/2));
P1=w1*r1+t1*u1;
P=P1+P2;P=3*10^4;
%想在已知P和D1,D2的情况下求解al,D1和D2是一组数据,最后画出三维图
MATLAB复杂隐函数求解!求大神指导!syms D1 D2 G1 G2 al E1 E2 Km w6G1=10^(-5);G2=10^(-5); E1=2*(10^11);E2=2*(10^11);Km=1.2;w6=2*(10^9);for D1=1.3:0.1:1.8;for D2=1.3:0.1:1.8;ac1=(G1^2)/((Km*0.01/2)^(2/(D1-1)));ac2=(G2^2)/((Km*0.01/2)^(2/(D2-1)
题目写的比较乱.
明确几点:
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范围内有两个解.