Recent content by henrybrent

Homework Statement A magnetic data set is believed to be dominated by a strong periodic tidal signal of known tidal period\Omega The field strength F(t) is assumed to follow the relation: F=a+b\cos\Omega t + c\sin\Omega t If the data were evenly spaced in time, then Fourier analysis would...

Ah, all it was was a simple case of adding contour(X(:,:,21),Y(:,:,21),wf(:,:,21));

I'm trying to generate contours that plot with the animation using contour(X,Y) but it's not quite working. I only need it to be in the x-y plane at 0km depth (see picture for greater clarification) % A script to plot two 3D waves. clear %Define the total time for simulation T_total=10; %Set...

Still haven't cracked it, haven't had time to work on it.

Still haven't managed to solve this. Will post what I've got soon when I'm home

Ok, thank you

u(x,t) = \sum\limits_{n=0}^{\infty} \sin \left( \frac{n\pi x_s}{L} \right) F(\omega_n)\sin(\frac{n\pi x}{L})\cos(\omega_nt) I think this is the equation we use. However, to plot the first few harmonics, we just use \sin \left( \frac{n\pi x_s}{L} \right) F(\omega_n)\sin(\frac{n\pi...

Almost, rather than think of it as two problems, I'm thinking of it as just 1. Write a script that sums the first 100? harmonics of the string when the source is applied, and plot the first 5 harmonics as well as the combination of 100. I should be able to use just the equation for the source...

The equation used to describe the standing wave is The code for the standing wave uses equation u(x,t) = A_nsin(\frac{n \pi x}{L})cos(\omega_n t) where \omega_n = \frac{n\pi v x}{L} the source is described as F( \omega_n) = e^\frac{ {(- \omega_n \tau)^2}}{4} When the source is...

I ran your code, it looks/runs great. My standing waves should be depending on at least the tuning parameter used in equation 1 on the sheet I attached. The standing waves result due to a source being applied (equation 1). The angular frequency is a variable that changes, yes, but I think this...

a_n = 1; L = 1; k(n) = @(n) n*pi/L; w_n = @(n,x,t) k(n).*x./t; x = 0:0.05:4; t = 1:0.01:20; for n = 1:4 for j=1:length(t) u = a_n*sin(k(n)* x).* cos(w_n(j,x,t(j)).*t(j)); %Formula for displacement pause (0.1) plot(u) axis([1,65,-1.5,1.5]) grid on end end Atm...

k = n*pi/L w_n = @(n,x,t) (k*x)*t; Getting a dimensional error

Ok, that makes sense. What is K in your formula? In mine, it was to represent n*pi/L, just a constant in my mind, but obviously that's wrong. In your code its still present as k=5?

Would Indexing instead of using nested for loops be a better/more efficient way of doing this?

So this is where I apply the loop over the initial loop for n = 1:4, to separate the harmonics. The code looks like a load of harmonics looping over and over again.