Conducting Radar Maritime Reflectivity

% Analyze for sea states 1 through 5 ss = 1:5; % Sea states

% Initialize outputs numSeaStates = numel(ss); hgtsd = zeros(1,numSeaStates); beta0 = zeros(1,numSeaStates); vw= zeros(1,numSeaStates);

% Obtain sea state properties for is = 1:numSeaStates [hgtsd(is),beta0(is),vw(is)] = searoughness(ss(is)); end

% Plot results helperPlotSeaRoughness(ss,hgtsd,beta0,vw);

[Image Analyst]

I formatted your code as code for you, but please define "Coducting" - I'm not familiar with that word.

If you have any more questions, then attach your data and code to read it in with the paperclip icon after you read this:

TUTORIAL: How to ask a question (on Answers) and get a fast answer

Especially ask a question, and give us the code for searoughness and helperPlotSeaRoughness.

[John D'Errico]

What is your question?

[I.Osman]

The code doesn't working

[Dyuman Joshi]

How is it not working? What are you trying to do?

Do you have the Radar toolbox that is required to run searoughness? What is helperPlotSeaRoughness()?

You have to specify things, as the mind reading toolbox has not been released yet.

Go through the link @Image Analyst has mentioned in their comment and edit your question accordingly.

Answer 1

[Sara James]

First, the above code requires Radar Toolbox to use the searoughness function.

The helperPlotSeaRoughness is not a documented Radar Toolbox function. It was a helper specifically created for the Maritime Radar Sea Clutter Modeling example. If you would like to use this function elsewhere, you can save off the local function from the example livescript into its own file. The code for the helper from that example is below.

function helperPlotSeaRoughness(ss,hgtsd,beta0,vw)
% Creates 3x1 plot of sea roughness outputs 

% Create figure
figure

% Plot standard deviation of sea wave height
subplot(3,1,1)
plot(ss,hgtsd,'-o','LineWidth',1.5)
ylabel([sprintf('Wave\nHeight ') '\sigma_h (m)'])
title('Sea Wave Roughness')
grid on; 

% Plot sea wave slope
subplot(3,1,2)
plot(ss,beta0,'-o','LineWidth',1.5)
ylabel([sprintf('Wave\nSlope ') '\beta_0 (deg)'])
grid on; 

% Plot wind velocity 
subplot(3,1,3)
plot(ss,vw,'-o','LineWidth',1.5)
xlabel('Sea State')
ylabel([sprintf('Wind\nVelocity ')  'v_w (m/s)'])
grid on; 
end