This document contains solutions to four MATLAB exercises:
1) Generate the first 100 terms of a recursively defined sequence.
2) Write an M-file function to calculate values of a damped harmonic oscillator equation based on user inputs.
3) The function is tested for two cases and the output is plotted from 0 to 10 seconds.
4) A function is defined and values from 0 to 2π are calculated and plotted to graph the function Y(x)=exp(-x)sin(8x).
This document contains solutions to four MATLAB exercises:
1) Generate the first 100 terms of a recursively defined sequence.
2) Write an M-file function to calculate values of a damped harmonic oscillator equation based on user inputs.
3) The function is tested for two cases and the output is plotted from 0 to 10 seconds.
4) A function is defined and values from 0 to 2π are calculated and plotted to graph the function Y(x)=exp(-x)sin(8x).
This document contains solutions to four MATLAB exercises:
1) Generate the first 100 terms of a recursively defined sequence.
2) Write an M-file function to calculate values of a damped harmonic oscillator equation based on user inputs.
3) The function is tested for two cases and the output is plotted from 0 to 10 seconds.
4) A function is defined and values from 0 to 2π are calculated and plotted to graph the function Y(x)=exp(-x)sin(8x).
This document contains solutions to four MATLAB exercises:
1) Generate the first 100 terms of a recursively defined sequence.
2) Write an M-file function to calculate values of a damped harmonic oscillator equation based on user inputs.
3) The function is tested for two cases and the output is plotted from 0 to 10 seconds.
4) A function is defined and values from 0 to 2π are calculated and plotted to graph the function Y(x)=exp(-x)sin(8x).
Download as DOCX, PDF, TXT or read online from Scribd
Download as docx, pdf, or txt
You are on page 1/ 7
Exercise 1: Use MATLAB to generate the first 100 terms in the sequence a (n) define recursively by
(𝑛 + 1) = 𝑝 ∗ (𝑛) ∗ (1 − (𝑛)) with p=2.9 and a(1) = 0.5
Solution: Code: p=2.9; a(1)=0.5; for n=1:100 a(n+1)=p*a(n)*(1-a(n)) end
Result: Solution: Code: for Tc=0:10:200 Tf=9/5*Tc+32 end
Result: Exercise 2: MATLAB M-file Function Consider the following equation (𝑡)=𝑦(0)√1−𝜁𝑒−𝜁𝜔𝑛𝑡sin(𝜔𝑛√1−𝜁2∗𝑡+𝜃) a) Write a MATLAB M-file function to obtain numerical values of y(t). Your function must take (0), 𝜁, 𝜔𝑛, t and 𝜃 as function inputs and (𝑡) as output argument. b) Obtain the plot for y(t) for 0<t<10 with an increment of 0.1, by considering the following two cases
Case 1: (0)=0.15 m, 𝜔𝑛=√2 rad/sec, 𝜁=3/(2√2) and 𝜃=0;
Case 2: (0)=0.15 m, 𝜔𝑛=√2 rad/sec, 𝜁=1/(2√2) and 𝜃=0;
Code:
%Lab #1
%MatLab M-file Function Exersice
%{
This program will ask users to input values for the following equation
The above output demonstrate that Mat Lab can be use to obtain numerical values of an equation such as (𝑡)=𝑦(0)√1−𝜁𝑒−𝜁𝜔𝑛𝑡sin(𝜔𝑛√1−𝜁2∗𝑡+𝜃) Evaluated at the users specified values.
Graph plot part b:
The above graph represents the output plot of the equation
y(t) = 𝑦(0)√1−𝜁𝑒−𝜁𝜔𝑛𝑡sin(𝜔𝑛√1−𝜁2∗𝑡+𝜃)
With the conditions of y0=0.15 m, ωn = √ rad/sec, ζ = 3/(2√ ) and θ = 0
Exercise 4: Plot the function Y(x)=exp(-x)sin(8x) for 0<=x<=2π Solution: Code: Function y = func(x) % function decleration
