Can't sort this code

6 views (last 30 days)
Tom
Tom on 25 Jan 2012
The original question has been removed by the user.
  4 Comments
Jan
Jan on 26 Jan 2012
Please take the time to formulate a complete question. Currently it is not clear, what the program does, if it fails or not, if it creates the wanted results, or if not, what the dfferences between the results and your expectations are.
Tom
Tom on 26 Jan 2012
Okay yeah sorry, I'll try to bear that in mind.

Sign in to comment.

Accepted Answer

Dr. Seis
Dr. Seis on 26 Jan 2012
You define fnj as MxN, but in your for loop you are assigning values to fnj(j,n) where:
j = 1 : N
and
n = 1 : M
I think you meant to do "j = 1 : M" and "n = 1 : N"
Also make the change suggested by proecsm, namely
Fn = sum(fnj,2)
  2 Comments
Tom
Tom on 26 Jan 2012
Yeah I think that was like that because I couldn't get my head around which way round j and n were supposed to go - plus I wanted to see if I'd get the right plot by switching them. J does have to go to N though and n goes to M for sure. I've swapped fnj(j,n) for fnj, n,j). I'll explain a bit more about this in an answer.
Dr. Seis
Dr. Seis on 26 Jan 2012
If you swap fnj(j,n) for fnj(n,j) - make sure you get both -, then you don't need to make the change to sum. Leave it as "Fn = sum(fnj);"

Sign in to comment.

More Answers (2)

bym
bym on 26 Jan 2012
You do not have 2 lines on the graph: it goes from left to right & back again. After some noodling, I think you want to replace
Fn = sum(fnj);
%plot(AnRec,Fn)
pn = Fn*(M./(ar*N));
plot(xjRec,pn,'LineWidth',1.7)
with the following, which makes more sense (or at the very least presents more pleasing plots :) )
Fn = sum(fnj,2);
plot(AnRec,Fn)
figure
pn = Fn*(M./(ar*N));
plot(xjRec,pn,'LineWidth',1.7)
  2 Comments
Tom
Tom on 26 Jan 2012
Thanks - so what does the added ',2' do?
Dr. Seis
Dr. Seis on 26 Jan 2012
Adding the 2 will sum all the values in a row.

Sign in to comment.


Tom
Tom on 26 Jan 2012
I have this as a similar problem, which I can use as a kind of template for this one...
'1024 samples of one period of a 1Hz sine wave are taken. By analysing the amplitude range from – 10 Pa to 10 Pa in 60 equal segments, estimate the probability density function of the sine wave.'
The solution is here - http://dl.dropbox.com/u/11341635/Tut10-11SolQ2.png No built in Matlab functions allowed!
It's basically sampling the old-fashioned way - creating an array of static amplitude points as a framework for the receiving part of the system. The incoming signal is sampled at a given sampling frequency, and each sample received by the system will be analysed to see which amplitude point to match it up to - hence the two FOR loops, with an IF in the middle. One FOR loop is to consider each sample one at a time, and the other FOR loop is to determine whether each incoming sample in close to the amplitude point at hand. If it is, a value of 1 is placed in the 2-D matrix, the size of which is determined by the no. of amplitude points and the sampling frequency.
Phew.

Categories

Find more on MATLAB in Help Center and File Exchange

Tags

Products

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!