کد های رایگان تجزیه چولسکی + پروژه کامل و رایگان + نرم افزار متلب
توضیحات:
در این پروژه به بررسی کد برنامه تجزیه چولسکی در متلب به رایگان پرداخته شده است. کد برنامه بطور کامل توصیف گردیده است.
کد های زیر را در ام فایل متلب کپی کرده و فایل اجرایی دوم را ران کنید.
function [F]=cholesky(A,option)
%Function to find the Cholesky factor of a Positive Definite matrix A
Licensed under Creative Commons: CC-NC-BY-SA 3.0
%A = positive definite matrix
%Option can be one of the following 'Lower','Upper'
%L = Cholesky factorizaton of A such that A=LL^T
%If option ='Lower', then it returns the Cholesky factored matrix L in
%lower triangular form
%If option ='Upper', then it returns the Cholesky factored matrix L in
%upper triangular form
%Test for positive definiteness (symmetricity need to satisfy)
%Check if the matrix is symmetric
if ~isequal(A,A'),
error('Input Matrix is not Symmetric');
end
if isPositiveDefinite(A),
[m,n]=size(A);
L=zeros(m,m);%Initialize to all zeros
row=1;col=1;
j=1;
for i=1:m,
a11=sqrt(A(1,1));
L(row,col)=a11;
if(m~=1), %Reached the last partition
L21=A(j+1:m,1)/a11;
L(row+1:end,col)=L21;
A=(A(j+1:m,j+1:m)-L21*L21');
[m,n]=size(A);
row=row+1;
col=col+1;
end
end
switch nargin
case 2
if strcmpi(option,'upper'),F=L';
else
if strcmpi(option,'lower'),F=L;
else error('Invalid option');
end
end
case 1
F=L;
otherwise
error('Not enough input arguments')
end
else
error('Given Matrix A is NOT Positive definite');
end
end
برنامه دوم:
function x=isPositiveDefinite(A)
%Function to check whether a given matrix A is positive definite
%Licensed under Creative Commons: CC-NC-BY-SA 3.0
%Returns x=1, if the input matrix is positive definite
%Returns x=0, if the input matrix is not positive definite
[m,~]=size(A);
%Test for positive definiteness
x=1; %Flag to check for positiveness
for i=1:m
subA=A(1:i,1:i); %Extract upper left kxk submatrix
if(det(subA)<=0); %Check if the determinent of the kxk submatrix is +ve
x=0;
break;
end
end
%For debug
%if x, display('Given Matrix is Positive definite');
%else display('Given Matrix is NOT positive definite'); end
end
دیگر پروژه های مرتبط:
شبیه سازی مقاله مرتبط با تجزیه چولسکی برای ماتریس های kernel
*************************
کانال قطعات موبایل اسنپ زیر نظر سایت متلب تولز افتتاح شد. عضو شوید تا از موارد زیر با قیمت ارزان مطلع شوید:
رفتن به لینک سایت قطعات موبایل
همچنین کانال تلگرامی ما را نیز دنبال نمائید:
یا snap_tools@ را در تلگرام جستجو نمایید.
دارای نماد اعتماد الکترونیکی
مجموعه: آموزش دستورات متلب, ریاضیات و محاسبات عددی