| ²é¿´: 494 | »Ø¸´: 1 | ||
Ìïɽ¶«¾èÖú¹ó±ö (ÖøÃûдÊÖ)
|
[ÇóÖú]
H*pinv(H)ΪʲôÔڷǶԽÇÔª³öÏÖ·ÇÁãÖµ£¿
|
|
ÈçÌ⣬´úÂëÈçÏ¡£¶àл¸÷λͬѧ°¡£¡£¡ clc; clear; N=8; M=5; A1=zeros(N,N); A2=zeros(N,N); B1=zeros(N,N); B2=zeros(N,N); C1=zeros(N,N); C2=zeros(N,N); H=zeros(M*N,M*N); Gr=zeros(M*N,M*N); t=2.8; %É趨zigzag±ß½çϵͳ¹þÃܶپØÕó for j=1:N-1 if(mod(j,2)==0) A1(j,j+1)=-t; A1(j+1,j)=-t; end if(mod(j,2)==1) A2(j,j+1)=-t; A2(j+1,j)=-t; end end for j=1:N B1(j,j)=-t; B2(j,j)=-t; end C2=B1'; C1=B2'; %ϵͳÕûÌå¹þÃܶپØÕó for j=1:M if(mod(j,2)==1) H((j-1)*N+1:j*N, (j-1)*N+1:j*N)=A1; else H((j-1)*N+1:j*N, (j-1)*N+1:j*N)=A2; end end for j=1:M if(mod(j,2)==1) if j==1 H((j-1)*N+1:j*N,j*N+1: (j+1)*N)=B1; else if j<M H((j-1)*N+1:j*N,j*N+1: (j+1)*N)=B1; H((j-1)*N+1:j*N, (j-2)*N+1: (j-1)*N)=C1; else H((j-1)*N+1:j*N, (j-2)*N+1: (j-1)*N)=C1; end end else H((j-1)*N+1:j*N,j*N+1: (j+1)*N)=B2; H((j-1)*N+1:j*N, (j-2)*N+1: (j-1)*N)=C2; end end Gr=pinv(H); HH=Gr*H; ÕâÀïÓ¦¸ÃµÃµ½µ¥Î»¾ØÕ󣬵«ÊÇÔÚHH·Ç¶Ô½ÇÔªµÃµ½·ÇÁãÖµ¡£ºÜÆæ¹Ö °¡£¡£¡Çó¸ßÈ˽â´ð£¡£¡ [ Last edited by Ìïɽ¶« on 2013-8-22 at 16:05 ] |
» ²ÂÄãϲ»¶
268Çóµ÷¼Á
ÒѾÓÐ9È˻ظ´
301Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
0703»¯Ñ§336·ÖÇóµ÷¼Á
ÒѾÓÐ5È˻ظ´
301Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
²ÄÁÏר˶306Ó¢Ò»Êý¶þ
ÒѾÓÐ9È˻ظ´
304Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
²ÄÁϹ¤³Ìר˶µ÷¼Á
ÒѾÓÐ5È˻ظ´
344Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
326Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
296Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´

alorsmail
Ìú¸Ëľ³æ (ÖøÃûдÊÖ)
- Ó¦Öú: 52 (³õÖÐÉú)
- ½ð±Ò: 7398.3
- ºì»¨: 41
- Ìû×Ó: 2104
- ÔÚÏß: 237Сʱ
- ³æºÅ: 2537936
- ×¢²á: 2013-07-09
- רҵ: ϵͳÉúÎïѧ
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xzhdty: ½ð±Ò+1, ¸Ðл²ÎÓë 2013-08-22 20:02:01
Ìïɽ¶«: ½ð±Ò+5, ¡ï¡ï¡ïºÜÓаïÖú, ¶àл»Ø¸´£¬Ò»»áÃþË÷Ï 2013-08-22 20:10:03
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
xzhdty: ½ð±Ò+1, ¸Ðл²ÎÓë 2013-08-22 20:02:01
Ìïɽ¶«: ½ð±Ò+5, ¡ï¡ï¡ïºÜÓаïÖú, ¶àл»Ø¸´£¬Ò»»áÃþË÷Ï 2013-08-22 20:10:03
|
ÕâÀïÓ¦¸ÃÊÇÎó²î£¬ Èç¹ûÄã¶ÔÔ HH ×öÒÔϱ任µÄ»°£¬¾Í»áµÃµ½µ¥Î»¾ØÕó threshold = 0.5; HH1 = double(HH>=threshold); µ±È»£¬Èç¹ûÄãÈ¡ threshold =0.1; µÄ»°£¬Äã»áµÃµ½ÁíÍâÁ½¸ö·ÇÁã ¶Ô½Ç Ñо¿Ò»ÏÂΪʲô threshold »áÓ°ÏìÄãµÄ½á¹û£¬ |
2Â¥2013-08-22 17:04:07













»Ø¸´´ËÂ¥