Purchase Solution

Gauss-Seidel iteration in Matlab

Not what you're looking for?

Ask Custom Question

Hello

I am taking a numerical methods class and I have limited Matlab skills to be able to complete this problem by myself.

The problem is to linearize the Lorenz equations and then to write a Matlab program without using any built in Matlab solvers such as ODE45, etc. to create the Lorenz butterfly using linearized Lorenz equations with Gauss-Seidel iteration.

I have the linearized Lorenz equations, constants, stopping criteria and initial conditions here:

stopping criteria: 1.e-6
initial conditions: x = 2, y = 2, z = 4
sigma = 10, rho = 28, beta = 8/3
delta_t = 1.e-4

linearized lorenz equations:

x^n+1 - x^n = delta_t * sigma*(y^n+1 - x^n+1)
[please rewrite so equation begins x^n+1 = ...]

y^n+1 = (1 + delta_t*x^n+1*(rho - z^n+1)/(1 + delta_t)

z^n+1 = (1 + delta_t * x^n+1 * y^n+1)/(1 + beta * delta_t)

The code MUST only use Gauss-Seidel iteration to run the iterations (no RK2 or other methods). The code may not use any built in Matlab functions or solvers or matrices to store answers. It needs to tell how much CPU time was used (incorporate tic toc?) What was also suggested was to use a file to dump all results to while the iterations run, making the program run faster. With the results of the iterations a 3D plot needs to be plotted to show the Lorenz butterfly (attractor) using delta_t = 1.e-4. Then a second plot needs to be created that shows the x on the vertical axis and t on the horizontal axis plotting both delta_t = 1.e-4 and delta_t = 1.e-5 so you can see how a slight variation in delta_t causes the two lines to depart after a short while

Purchase this Solution
Solution Summary

This provides an example of creating a Matlab code for Gauss-Seidel iteration.

Solution Preview

The case with dt = 1e-4 runs 193 seconds
on my PC ( 6 years old, one 2GHz processor, 2GB memory),
so the case with dt= 1e-5 should probably run 10 times as much.
(!!! Do not forget to modify the values of dt1 and dt2 at the start of ...

Purchase this Solution
Free BrainMass Quizzes
Exponential Expressions

In this quiz, you will have a chance to practice basic terminology of exponential expressions and how to evaluate them.

Graphs and Functions

This quiz helps you easily identify a function and test your understanding of ranges, domains , function inverses and transformations.

Multiplying Complex Numbers

This is a short quiz to check your understanding of multiplication of complex numbers in rectangular form.

Geometry - Real Life Application Problems

Understanding of how geometry applies to in real-world contexts

Solving quadratic inequalities

This quiz test you on how well you are familiar with solving quadratic inequalities.

View More Free Quizzes
Related BrainMass Solutions

  • Jacobi and Gauss-Seidel iteration

    Start with P0 = 0 and use Jacobi iteration to find..... (Complete problem found in attachment) Attached This shows how to work with Jacobi and Gauss-Seidel iteration.

  • Linearizing Lorenz Equations using the Implicit Euler Method

    I need them linearized so that I can use Gauss-Seidel iteration in Matlab to create the butterfly effect (we are not to use ode45 or other solvers, hence the need for GS iteration and linearized equations).

  • Gaussian elimination and Gauss-Siedel

    (ii) Using the Gauss Seidel Method: The first step is to solve the first equation for W and the second for S when the system becomes: 0.25W + 0.21I = 1 0.015W + 0.04I = 1 If we assume the initial solution to be W = 0, I = 0, plugging in these

  • One Dimensional Dirichlet

    On each level of multigrid use Gauss-Seidel algorithm to solve the system of equations.

  • Gauss-Markov Simulated Annealing

    49263 Gauss-Markov Simulated Annealing Matlab code for Gauss-Markov Simulated Annealing This is a Matlab program for Gauss-Markov Simulated Annealing A Matlab code for Gauss-Markov simulated annealing.

  • Model Building and Parameter Estimation

    Write a MATLAB function that, given the time vector t and a value of a, returns a vector containing the gradient. Implement the Gauss-Newton algorithm as a MATLAB function that returns the estimated value of a and the resulting model y(t).

  • Solutions of Linear Equations

    10891 Solutions of Linear Equations I would like a short explanation of Gaussian Elimination with partial pivoting and Gauss-Seidel. Also, explain when each applies or when one is better than the other. Please include some examples.

  • Fixed Point Iteration Problem Using Matlab

    58834 Fixed Point Iteration Problem Using Matlab Given the function f(x) =1 + 0.5 sin x, show that it has a unique fixed point x* in the interval I=[0,2] and that the iterations xn+1=f(xn) converge to the fixed point for any x0 Є I.

  • Non-linear scalar equation(Matlab only)

    All the zero crossings should be in the plot. How many are there? 2) Write a program in Matlab that computes the roots using the fixed point iteration x_n+1=-sin(2x_n) +(5/4)x_n -(3/4).

View More Related Solutions