zx199910 2021-01-30 23:44
浏览 824

matlab求解含三个未知数的三个方程组成的方程组求不出正确解

s=tf('s')
Gop1=1.3/(1+s*2)*exp(-1.8*s)
GCL1=Gop1/(1+Gop1);
y1=step(GCL1,[0:0.01:50]);
figure(1)
step(GCL1,[0:0.01:50])


%从图中可以读出: 
%pusai=(0.846-0.621)/(0.846-0.565)=0.8007;
%w=2*3.14/6.3;
%m=0.2567;

syms x y z real
w=2*3.14/6.3;
m=0.2567;
s=-0.2567*w+j*w;
Gs=x/(1+s*y)*exp(-z*s)
Im=imag(Gs)
Re=real(Gs)
eqns=[x/(1+x)==0.565,(x*sin((314*z)/315)*exp((4609615792312965*z)/18014398509481984)*((4609615792312965*y)/18014398509481984 - 1))/(((4609615792312965*y)/18014398509481984 - 1)^2 + (98596*y^2)/99225) - (314*x*y*cos((314*z)/315)*exp((4609615792312965*z)/18014398509481984))/(315*(((4609615792312965*y)/18014398509481984 - 1)^2 + (98596*y^2)/99225))== 0,- (314*x*y*sin((314*z)/315)*exp((4609615792312965*z)/18014398509481984))/(315*(((4609615792312965*y)/18014398509481984 - 1)^2 + (98596*y^2)/99225)) - (x*cos((314*z)/315)*exp((4609615792312965*z)/18014398509481984)*((4609615792312965*y)/18014398509481984 - 1))/(((4609615792312965*y)/18014398509481984 - 1)^2 + (98596*y^2)/99225)== -1];
xyz=vpasolve(eqns)
xyz.x
xyz.y
xyz.z
 

 

以上是全部的代码,显然这个方程组中的x很好求出来,但是跑完程序以后的x不对...这个方程组中的另外两个方程是根据虚部=0,实部=-1得到的

  • 写回答

0条回答 默认 最新

    报告相同问题?

    悬赏问题

    • ¥15 drone 推送镜像时候 purge: true 推送完毕后没有删除对应的镜像,手动拷贝到服务器执行结果正确在样才能让指令自动执行成功删除对应镜像,如何解决?
    • ¥15 求daily translation(DT)偏差订正方法的代码
    • ¥15 js调用html页面需要隐藏某个按钮
    • ¥15 ads仿真结果在圆图上是怎么读数的
    • ¥20 Cotex M3的调试和程序执行方式是什么样的?
    • ¥20 java项目连接sqlserver时报ssl相关错误
    • ¥15 一道python难题3
    • ¥15 牛顿斯科特系数表表示
    • ¥15 arduino 步进电机
    • ¥20 程序进入HardFault_Handler