Advertisement
szymcio93

townini2

Oct 21st, 2014
157
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.21 KB | None | 0 0
  1. PIERWSZE:
  2.  
  3. N=10 %wybieracie wymiar, praktycznie oZnacza to liczbe rezystorow bez R0
  4.  
  5. E=randn(10,1,N+1); %losuje wartosci R i I, normalnie powinny byc jakos ogarniete na tyle zeby nie byly np. rowne 0 lub ujemne
  6. R=randn(10,1,N+1);
  7.  
  8. E(N+1)=0; %dodaje dodatkowe zero co pomaga w pozniejszych obliczeniach
  9. R(N+1)=0; %tutaj trza dodac nie 0 jeno R0
  10.  
  11. for i=1:N
  12. MATR_E(i)=E(i)+E(i+1); %wypelniam wektor wg wyliczonej zaleznosci
  13. end
  14. MATR_E=MATR_E';
  15.  
  16.  
  17. MATR_R=zeros(N,N);
  18.  
  19. for i=1:N
  20. MATR_R(i,i)=R(i)+R(i+1);
  21. if (i<N) %macierz rezystancji teZwynikajaca z Kirchoffa
  22. MATR_R(i+1,i)=-R(i+1);
  23. MATR_R(i,i+1)=-R(i+1);
  24. end
  25. end
  26.  
  27.  
  28. MATR_R(1:N,N+1)=MATR_E; %dodaje wektor E jako ostatnia kolumna macierzy R, aby korzystac z gaussa-jordana
  29.  
  30.  
  31.  
  32. for i=1:N
  33. MATR_R(i,:)=MATR_R(i,:)/MATR_R(i,i);
  34. if (i<N)
  35. MATR_R(i+1,:)=MATR_R(i+1,:)-MATR_R(i,:)*MATR_R(i+1,i);
  36. end
  37. if (i>1)
  38. for k=1:N-1 %gauss jordan
  39. if (i-k>0)
  40. MATR_R(i-k,:)=MATR_R(i-k,:)-MATR_R(i,:)*MATR_R(i-k,i);
  41. end
  42. end
  43. end
  44. end
  45.  
  46. MATR_R %wypisanie koncowe
  47.  
  48. MATR_I(1:N,1)=MATR_R(1:N,N+1) %ostatnia kolumna pokazana jako wektor I (nasze rozwiazanie)
  49.  
  50. --------------------------------------------------------------------------------------------------
  51. DRUGIE:
  52.  
  53. tic
  54. N = 10;
  55. A = rand(10);
  56. b = rand(10);
  57. b = b(1:N,1);
  58. [L,U] = lu(A);
  59. y = inv(L)*b;
  60. x = inv(U)*y;
  61. toc
  62.  
  63. DRUGIE2:
  64.  
  65.  
  66. A=[6 3 6;2 3 3;1 2 2];
  67. b=[30 17 11];
  68.  
  69. % [N, M]= size(A);
  70. % L=eye(N);
  71. % U=zeros(N);
  72. % for i=1:N
  73. % for j=i:N
  74. % U(i,j)=A(i,j);
  75. % for k=1:(i-1)
  76. % U(i,j)=U(i,j)-L(i,k)*U(k,j);
  77. % end
  78. % end
  79. % for j=i+1:N
  80. % L(j,i)=A(j,i);
  81. % for k=1:(i-1);
  82. % L(j,i)=L(j,i)-L(j,k)*U(k,i);
  83. % end
  84. % L(j,i)=L(j,i)/U(i,i);
  85. % end
  86. % end
  87.  
  88. --------------------------------------------------------------------------------
  89. CZWARTE:
  90.  
  91. iks=[-1,-1,-1,-1,1 ,-1,1 ,-1,1, 1,1,-1,1, 1,-1, 1 ,-1 ,1 ,-1 ,-1 ,1]
  92. iks(22:42)=iks
  93. iks(43:63)=iks(1:21)
  94.  
  95. L=32
  96.  
  97. for i = 1:32
  98. for j = 1:32
  99. MATR(i,j)=iks(L+i-j)
  100. end
  101. end
  102.  
  103. y=MATR*h'
  104.  
  105. y0=awgn(y,0)
  106. y3=awgn(y,3)
  107. y5=awgn(y,5)
  108. y10=awgn(y,10)
  109. y100=awgn(y,100)
  110.  
  111. -----------------------------------------------------------------------------------------
  112.  
  113. PIĄTE:
  114.  
  115. for i = 1:128
  116. for j = i+1: 128
  117. r(i)=data(j)*data(j-i);
  118. end
  119. end
  120.  
  121. plot(r)
  122. autocorr(r)
  123.  
  124. max=r(i)
  125. T=i
  126.  
  127. for i=2:1577
  128. if(r(i)>max)
  129. max=r(i)
  130. T=i
  131. end
  132. end
  133.  
  134. for i= 1: 10
  135. for j = 1 : 10
  136. R(i,j)=r(i+j-1)
  137. end
  138. end
  139.  
  140. r10=r(1:10)
  141. r10=r10'
  142. a=-inv(R)*r10
  143.  
  144. PIĄTE2:
  145. K=128;
  146.  
  147. r=xcorr(data);
  148.  
  149. r=r(1:128);
  150.  
  151. plot(r);
  152.  
  153. [m,T]=max(r); %m-warto¶ć T-miejsce
  154.  
  155. P=10;
  156. R=zeros(10,10);
  157.  
  158. for i=1:P
  159.  
  160. R(:,i)=r(i:P+i-1);
  161.  
  162. end
  163. r=r(1:10);
  164.  
  165. a=inv(R)*r;
  166.  
  167. M=24000;
  168.  
  169. for i=1:24000
  170.  
  171. if mod(i,T)==0
  172.  
  173. e(i)=1;
  174.  
  175. end
  176. end
  177. y(1:10)=0;
  178.  
  179. for i=11:24000
  180.  
  181. y(i)=e(i)-a(1)*y(i-1)-a(2)*y(i-2)-a(3)*y(i-3)-a(4)*y(i-4)-a(5)*y(i-5)-a(6)*y(i-6)-a(7)*y(i-7)-a(8)*y(i-8)-a(9)*y(i-9)-a(10)*y(i-10);
  182.  
  183. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement