查看: 296|回复: 0|关注: 0

[未答复] slove 未定義x出現x

[复制链接]

新手

5 麦片

财富积分


050


1

主题

1

帖子

0

最佳答案
发表于 2019-2-11 14:45:32 | 显示全部楼层 |阅读模式
本帖最后由 PeterSu 于 2019-2-11 14:53 编辑

最近小弟使用 利用slove求解 矩陣方程 之解在未定一x的情況下 但結果出現了x
1為其中一個出現的結果
2為我的工作環境

我一開始以為是非常小的小數
所以我後來利用vpa去求解近似的小數
但出現的x卻沒有消失

請問有人遇過類似的問題嗎

代碼如下
c=1;
m=4;
n=3;
syms z u r R ee

z=1
u=2
r=2
%B
B=eye(m)
B=B*z

%A
A = cell(c,1);
E = cell(c,1);
Q=eye(m)
W=zeros(m)
W(1,1)=-1
W(1,m)=1

for i= 1:c
    E{i}=zeros(m);
  for j=2:m
     E{i}(j,j)=i;

    for k = 1:i
        if c>=i
        end
        E{i}(k+1,k+1)=k;  
    end   
  end  
end
eval(['A' num2str(0) ' = W*r +Q*(-z)']);
for i= 1:c
    A{i}=E{i}*(-u)+W*r+Q*(-z);

end

%C
C = cell(c,1);

for i= 1:c
    C{i}=zeros(m);

  for j=1:(m-1)
      C{i}((1+j),j)=i;
    for k = 1:i
        if c>=i
        end
        C{i}((1+k),k)=k;

    end   
  end  
end

for i = 1:c
   C{i}=C{i}*u ;
end



R= sym('r',[m m]);
eq=B+R*A{c}+(R^2)*C{c};
T=zeros(m);
sol = solve(eq==T,R,'Real',true);

矩陣

矩陣

工作環境

工作環境
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

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