matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/05 10:55:58
matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(

matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(
matlab solve
global b
global n
global y
t1=0.1;
t2=0.001;
t3=0.005;
t4=0.01;
t5=0.05;
t6=0.1;
t7=0.002;
t8=0.006;
t9=0.07;
t10=0.1;
shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(t10-y)^(b-1))-(b-1)*((t1-y)^(-1)+(t2-y)^(-1)+(t3-y)^(-1)+(t4-y)^(-1)+(t5-y)^(-1)+(t6-y)^(-1)+(t7-y)^(-1)+(t8-y)^(-1)+(t9-y)^(-1)+(t10-y)^(-1))=0';
shizi2='-j/n+1/n^2*((t1-y)^b+(t2-y)^b+(t3-y)^b+(t4-y)^b+(t5-y)^b+(t6-y)^b+(t7-y)^b+(t8-y)^b+(t9-y)^b+(t10-y)^b)=0';
shizi3='n/b+(log(t1-y)+log(t2-y)+log(t3-y)+log(t4-y)+log(t5-y)+log(t6-y)+log(t7-y)+log(t8-y)+log(t9-y)+log(t10-y))-1/n*((t1-y)^b*log(t1-y)+(t2-y)^b*log(t2-y)+(t3-y)^b*log(t3-y)+(t4-y)^b*log(t4-y)+(t5-y)^b*log(t5-y)+(t6-y)^b*log(t6-y)+(t7-y)^b*log(t7-y)+(t8-y)^b*log(t8-y)+(t9-y)^b*log(t9-y)+(t10-y)^b*log(t10-y))=0';
[b,y,n]=solve(shizi1,shizi2,shizi3,b,y,n);
哪儿错了?
错误提示:
Error using mupadengine/feval (line 157)
MuPAD error:Error:']' expected.[line 1,col 35]
Error in solve (line 160)
sol = eng.feval('symobj::solvefull',eqns,vars);
Error in disanti (line 17)
[a,d,f]=solve(shizi1,shizi2,shizi3,b,y,

matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+(
将前3行全局变量定义global b;global n;global y;
改为符号变量定义:syms b n y.
希望你的shizi表达式有意义.我的机子求解你的问题一直busy

没有报错啊,结果是:
ans =

-((COH0^2*COHadd^2*Kb^4 + 2*COH0^2*COHadd*Kb^3 + 4*COH0^2*Cblg*Kb^3*N1 + 4*COH0^2*Cblg*Kb^3*N2 + 4*COH0^2*Cblg*Kb^3*N3 + 4*COH0^2*Cblg*Kb^3*N4 + 4*COH0^2*Cblg*Kb^3*N5 +...

全部展开

没有报错啊,结果是:
ans =

-((COH0^2*COHadd^2*Kb^4 + 2*COH0^2*COHadd*Kb^3 + 4*COH0^2*Cblg*Kb^3*N1 + 4*COH0^2*Cblg*Kb^3*N2 + 4*COH0^2*Cblg*Kb^3*N3 + 4*COH0^2*Cblg*Kb^3*N4 + 4*COH0^2*Cblg*Kb^3*N5 + 4*COH0^2*Cblg*Kb^3*N6 + 4*COH0^2*Cblg*Kb^3*N7 + COH0^2*Kb^2 + 2*COH0*COHadd^2*Kb^3 - 2*COH0*COHadd*Cblg*Kb^3*N1 - 2*COH0*COHadd*Cblg*Kb^3*N2 - 2*COH0*COHadd*Cblg*Kb^3*N3 - 2*COH0*COHadd*Cblg*Kb^3*N4 - 2*COH0*COHadd*Cblg*Kb^3*N5 - 2*COH0*COHadd*Cblg*Kb^3*N6 - 2*COH0*COHadd*Cblg*Kb^3*N7 + 4*COH0*COHadd*Kb^2 + 6*COH0*Cblg*Kb^2*N1 + 6*COH0*Cblg*Kb^2*N2 + 6*COH0*Cblg*Kb^2*N3 + 6*COH0*Cblg*Kb^2*N4 + 6*COH0*Cblg*Kb^2*N5 + 6*COH0*Cblg*Kb^2*N6 + 6*COH0*Cblg*Kb^2*N7 + 2*COH0*Kb + COHadd^2*Kb^2 - 2*COHadd*Cblg*Kb^2*N1 - 2*COHadd*Cblg*Kb^2*N2 - 2*COHadd*Cblg*Kb^2*N3 - 2*COHadd*Cblg*Kb^2*N4 - 2*COHadd*Cblg*Kb^2*N5 - 2*COHadd*Cblg*Kb^2*N6 - 2*COHadd*Cblg*Kb^2*N7 + 2*COHadd*Kb + Cblg^2*Kb^2*N1^2 + 2*Cblg^2*Kb^2*N1*N2 + 2*Cblg^2*Kb^2*N1*N3 + 2*Cblg^2*Kb^2*N1*N4 + 2*Cblg^2*Kb^2*N1*N5 + 2*Cblg^2*Kb^2*N1*N6 + 2*Cblg^2*Kb^2*N1*N7 + Cblg^2*Kb^2*N2^2 + 2*Cblg^2*Kb^2*N2*N3 + 2*Cblg^2*Kb^2*N2*N4 + 2*Cblg^2*Kb^2*N2*N5 + 2*Cblg^2*Kb^2*N2*N6 + 2*Cblg^2*Kb^2*N2*N7 + Cblg^2*Kb^2*N3^2 + 2*Cblg^2*Kb^2*N3*N4 + 2*Cblg^2*Kb^2*N3*N5 + 2*Cblg^2*Kb^2*N3*N6 + 2*Cblg^2*Kb^2*N3*N7 + Cblg^2*Kb^2*N4^2 + 2*Cblg^2*Kb^2*N4*N5 + 2*Cblg^2*Kb^2*N4*N6 + 2*Cblg^2*Kb^2*N4*N7 + Cblg^2*Kb^2*N5^2 + 2*Cblg^2*Kb^2*N5*N6 + 2*Cblg^2*Kb^2*N5*N7 + Cblg^2*Kb^2*N6^2 + 2*Cblg^2*Kb^2*N6*N7 + Cblg^2*Kb^2*N7^2 + 2*Cblg*Kb*N1 + 2*Cblg*Kb*N2 + 2*Cblg*Kb*N3 + 2*Cblg*Kb*N4 + 2*Cblg*Kb*N5 + 2*Cblg*Kb*N6 + 2*Cblg*Kb*N7 + 1)^(1/2) + COH0*Kb - COHadd*Kb + Cblg*Kb*N1 + Cblg*Kb*N2 + Cblg*Kb*N3 + Cblg*Kb*N4 + Cblg*Kb*N5 + Cblg*Kb*N6 + Cblg*Kb*N7 - COH0*COHadd*Kb^2 + 1)/(2*(COH0*Kb^2 + Kb))
-(COH0*Kb - (COH0^2*COHadd^2*Kb^4 + 2*COH0^2*COHadd*Kb^3 + 4*COH0^2*Cblg*Kb^3*N1 + 4*COH0^2*Cblg*Kb^3*N2 + 4*COH0^2*Cblg*Kb^3*N3 + 4*COH0^2*Cblg*Kb^3*N4 + 4*COH0^2*Cblg*Kb^3*N5 + 4*COH0^2*Cblg*Kb^3*N6 + 4*COH0^2*Cblg*Kb^3*N7 + COH0^2*Kb^2 + 2*COH0*COHadd^2*Kb^3 - 2*COH0*COHadd*Cblg*Kb^3*N1 - 2*COH0*COHadd*Cblg*Kb^3*N2 - 2*COH0*COHadd*Cblg*Kb^3*N3 - 2*COH0*COHadd*Cblg*Kb^3*N4 - 2*COH0*COHadd*Cblg*Kb^3*N5 - 2*COH0*COHadd*Cblg*Kb^3*N6 - 2*COH0*COHadd*Cblg*Kb^3*N7 + 4*COH0*COHadd*Kb^2 + 6*COH0*Cblg*Kb^2*N1 + 6*COH0*Cblg*Kb^2*N2 + 6*COH0*Cblg*Kb^2*N3 + 6*COH0*Cblg*Kb^2*N4 + 6*COH0*Cblg*Kb^2*N5 + 6*COH0*Cblg*Kb^2*N6 + 6*COH0*Cblg*Kb^2*N7 + 2*COH0*Kb + COHadd^2*Kb^2 - 2*COHadd*Cblg*Kb^2*N1 - 2*COHadd*Cblg*Kb^2*N2 - 2*COHadd*Cblg*Kb^2*N3 - 2*COHadd*Cblg*Kb^2*N4 - 2*COHadd*Cblg*Kb^2*N5 - 2*COHadd*Cblg*Kb^2*N6 - 2*COHadd*Cblg*Kb^2*N7 + 2*COHadd*Kb + Cblg^2*Kb^2*N1^2 + 2*Cblg^2*Kb^2*N1*N2 + 2*Cblg^2*Kb^2*N1*N3 + 2*Cblg^2*Kb^2*N1*N4 + 2*Cblg^2*Kb^2*N1*N5 + 2*Cblg^2*Kb^2*N1*N6 + 2*Cblg^2*Kb^2*N1*N7 + Cblg^2*Kb^2*N2^2 + 2*Cblg^2*Kb^2*N2*N3 + 2*Cblg^2*Kb^2*N2*N4 + 2*Cblg^2*Kb^2*N2*N5 + 2*Cblg^2*Kb^2*N2*N6 + 2*Cblg^2*Kb^2*N2*N7 + Cblg^2*Kb^2*N3^2 + 2*Cblg^2*Kb^2*N3*N4 + 2*Cblg^2*Kb^2*N3*N5 + 2*Cblg^2*Kb^2*N3*N6 + 2*Cblg^2*Kb^2*N3*N7 + Cblg^2*Kb^2*N4^2 + 2*Cblg^2*Kb^2*N4*N5 + 2*Cblg^2*Kb^2*N4*N6 + 2*Cblg^2*Kb^2*N4*N7 + Cblg^2*Kb^2*N5^2 + 2*Cblg^2*Kb^2*N5*N6 + 2*Cblg^2*Kb^2*N5*N7 + Cblg^2*Kb^2*N6^2 + 2*Cblg^2*Kb^2*N6*N7 + Cblg^2*Kb^2*N7^2 + 2*Cblg*Kb*N1 + 2*Cblg*Kb*N2 + 2*Cblg*Kb*N3 + 2*Cblg*Kb*N4 + 2*Cblg*Kb*N5 + 2*Cblg*Kb*N6 + 2*Cblg*Kb*N7 + 1)^(1/2) - COHadd*Kb + Cblg*Kb*N1 + Cblg*Kb*N2 + Cblg*Kb*N3 + Cblg*Kb*N4 + Cblg*Kb*N5 + Cblg*Kb*N6 + Cblg*Kb*N7 - COH0*COHadd*Kb^2 + 1)/(2*(COH0*Kb^2 + Kb))
接着用subs函数代入数据就好了
望采纳,谢谢!

收起

MATLAB,solve函数, matlab global定义的变量是什么类型 matlab怎么用solve求解析解 为什么solve函数后matlab无反应 global global 写一篇英语短文 how to solve the problem of global warming Matlab中solve解非线性方程组用的是什么算法?能否查看solve的源程序? matlab 中的 solve怎么用 为什么我老是出错 matlab解方程组,通常用solve还是用矩阵? matlab中solve函数不能用的原因 matlab解非线性方程组解不出来:代码如下:global b T temp;b=1.2;T=15;temp=1;syms a R;[a,R]=solve('a-(R-b)/2+0.5*((R-b)^2-0.281*R*T)^0.5=0','13*temp-13*(a+b)*a-0.912*T*R=0')运行结果如下:a =-(0.25*(0.281*R*T + (R - 2.0*a)^2 - 1.0 matlab solve global bglobal nglobal yt1=0.1;t2=0.001;t3=0.005;t4=0.01;t5=0.05;t6=0.1;t7=0.002;t8=0.006;t9=0.07;t10=0.1;shizi1='b/n*((t1-y)^(b-1)+(t2-y)^(b-1)+(t3-y)^(b-1)+(t4-y)^(b-1)+(t5-y)^(b-1)+(t6-y)^(b-1)+(t7-y)^(b-1)+(t8-y)^(b-1)+(t9-y)^(b-1)+( matlab solve函数 xmaxr=solve(dydx,x)其中dydx是一个x和r的多项式 matlab中solve函数solve函数最多可以解多少维的方程组,上千维的可以吗, matlab中,为什么solve不能解三解方程?比如:solve('cos(3a)=1','a') 有没有其它的解法? 用MATLAB如何解指数方程matlab中用solve解 下面的方程解不了,solve('x^0.65-0.1875*x^-0.35+1.991-1.991*x=0.8125','x') matlab solve 问题clear all;clc;solve('a+5=0','a');然后matlab 显示ans 为,没有具体解.不知道这是为什么,难道是我solve 理解有问题?