12.851 total contributions since 2012

Hic sunt dracones!

PROFESSIONAL: Physician (M.D.): Diplomate of the American Board of Internal Medicine; M.Sc. Biomedical Engineering: Instrumentation, Signal Processing, Control, System Identification, Parameter Estimation

NON-PROFESSIONAL: Amateur Extra Class Amateur Radio Operator; Private Pilot, Airplane Single Engine Land, Instrument Rating Airplane; Gamer

NOTE: I do not respond to emails or personal messages, unless they are about my File Exchange contributions.

Answered

Error trying to plot circle

You apparently want to plot three different circles with different radii. This works: A = [10 8 2; 10 8 1; -2 9 1; 2 -1 1] ...

Error trying to plot circle

You apparently want to plot three different circles with different radii. This works: A = [10 8 2; 10 8 1; -2 9 1; 2 -1 1] ...

ongeveer een uur ago | 0

| accepted

Answered

Marking a specific point in a graph made out of a loop

Try this: b=400; %mm d=500; %mm Asv = [3000 5000 3000 5000]; %mm^2 fckv = [30 30 90 90]; %Mpa tic for k = 1:numel(fc...

Marking a specific point in a graph made out of a loop

Try this: b=400; %mm d=500; %mm Asv = [3000 5000 3000 5000]; %mm^2 fckv = [30 30 90 90]; %Mpa tic for k = 1:numel(fc...

ongeveer 19 uur ago | 1

Answered

Calculating the area under a curve using cell arrays

Try this: D = load('pwp_diff.mat'); dPWP = D.PWP_diff; tv = 0:size(dPWP,2)-1; tvi = linspace(min(tv), max(tv), numel(tv)*10...

Calculating the area under a curve using cell arrays

Try this: D = load('pwp_diff.mat'); dPWP = D.PWP_diff; tv = 0:size(dPWP,2)-1; tvi = linspace(min(tv), max(tv), numel(tv)*10...

1 dag ago | 0

| accepted

Answered

How to retrieve first 4 digits after decimal point?

One approach: x = pi AfterDecimal = rem(x,1) FirstFour = fix(AfterDecimal*1E+4)*1E-4 producing: x = 3.141592653589793...

How to retrieve first 4 digits after decimal point?

One approach: x = pi AfterDecimal = rem(x,1) FirstFour = fix(AfterDecimal*1E+4)*1E-4 producing: x = 3.141592653589793...

1 dag ago | 0

| accepted

Answered

To arrange matrix based off criteria in column

Try this: Col1 = randperm(300); % Create Column 1 Col2 = randi(9, 1, 300); ...

To arrange matrix based off criteria in column

Try this: Col1 = randperm(300); % Create Column 1 Col2 = randi(9, 1, 300); ...

1 dag ago | 0

| accepted

Answered

Least squares fit/line fit for 3D data

For a linear regression, this is straightforward: B = [x(:) y(:) ones(size(x(:)))] \ z(:); % Linear Paramet...

Least squares fit/line fit for 3D data

For a linear regression, this is straightforward: B = [x(:) y(:) ones(size(x(:)))] \ z(:); % Linear Paramet...

1 dag ago | 0

Answered

How to find the time it takes till the first upward deflection of square wave pulse?

See if the midcross function willl do what you want.

How to find the time it takes till the first upward deflection of square wave pulse?

See if the midcross function willl do what you want.

2 dagen ago | 0

Answered

How to plot this?

I am not certain what you want to plot, since you did not specify that. Try this: tideTime = linspace(0, 4*pi); ...

How to plot this?

I am not certain what you want to plot, since you did not specify that. Try this: tideTime = linspace(0, 4*pi); ...

2 dagen ago | 0

Answered

pzmap(G+G) produces incorrect plot

The ‘+’ operator connects the two ‘G’ models in parallel. They do appear to have pole-zero cancellation as the result: s = tf...

pzmap(G+G) produces incorrect plot

The ‘+’ operator connects the two ‘G’ models in parallel. They do appear to have pole-zero cancellation as the result: s = tf...

2 dagen ago | 0

| accepted

Answered

Dividing the values of a vector

Try this: A=[1 1.1; 1.1 1]; l = [2; 1]/3; % Create As Column Vector w=1; B=[1.09 ...

Dividing the values of a vector

Try this: A=[1 1.1; 1.1 1]; l = [2; 1]/3; % Create As Column Vector w=1; B=[1.09 ...

2 dagen ago | 0

Answered

Creating a Code to determine Lift Curve Slope

There are several problems with your code. You are only saving the last iteration of ‘A’, and since it does not appear that ‘...

Creating a Code to determine Lift Curve Slope

There are several problems with your code. You are only saving the last iteration of ‘A’, and since it does not appear that ‘...

2 dagen ago | 0

| accepted

Answered

Parallel form filter design

I am not quite certain what you want to do. To express your elliptic filter as second-order-sections, the Signal Processing T...

Parallel form filter design

I am not quite certain what you want to do. To express your elliptic filter as second-order-sections, the Signal Processing T...

3 dagen ago | 0

Answered

Getting a for loop to save a plot every step

Use the same equation and pass the different parameters as extra parameters in a loop. See the documentation secton on Passing ...

Getting a for loop to save a plot every step

Use the same equation and pass the different parameters as extra parameters in a loop. See the documentation secton on Passing ...

3 dagen ago | 0

Answered

how to plot a conventional DSB AM signal

I am not certain how much help you need for this. For a start, see the documentation for the plot function. I am not certain i...

how to plot a conventional DSB AM signal

I am not certain how much help you need for this. For a start, see the documentation for the plot function. I am not certain i...

3 dagen ago | 0

Answered

How to formulate and use non linier curve-fitting

Try this: D = dlmread('Example.txt', '\t', 1, 0); objfcn = @(b,x) exp(b(1) - b(2).*x(:,1).*x(:,2)./x(:,3)); ...

How to formulate and use non linier curve-fitting

Try this: D = dlmread('Example.txt', '\t', 1, 0); objfcn = @(b,x) exp(b(1) - b(2).*x(:,1).*x(:,2)./x(:,3)); ...

3 dagen ago | 0

| accepted

Answered

Solving for a system

The easiest way is to loop through the values of ‘ax’: lambda=1.55e-6; %Field in vertical direction axv=linspace(1e-6,2e-6,1...

Solving for a system

The easiest way is to loop through the values of ‘ax’: lambda=1.55e-6; %Field in vertical direction axv=linspace(1e-6,2e-6,1...

3 dagen ago | 0

| accepted

Answered

How can I show (x,y) of a point in fsurf function?

Try this: Rx = eye(2); Rdx = [2; 4.5]; wo = Rx\Rdx; var_d = 24.4; syms w1 w2 vec_w = [w1 ; w2]; J(w1,w2) = var_d - 2*vec...

How can I show (x,y) of a point in fsurf function?

Try this: Rx = eye(2); Rdx = [2; 4.5]; wo = Rx\Rdx; var_d = 24.4; syms w1 w2 vec_w = [w1 ; w2]; J(w1,w2) = var_d - 2*vec...

4 dagen ago | 1

| accepted

Answered

how to find kurtosis in freqeuncy domain

Your best option is likey to use a digital bandpass filter (for 0-2 Hz use a lowpass filter) to isolate those frequency ranges i...

how to find kurtosis in freqeuncy domain

Your best option is likey to use a digital bandpass filter (for 0-2 Hz use a lowpass filter) to isolate those frequency ranges i...

4 dagen ago | 0

Answered

Solving an Equation and then plotting a graph

In these two lines: sigmaSteel=@(c) Es*epss(c) .* (epss(c)<=epsy) + fy .* (epss(c)>epsy & epss(c)<=epssh) + (fsu+(fy-fsu)*abs(...

Solving an Equation and then plotting a graph

In these two lines: sigmaSteel=@(c) Es*epss(c) .* (epss(c)<=epsy) + fy .* (epss(c)>epsy & epss(c)<=epssh) + (fsu+(fy-fsu)*abs(...

4 dagen ago | 1

| accepted

Answered

Selecting numbers in matrix that are repeated

Try this: A = randi(1000, 60, 5); % Create Matrix [Au,~,ix] = unique(A,'stable'); ...

Selecting numbers in matrix that are repeated

Try this: A = randi(1000, 60, 5); % Create Matrix [Au,~,ix] = unique(A,'stable'); ...

4 dagen ago | 0

| accepted

Answered

Passing parameters in boundary value problem using BVP4C

Define ‘a’ in your workspace, then use this bvp4c call: sol = bvp4c(@(x,y)twoode(x,y,a),@(ya,yb)twobc(ya,yb,a),solinit); You ...

Passing parameters in boundary value problem using BVP4C

Define ‘a’ in your workspace, then use this bvp4c call: sol = bvp4c(@(x,y)twoode(x,y,a),@(ya,yb)twobc(ya,yb,a),solinit); You ...

4 dagen ago | 1

| accepted

Answered

Unit step function explanation

Use an anonymous function to shorten it even further: u = @(t) t>0; It works by using ‘logical indexing’, so that here, if ‘t...

Unit step function explanation

Use an anonymous function to shorten it even further: u = @(t) t>0; It works by using ‘logical indexing’, so that here, if ‘t...

5 dagen ago | 1

| accepted

Answered

How to plot a surface of vector z in a 2D grid made of vectors x and y?

I am not certain what you want the plots to look like, since I am not at all familiar with what you are doing. Try this: fi...

How to plot a surface of vector z in a 2D grid made of vectors x and y?

I am not certain what you want the plots to look like, since I am not at all familiar with what you are doing. Try this: fi...

5 dagen ago | 0

Answered

I cant figure out how to make the plot show up

Try this: stage1BurnTimeStart=0; stage2BurnTimeStart=165; stage3BurnTimeStart=360; stage3BurnTimeEnd=500; tv = 6.5:stage3B...

I cant figure out how to make the plot show up

Try this: stage1BurnTimeStart=0; stage2BurnTimeStart=165; stage3BurnTimeStart=360; stage3BurnTimeEnd=500; tv = 6.5:stage3B...

5 dagen ago | 0

| accepted

Answered

How to index one array with the index of another array that meets a certain condition

Try this: cloud_type = [0 0 6 8 8 5 1 8 2 2 8]; radius = [20 13 14 25 30 10 5 27 13 14 25]; out = radius(cloud_type==8) p...

How to index one array with the index of another array that meets a certain condition

Try this: cloud_type = [0 0 6 8 8 5 1 8 2 2 8]; radius = [20 13 14 25 30 10 5 27 13 14 25]; out = radius(cloud_type==8) p...

5 dagen ago | 1

| accepted

Answered

Error :Nonscalar arrays of function handles are not allowed; use cell arrays instead.

The loop is not necessary. This required a few corrections (a typographical error and to vectorise the exponentiation), and n...

Error :Nonscalar arrays of function handles are not allowed; use cell arrays instead.

The loop is not necessary. This required a few corrections (a typographical error and to vectorise the exponentiation), and n...

5 dagen ago | 0

| accepted

Answered

How to write Continue Long Statements on Multiple Lines ?

Use repmat to avoid manually typing all 160 of those ‘%q’ repeats! Try this: formatSpec = ['%f',repmat('%q',1,160), '%[^\n\...

How to write Continue Long Statements on Multiple Lines ?

Use repmat to avoid manually typing all 160 of those ‘%q’ repeats! Try this: formatSpec = ['%f',repmat('%q',1,160), '%[^\n\...

5 dagen ago | 0

| accepted

Answered

Positive or zero intercept for linear function

Linearising the exponiential function may not be the best option. I am not certain what you are modeling, however using a nonli...

Positive or zero intercept for linear function

Linearising the exponiential function may not be the best option. I am not certain what you are modeling, however using a nonli...

5 dagen ago | 0

Answered

for loop with different range of number

You can use a for loop, however I would do this instead: r = 1:500; % Create ‘r’ To T...

for loop with different range of number

You can use a for loop, however I would do this instead: r = 1:500; % Create ‘r’ To T...

6 dagen ago | 0

Answered

How do integrate this function over t=0 to t=165?

I am not certain what you want to integrate, or what you want the final result to be. Try this: mi=2290000; mf=130000; tb...

How do integrate this function over t=0 to t=165?

I am not certain what you want to integrate, or what you want the final result to be. Try this: mi=2290000; mf=130000; tb...

6 dagen ago | 0

| accepted