Fdtxjr

Why is delta-v is the most useful quantity for planning space travel?

Should my PhD thesis be submitted under my legal name?

Visiting the UK as unmarried couple

Meta programming: Declare a new struct on the fly

How can I successfully establish a nationwide combat training program for a large country?

Pronouncing Homer as in modern Greek

Superhero words!

Lightning Web Component - do I need to track changes for every single input field in a form

Should a half Jewish man be discouraged from marrying a Jewess?

Indicating multiple different modes of speech (fantasy language or telepathy)

Can the harmonic series explain the origin of the major scale?

Can I Retrieve Email Addresses from BCC?

How will losing mobility of one hand affect my career as a programmer?

Organic chemistry Iodoform Reaction

Female=gender counterpart?

Giant Toughroad SLR 2 for 200 miles in two days, will it make it?

Can a Gentile theist be saved?

Why isn't KTEX's runway designation 10/28 instead of 9/27?

Resetting two CD4017 counters simultaneously, only one resets

How to color a zone in Tikz

Is there an Impartial Brexit Deal comparison site?

Was the picture area of a CRT a parallelogram (instead of a true rectangle)?

"lassen" in meaning "sich fassen"

Calculating the number of days between 2 dates in Excel

How to solve a matrix using Cholesky Decompositon on Matlab

Implement a program in Matlab for LU decomposition with pivotingSignificant improvement when I use lsqnonlin function with wrong sized XHow to solve linear equation using Matlab?How to solve over-determined linear system of equations?LU Decomposition vs. QR Decomposition for similar problemsusing MATLAB to solve a linear programHow to solve Newton's method in matlab?Solve intregal equation using MATLABIdentify a state space model from measured inputs and outputsUsing MATLAB to solve Poisson matrix equation

0

$begingroup$

I'm trying to solve a system of linear equations on MATLAB, I have written a code that solves the problem using Gaussian Elimination. But I was wondering how I could modify this to use other methods of matrix decomposition, such as Cholesky Decomposition?

%set up matrices
A=[-900 400 0 0 0; 500 -900 400 0 0; 0 500 -900 400 0; 0 0 500 -900 400; 0 0 0 400 -900];
b=[0;0;0;0;0];

%decompose A matrix
[L, U] = lu(A);

%set up series of yin's for plot
yin=0:0.01:0.25;
yout=zeros(length(yin),1);
xout=zeros(length(yin),1);

%create b based on yin, and solve
for i=1:1:length(yin)
 b(1)=-500*yin(i);
 d = Lb;
 y = Ud;
 %save yout and xout values before next iteration
 yout(i)=y(5);
 xout(i)=y(1)*4;
end

matlab

share|cite|improve this question

asked Mar 16 at 22:48

eenzeenz

11

$endgroup$

add a comment | 

0

$begingroup$

I'm trying to solve a system of linear equations on MATLAB, I have written a code that solves the problem using Gaussian Elimination. But I was wondering how I could modify this to use other methods of matrix decomposition, such as Cholesky Decomposition?

%set up matrices
A=[-900 400 0 0 0; 500 -900 400 0 0; 0 500 -900 400 0; 0 0 500 -900 400; 0 0 0 400 -900];
b=[0;0;0;0;0];

%decompose A matrix
[L, U] = lu(A);

%set up series of yin's for plot
yin=0:0.01:0.25;
yout=zeros(length(yin),1);
xout=zeros(length(yin),1);

%create b based on yin, and solve
for i=1:1:length(yin)
 b(1)=-500*yin(i);
 d = Lb;
 y = Ud;
 %save yout and xout values before next iteration
 yout(i)=y(5);
 xout(i)=y(1)*4;
end

matlab

share|cite|improve this question

asked Mar 16 at 22:48

eenzeenz

11

$endgroup$

add a comment | 

$begingroup$

I'm trying to solve a system of linear equations on MATLAB, I have written a code that solves the problem using Gaussian Elimination. But I was wondering how I could modify this to use other methods of matrix decomposition, such as Cholesky Decomposition?

%set up matrices
A=[-900 400 0 0 0; 500 -900 400 0 0; 0 500 -900 400 0; 0 0 500 -900 400; 0 0 0 400 -900];
b=[0;0;0;0;0];

%decompose A matrix
[L, U] = lu(A);

%set up series of yin's for plot
yin=0:0.01:0.25;
yout=zeros(length(yin),1);
xout=zeros(length(yin),1);

%create b based on yin, and solve
for i=1:1:length(yin)
 b(1)=-500*yin(i);
 d = Lb;
 y = Ud;
 %save yout and xout values before next iteration
 yout(i)=y(5);
 xout(i)=y(1)*4;
end

matlab

share|cite|improve this question

asked Mar 16 at 22:48

eenzeenz

11

$endgroup$

%set up matrices
A=[-900 400 0 0 0; 500 -900 400 0 0; 0 500 -900 400 0; 0 0 500 -900 400; 0 0 0 400 -900];
b=[0;0;0;0;0];

%decompose A matrix
[L, U] = lu(A);

%set up series of yin's for plot
yin=0:0.01:0.25;
yout=zeros(length(yin),1);
xout=zeros(length(yin),1);

%create b based on yin, and solve
for i=1:1:length(yin)
 b(1)=-500*yin(i);
 d = Lb;
 y = Ud;
 %save yout and xout values before next iteration
 yout(i)=y(5);
 xout(i)=y(1)*4;
end

share|cite|improve this question

asked Mar 16 at 22:48

eenzeenz

11

share|cite|improve this question

asked Mar 16 at 22:48

eenzeenz

11

asked Mar 16 at 22:48

eenzeenz

11

asked Mar 16 at 22:48

eenzeenz

11

1 Answer
1

active

oldest

votes

0

$begingroup$

Given a self-adjoint positive-definite square matrix $A$, $operatornamechol(A)$ returns an upper triangular matrix $U$ such that $A = U^* U$

To turn this into code, replace

[L, U] = lu(A);

by

U = chol(A);
L = U';

share|cite|improve this answer

answered yesterday

Alex VongAlex Vong

1,309819

$endgroup$

add a comment | 

Your Answer

StackExchange.ifUsing("editor", function ()
return StackExchange.using("mathjaxEditing", function ()
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
);
);
, "mathjax-editing");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "69"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3150943%2fhow-to-solve-a-matrix-using-cholesky-decompositon-on-matlab%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

0

$begingroup$

Given a self-adjoint positive-definite square matrix $A$, $operatornamechol(A)$ returns an upper triangular matrix $U$ such that $A = U^* U$

To turn this into code, replace

[L, U] = lu(A);

by

U = chol(A);
L = U';

share|cite|improve this answer

answered yesterday

Alex VongAlex Vong

1,309819

$endgroup$

add a comment | 

0

$begingroup$

Given a self-adjoint positive-definite square matrix $A$, $operatornamechol(A)$ returns an upper triangular matrix $U$ such that $A = U^* U$

To turn this into code, replace

[L, U] = lu(A);

by

U = chol(A);
L = U';

share|cite|improve this answer

answered yesterday

Alex VongAlex Vong

1,309819

$endgroup$

add a comment | 

$begingroup$

Given a self-adjoint positive-definite square matrix $A$, $operatornamechol(A)$ returns an upper triangular matrix $U$ such that $A = U^* U$

To turn this into code, replace

[L, U] = lu(A);

by

U = chol(A);
L = U';

share|cite|improve this answer

answered yesterday

Alex VongAlex Vong

1,309819

$endgroup$

[L, U] = lu(A);

U = chol(A);
L = U';

share|cite|improve this answer

answered yesterday

Alex VongAlex Vong

1,309819

answered yesterday

Alex VongAlex Vong

1,309819

answered yesterday

Alex VongAlex Vong

1,309819