查看: 144|回复: 3|关注: 0

[已答复] 多元方程求解,传递函数

[复制链接]

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
发表于 2020-6-28 13:29:17 | 显示全部楼层 |阅读模式
       由已知传递函数G,G1和已知的公式g,利用待定系数法求未知传递函数G2,对于一般简单的系数为整数的传递函数利用下面的代码可以求出,但是对于像下文中这种系数都不是整数的,点击运行后一直处于求解状态,这个是怎么回事?
对于待定系数解方程组察看链接
多元方程求解
https://www.ilovematlab.cn/thread-600067-1-1.html
(出处: MATLAB中文论坛)






clc;
clear;

syms s x1 x2 x3 x4 x5 x6 x7 x8 x9 x10


% G1
numb1=[0.995674859669345 667.446369410333 236687.454351649];
denb1=[1 669.431861889998 242335.214127284];
g1=tf(numb1,denb1);
a=poly2sym(numb1,s)/poly2sym(denb1,s);
G1=vpa(a,15);


% %12.12  99.95%总G
numb=[0.0589876439904560,-294.729931998210,2079737.14084625,-8467097458.52300,32043389853130.9,-1.10500292141589e+17,1.67899243282274e+20,-2.72761702273197e+23,1.03729675923455e+27,-3.32746023363210e+30,-3.63887254133730e+33,1.86700553885761e+37,1.28873501289187e+40];
denb=[1,758.759924546565,63867789.8129226,45789820518.3603,1.38028265979850e+15,9.38739501214050e+17,1.19963549759688e+22,7.78528800328017e+24,3.85331119386103e+28,2.40413552592136e+31,3.08685181495285e+34,1.86193854863994e+37,1.53317071689829e+38];
% G=tf(numb,denb);
b=poly2sym(numb,s)/poly2sym(denb,s);
G=vpa(b,15);


G2=(x5*s^4+x4*s^3+x3*s^2+x2*s+x1)/(x10*s^4+x9*s^3+x8*s^2+x7*s+x6);


g=(G2+2*G1*G2+G1^2*G2)/(1-G1^2*G2^2);
[n,d]=numden(g);


n2=collect(n);d2=collect(d);
n3=coeffs(n2,s)/1e70;d3=coeffs(d2,s)/1e70;
n4=simplify(n3);d4=simplify(d3);
k=n4./d4-numb./denb;

for i=1:13
y(i)=k(i);
end
y=[y(1),y(2),y(3),y(4),y(5),y(6),y(7),y(8),y(9),y(10),y(11),y(12),y(13)]';
s=solve(y==0,[x1,x2,x3,x4,x5,x6,x7,x8,x9,x10])

回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

6

主题

1440

帖子

307

最佳答案
  • 关注者: 70
发表于 2020-6-28 15:07:30 | 显示全部楼层
我这边求出无解,试试其他函数,或者1STOPT软件
1.PNG
回复此楼 已获打赏: 0 积分

举报

新手

16 麦片

财富积分


050


9

主题

20

帖子

0

最佳答案
 楼主| 发表于 2020-6-28 15:14:10 | 显示全部楼层
20141303 发表于 2020-6-28 15:07
我这边求出无解,试试其他函数,或者1STOPT软件

s..x1 返回空解吗?
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

6

主题

1440

帖子

307

最佳答案
  • 关注者: 70
发表于 2020-6-28 16:19:50 | 显示全部楼层
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表