i need change this to MATLAB CODE , its Necessary

x0 = [1 1;1 1];
A = [32 55; 1.2 1.5; 1.9 3.2; 5.2 4.3; 2.7 3.3;3.4 2.1];
Q = [5 3.1;2.2 4];
C = rand(5,5);
disp(' ')
disp('do')
load handel
disp('yourow')
sound(y, Fs)
disp('ndamnw')
disp('o')
try
    A = A+Q
    kron(13)
    solver(X_k_+1)
    type('Kronecker')
end
disp('rk')

Anyone on here want to run this to validate that I did it right?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Martin 2014-7-31

i write it myself :

clc close all clear all

x0=[1 1 ;1 1 ]; % A=[3.2 5.5; 1.2 1.5; 1.9 3.2 ; 5.2 4.3; 2.7 3.3 ; 3.4 2.1]; A=[3.2 5.5; 1.2 1.5; 1.9 3.2 ; 5.2 4.3];

Q=[5 3.1 ; 2.2 4]; % C=[2.2 3.2 1.3 5.1 4.8 5.9 % 4 2 3 8 5 1 % 1.1 1.5 2.1 2.2 3.7 3.8 % 2.4 3.4 4.4 5.5 3.5 3 % 1 2 4 5 3 2 % 1.2 2.3 3.5 5.3 4 1.8]; C=rand(4,4); H=1;X=[];x=x0; for k=1:100 I=eye(length(x)); x=Q+A'*inv((kron(I,x)-C) )*A; X=[X x]; end

its for c => 4*4 AND A=> 4*2 but i need c => 6*4 AND A=> 6*2

Answer 2

Joseph Cheng 2014-7-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/143629-i-need-change-this-to-matlab-code-its-necessary#answer_146678

the functions you'll need to use are kron() for the kronecker and ctranspose() or the symbol ' and eye() for the identity matrix. a for loop of the equation will work through your iterations.