[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [obm-l] Multiplicação de matrizes no determinante



pelo que entendi, essa resposta tá invertendo matriz não quadrada (B) e 
fazendo uma divisão de matrizes... é isso mesmo? 


Em (18:04:45), obm-l@mat.puc-rio.br escreveu: 


>desse jeito nao ta certo nao 
> det(1+AB)=det(1+AB)B^-1/B^-1=det(B^-1+A)/B^-1 
> agora multiplica no lado esquerdo por B 
> det(BB^-1+BA)/BB^-1 
> BB^-1=I 
> det(I+BA)=det(I+AB) 
> On 6/4/07, edneiramaral  wrote: 
> Há um tempo atrás tinha mandado esse problema na lista... Encontrei a 
>resposta e queria compartilhar com vcs. 
> 
>Basicamente, o problema pode ser reduzido a mostrar que: 
>det (I + AB) = det (I + BA) 
>qnd A e B não são quadradas. Digamos: 
>dim(A) = M x N 
>dim(b) = N x M 
> 
>Usei essa dica: 
> 
>http://www.mathematik.uni-bielefeld.de/~sillke/PUZZLES/det-i-ab 
> 
>e a propriedade de determinantes de matrizes particionadas (matrizes 
>definidas por partes): 
>det( [A B] ) = det( [A 0] ) = det(A).det(C) 
> [0 C] [B C] 
> 
>Valeu! 
> 
>"Marcelo Salhab Brogliato" wrote: 
>Opa, 
>é verdade! vou pensar melhor aqui.. 
>qualquer ideia eu mando amanha!! 
>abracos, 
>Salhab 
> 
>On 4/30/07, edneiramaral edneiramaral@ig.com.br > wrote: 
>> Acrescentando mais um dado, que existe no problema que estou trabalhando: 
>> R é tal que 
>> Rij = conj(Rji) 
>> 
>> Resposta ao Salhab: 
>> 
>> Se bem entendi sua idéia, eu cheguei a pensar algo parecido, mas parei pq 
> 
>as 
>> matrizes F e H não são quadradas e nesse caso o determinante não está 
>> definido, correto? 
>> 
>> Consigo por exemplo passar o R pro final e ficar com: 
>> 
>> det( I + R.H.F.F*.HH.F.F*.H*.R) 
>> 
>> (usando ainda que, no caso a que estou aplicando, R tem a propriedade 
>acima) 
>> 
>> Mas só faço isso pq tanto R (MxM) quanto as matrizes H.F.F*.H* e 
>H*.F*.F.H 
>> são quadradas. Mas passar o R pro meio da multiplicação eu não consigo 
> 
>> porque H.F ou F*.H* não são quadradas. 
>> 
>> Obrigado, 
>> Ednei Amaral 
>> 
>> 
>> Em (14:42:47), obm-l@mat.puc-rio.br escreveu: 
>> 
>> 
> 
>> >Olá, 
>> > 
>> >queremos mostrar que: 
>> >det( I + R.H.F.F*.H*) = det( I + F*.H*.R.H.F) 
>> > 
>> >sabemos que det(F) = conj[det(F*)] ... onde conj é conjugado do numero 
> 
>> >complexo 
>> > 
>> >assim: 
>> >det(H*)/conj[det(H)] . det(F*)/conj[det(F)] . det( I + R.H.F.F*.H*) = 
>> >det( F*H* + F*H*RHFF*H*) . 1/conj[det(H)] . 1/conj[det(F)] = det( I + 
>> >F*H*RHF) . det(H*)/conj[det(H)] . det(F*)/conj[det(F)] = det(I + 
> 
>> >F*H*RHF) 
>> > 
>> >o q eu fiz foi o seguinte.. multipliquei por det(H*)/conj[det(H)], que 
>> >é igual a 1, joguei o det(H*) pra dentro... coloquei dps o H* em 
>> >evidencia pela direita... tirei o det(H*) e simplifiquei... fiz isso 
> 
>> >com F e H.. 
>> > 
>> >espero que tenha dado pra entender 
>> > 
>> >abracos, 
>> >Salhab 
>> > 
>> >On 4/30/07, edneiramaral wrote: 
>> >> Olá, 
>> >> 
> 
>> >> estava verificando um resultado apresentado numa demonstaração e 
>cheguei 
>> a 
>> >> um resultado semelhante ao que queria provar. Na prática, o resultado 
é 
>> >> igual ao que cheguei (conforme verifiquei com alguns testes 
numéricos), 
> 
>> >> porém a forma apresentada está diferente. 
>> >> 
>> >> Gostaria de saber como faço para mostrar a seguinte igualdade: 
>> >> 
>> >> det( I + R.H.F.F*.H*) = det( I + F*.H*.R.H.F) 
> 
>> >> 
>> >> onde 
>> >> . significa multiplicação 
>> >> * significa conjungado transposto da matriz (hermitiano) 
>> >> H é matriz M x N 
>> >> R é matriz M x M 
> 
>> >> F é matriz N X P 
>> >> I é matriz identidade de tamanho compatível com a outra parcela da 
soma 
>> >> 
>> >> Obrigado, 
>> >> Ednei Amaral 
>> >> 
>> >> 
> 
>> >> 
>> > 
>> 
>========================================================================= 
>> >Instruções para entrar na lista, sair da lista e usar a lista em 
>> > 
>http://www.mat.puc-rio.br/~nicolau/olimp/obm-l.html 
>> 
>========================================================================= 
>> > 
>> >---------- 
>> 
>> 
>> 
> 
>----------