Теория информации. Лабораторный практикум в MATLAB. Учебное пособие

disp(P) %Ɂɚɞɚɧɢɟ ɜɟɤɬɨɪɚ ɧɚɱɚɥɶɧɵɯ ɜɟɪɨɹɬɧɨɫɬɟɣ ɫɨɫɬɨɹɧɢɣ p0 p=[1 0 0]; disp('ȼɟɤɬɨɪ ɧɚɱɚɥɶɧɵɯ ɜɟɪɨɹɬɧɨɫɬɟɣ ɫɨɫɬɨɹɧɢɣ p0:') disp(p) %Ɂɚɞɚɧɢɟ ɱɢɫɥɚ ɲɚɝɨɜ Ʉ K=15; disp('ɑɢɫɥɨ ɲɚɝɨɜ K:') disp(K) %ȼɵɱɢɫɥɟɧɢɹ ɩɨ ɲɚɝɚɦ for k=1:K formatSpec = '%i-ɣ ɲɚɝ:\n'; fprintf(formatSpec,k); disp('---------') pk=p*P^k; disp('ȼɟɤɬɨɪ ɜɟɪɨɹɬɧɨɫɬɟɣ ɫɨɫɬɨɹɧɢɣ:') disp(pk) disp('Ɇɧɨɝɨɲɚɝɨɜɚɹ ɦɚɬɪɢɰɚ ɩɟɪɟɯɨɞɧɵɯ ɜɟɪɨɹɬɧɨɫɬɟɣ:') disp(P^k) end %ȼɵɱɢɫɥɟɧɢɟ ɮɢɧɚɥɶɧɵɯ ɜɟɪɨɹɬɧɨɫɬɟɣ, %ɷɧɬɪɨɩɢɢ ɢ ɢɡɛɵɬɨɱɧɨɫɬɢ ɢɫɬɨɱɧɢɤɚ disp('ȼɵɱɢɫɥɟɧɢɟ ɮɢɧɚɥɶɧɵɯ ɜɟɪɨɹɬɧɨɫɬɟɣ:') disp('---------------------------------') %Ɉɩɪɟɞɟɥɟɧɢɹ ɩɨɪɹɞɤɚ ɦɚɬɪɢɰɵ Ɋ [~,m]=size(P); %Ɂɚɞɚɧɢɟ ɦɚɬɪɢɰɵ ɥɢɧɟɣɧɨɣ ɫɢɫɬɟɦɵ A=P'-eye(m); A(m,:)=ones(1,m); disp('Ɇɚɬɪɢɰɚ ɥɢɧɟɣɧɨɣ ɫɢɫɬɟɦɵ:') disp(A) %Ɂɚɞɚɧɢɟ ɜɟɤɬɨɪɚ ɫɜɨɛɨɞɧɵɯ ɱɥɟɧɨɜ b=[zeros(1,m-1) 1]'; disp('ȼɟɤɬɨɪ ɫɜɨɛɨɞɧɵɯ ɱɥɟɧɨɜ:') disp(b) %Ɋɟɲɟɧɢɟ ɫɢɫɬɟɦɵ ɥɢɧɟɣɧɵɯ ɭɪɚɜɧɟɧɢɣ pf=linsolve(A,b); disp('ȼɟɤɬɨɪ ɮɢɧɚɥɶɧɵɯ ɜɟɪɨɹɬɧɨɫɬɟɣ:') disp(pf) disp('ȼɵɱɢɫɥɟɧɢɟ ɷɧɬɪɨɩɢɢ ɢɫɬɨɱɧɢɤɚ:') disp('-------------------------------') disp('ɑɚɫɬɧɵɟ ɭɫɥɨɜɧɵɟ ɷɧɬɪɨɩɢɢ, ɛɢɬ/ɫɢɦɜɨɥ::')