Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- PIERWSZE:
- N=10 %wybieracie wymiar, praktycznie oZnacza to liczbe rezystorow bez R0
- 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
- R=randn(10,1,N+1);
- E(N+1)=0; %dodaje dodatkowe zero co pomaga w pozniejszych obliczeniach
- R(N+1)=0; %tutaj trza dodac nie 0 jeno R0
- for i=1:N
- MATR_E(i)=E(i)+E(i+1); %wypelniam wektor wg wyliczonej zaleznosci
- end
- MATR_E=MATR_E';
- MATR_R=zeros(N,N);
- for i=1:N
- MATR_R(i,i)=R(i)+R(i+1);
- if (i<N) %macierz rezystancji teZwynikajaca z Kirchoffa
- MATR_R(i+1,i)=-R(i+1);
- MATR_R(i,i+1)=-R(i+1);
- end
- end
- MATR_R(1:N,N+1)=MATR_E; %dodaje wektor E jako ostatnia kolumna macierzy R, aby korzystac z gaussa-jordana
- for i=1:N
- MATR_R(i,:)=MATR_R(i,:)/MATR_R(i,i);
- if (i<N)
- MATR_R(i+1,:)=MATR_R(i+1,:)-MATR_R(i,:)*MATR_R(i+1,i);
- end
- if (i>1)
- for k=1:N-1 %gauss jordan
- if (i-k>0)
- MATR_R(i-k,:)=MATR_R(i-k,:)-MATR_R(i,:)*MATR_R(i-k,i);
- end
- end
- end
- end
- MATR_R %wypisanie koncowe
- MATR_I(1:N,1)=MATR_R(1:N,N+1) %ostatnia kolumna pokazana jako wektor I (nasze rozwiazanie)
- --------------------------------------------------------------------------------------------------
- DRUGIE:
- tic
- N = 10;
- A = rand(10);
- b = rand(10);
- b = b(1:N,1);
- [L,U] = lu(A);
- y = inv(L)*b;
- x = inv(U)*y;
- toc
- DRUGIE2:
- A=[6 3 6;2 3 3;1 2 2];
- b=[30 17 11];
- % [N, M]= size(A);
- % L=eye(N);
- % U=zeros(N);
- % for i=1:N
- % for j=i:N
- % U(i,j)=A(i,j);
- % for k=1:(i-1)
- % U(i,j)=U(i,j)-L(i,k)*U(k,j);
- % end
- % end
- % for j=i+1:N
- % L(j,i)=A(j,i);
- % for k=1:(i-1);
- % L(j,i)=L(j,i)-L(j,k)*U(k,i);
- % end
- % L(j,i)=L(j,i)/U(i,i);
- % end
- % end
- --------------------------------------------------------------------------------
- CZWARTE:
- iks=[-1,-1,-1,-1,1 ,-1,1 ,-1,1, 1,1,-1,1, 1,-1, 1 ,-1 ,1 ,-1 ,-1 ,1]
- iks(22:42)=iks
- iks(43:63)=iks(1:21)
- L=32
- for i = 1:32
- for j = 1:32
- MATR(i,j)=iks(L+i-j)
- end
- end
- y=MATR*h'
- y0=awgn(y,0)
- y3=awgn(y,3)
- y5=awgn(y,5)
- y10=awgn(y,10)
- y100=awgn(y,100)
- -----------------------------------------------------------------------------------------
- PIĄTE:
- for i = 1:128
- for j = i+1: 128
- r(i)=data(j)*data(j-i);
- end
- end
- plot(r)
- autocorr(r)
- max=r(i)
- T=i
- for i=2:1577
- if(r(i)>max)
- max=r(i)
- T=i
- end
- end
- for i= 1: 10
- for j = 1 : 10
- R(i,j)=r(i+j-1)
- end
- end
- r10=r(1:10)
- r10=r10'
- a=-inv(R)*r10
- PIĄTE2:
- K=128;
- r=xcorr(data);
- r=r(1:128);
- plot(r);
- [m,T]=max(r); %m-warto¶ć T-miejsce
- P=10;
- R=zeros(10,10);
- for i=1:P
- R(:,i)=r(i:P+i-1);
- end
- r=r(1:10);
- a=inv(R)*r;
- M=24000;
- for i=1:24000
- if mod(i,T)==0
- e(i)=1;
- end
- end
- y(1:10)=0;
- for i=11:24000
- 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);
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement