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

Obtaining line equation having multiple points

The polyfit (and polyval) functions may be what you want. There are other options, however these are likely the easiest to us...

Obtaining line equation having multiple points

The polyfit (and polyval) functions may be what you want. There are other options, however these are likely the easiest to us...

ungefär 16 timmar ago | 0

Answered

how to fix the error ''Conversion to function_handle from double is not possible.'' for the below 3Eq euler matlab code

First, name the variables slightly differently from the function names, second define the initial conditions to conform with tho...

how to fix the error ''Conversion to function_handle from double is not possible.'' for the below 3Eq euler matlab code

First, name the variables slightly differently from the function names, second define the initial conditions to conform with tho...

ungefär 19 timmar ago | 1

Answered

how to color in between curves the two graphs?

It’s not possible to run images of code, especially images of incomplete code. Try this example: x1 = linspace(-0.1, -10, 3...

how to color in between curves the two graphs?

It’s not possible to run images of code, especially images of incomplete code. Try this example: x1 = linspace(-0.1, -10, 3...

ungefär 20 timmar ago | 0

Answered

Show numbers with Grouped Bars

Try this: x = [98 96; 142 52; 42 42]; hBar = bar(x); for k1 = 1:size(x,2) ctr(k1,:) = bsxfun(@plus, 1:numel(hBar(k1).XD...

Show numbers with Grouped Bars

Try this: x = [98 96; 142 52; 42 42]; hBar = bar(x); for k1 = 1:size(x,2) ctr(k1,:) = bsxfun(@plus, 1:numel(hBar(k1).XD...

en dag ago | 0

| accepted

Answered

How to fit this curve?

The best model is the mathematical expression of the process that created those data. A power law fit appears to provide a re...

How to fit this curve?

The best model is the mathematical expression of the process that created those data. A power law fit appears to provide a re...

en dag ago | 1

| accepted

Answered

Import text files with extensions other than .txt

Use detectImportOptions to tell readtable to read the file as a text file. One recent example from an earlier Answer: opts ...

Import text files with extensions other than .txt

Use detectImportOptions to tell readtable to read the file as a text file. One recent example from an earlier Answer: opts ...

en dag ago | 1

| accepted

Answered

± use in MATLAB

One option: CL = [(0.331+0.008) + (0.15+0.002)*ALPHA (0.331-0.008) + (0.15-0.002)*ALPHA]; Experiment to get differen...

± use in MATLAB

One option: CL = [(0.331+0.008) + (0.15+0.002)*ALPHA (0.331-0.008) + (0.15-0.002)*ALPHA]; Experiment to get differen...

en dag ago | 0

Answered

when i plot my figure doesnt show line

Use a colon to create the vectors instead of a semicolon, and be sure the vectors are equal length: x = 1:5; y = 5:9; plot...

when i plot my figure doesnt show line

Use a colon to create the vectors instead of a semicolon, and be sure the vectors are equal length: x = 1:5; y = 5:9; plot...

en dag ago | 0

Answered

How to create a Symbolic vector

Try this:: n = inputdlg('Enter the length of the vector (integer): '); n = str2double(n); syms x x = sym('x',[1 n]) pro...

How to create a Symbolic vector

Try this:: n = inputdlg('Enter the length of the vector (integer): '); n = str2double(n); syms x x = sym('x',[1 n]) pro...

2 dagar ago | 0

Answered

1D interpolation by y = ax+b

To set the start point at the origin (0,0), calculate the slope as: slope = x(:) \ y(:); for example: x = 0:0.1:10; y = ra...

1D interpolation by y = ax+b

To set the start point at the origin (0,0), calculate the slope as: slope = x(:) \ y(:); for example: x = 0:0.1:10; y = ra...

2 dagar ago | 0

Answered

A little question of [] and for loop

The full code should actually be: A = []; for i=1:50 A=[A i]; end It works by concatenating the value of ‘i’ to existi...

A little question of [] and for loop

The full code should actually be: A = []; for i=1:50 A=[A i]; end It works by concatenating the value of ‘i’ to existi...

2 dagar ago | 1

Answered

Plot surface from non-anonymus function

See if this does what you want — function z = test(w) x = w{1}; y = w{2}; z = 2*x.^2 + 2*y.^2 +10*(1.1*y.^2 +...

Plot surface from non-anonymus function

See if this does what you want — function z = test(w) x = w{1}; y = w{2}; z = 2*x.^2 + 2*y.^2 +10*(1.1*y.^2 +...

2 dagar ago | 0

| accepted

Answered

Find distance between one point and a subset of other points.

See if the pdist function will do what you want. (If so, squareform may also be useful.)

Find distance between one point and a subset of other points.

See if the pdist function will do what you want. (If so, squareform may also be useful.)

3 dagar ago | 0

Answered

detrending timeseries by removing the linear function fitted by samples?

Considering that the objective is to remove a linear trend, perhaps the easiest way would be: t = 11:22; p = polyfit(x(t),tac...

detrending timeseries by removing the linear function fitted by samples?

Considering that the objective is to remove a linear trend, perhaps the easiest way would be: t = 11:22; p = polyfit(x(t),tac...

3 dagar ago | 0

| accepted

Answered

Find peak and plot graph from peak

I am not certain what you want to do. If you want to get the peaks and their indices, request those outputs: [pks,locs] = f...

Find peak and plot graph from peak

I am not certain what you want to do. If you want to get the peaks and their indices, request those outputs: [pks,locs] = f...

3 dagar ago | 0

| accepted

Answered

How can I define a gradient function that can update k values after each iteration to optimize parameters in ODEs?

One option is the approach in Parameter Estimation for a System of Differential Equations. You will need to modify the ‘kinetic...

How can I define a gradient function that can update k values after each iteration to optimize parameters in ODEs?

One option is the approach in Parameter Estimation for a System of Differential Equations. You will need to modify the ‘kinetic...

3 dagar ago | 0

Answered

How to draw bar graph of different colors ?

Try this: y = [280 260 240]; figure c = bar(y); set(gca, 'XTickLabel',{'CS0','CS10','CS20'}) xlabel('Wt% of cenospheres'...

How to draw bar graph of different colors ?

Try this: y = [280 260 240]; figure c = bar(y); set(gca, 'XTickLabel',{'CS0','CS10','CS20'}) xlabel('Wt% of cenospheres'...

3 dagar ago | 1

| accepted

Answered

Problem to get the right format of a Surf Plot

It looks as though ‘X’ is not sorted. Without the data a guess would be this possible solution: BB = sortrows (BB,1); See th...

Problem to get the right format of a Surf Plot

It looks as though ‘X’ is not sorted. Without the data a guess would be this possible solution: BB = sortrows (BB,1); See th...

3 dagar ago | 0

| accepted

Answered

How to linearly fit a semilog graph

Try this: D1 = readmatrix('Data.xlsx'); x = D1(:,1); y = D1(:,2); [ymax,idx] = max(y); idxrng = idx:numel(y); objfcn ...

How to linearly fit a semilog graph

Try this: D1 = readmatrix('Data.xlsx'); x = D1(:,1); y = D1(:,2); [ymax,idx] = max(y); idxrng = idx:numel(y); objfcn ...

3 dagar ago | 0

Answered

I don't understand why the elements are not the same on both sides

The reason is that ‘t’ needs a subscript: x(n+1) = x(n) + h * ((8*(x(n))/(t(n)+1))-(t+1)*(x(n)^2)); ...

I don't understand why the elements are not the same on both sides

The reason is that ‘t’ needs a subscript: x(n+1) = x(n) + h * ((8*(x(n))/(t(n)+1))-(t+1)*(x(n)^2)); ...

3 dagar ago | 0

| accepted

Answered

Adding a new column to a table

One approach: T1 = table(randi([330 470],10,1), 'VariableNames',{'QTcorrected'}); % Create Data QTc = @(x) 1*...

Adding a new column to a table

One approach: T1 = table(randi([330 470],10,1), 'VariableNames',{'QTcorrected'}); % Create Data QTc = @(x) 1*...

4 dagar ago | 0

| accepted

Answered

How to Segment EMG data into time intervals

One option is the Signal Processing Toolbox buffer function. It uses samples as one of its arguments, however that is straigh...

How to Segment EMG data into time intervals

One option is the Signal Processing Toolbox buffer function. It uses samples as one of its arguments, however that is straigh...

4 dagar ago | 0

| accepted

Answered

Coupled rate ODEs with ode45

The concentrations change appropriately, however they don’t change much and the concentrations are vanishingly small. That’s th...

Coupled rate ODEs with ode45

The concentrations change appropriately, however they don’t change much and the concentrations are vanishingly small. That’s th...

5 dagar ago | 0

| accepted

Answered

Figure line won't connect/display all data points

Isolated points only plot if specified as markers. Try this: D1 = load('midpoint.mat'); midpoint = D1.midpoint_overlap_seq...

Figure line won't connect/display all data points

Isolated points only plot if specified as markers. Try this: D1 = load('midpoint.mat'); midpoint = D1.midpoint_overlap_seq...

5 dagar ago | 0

| accepted

Answered

How to solve this equation?

Probably the easiest approach: syms h z2 Eqn = h^0.66+0.0805*h == 0.045; [h,prms,cndx] = solve(Eqn,h, 'ReturnConditions',1)...

How to solve this equation?

Probably the easiest approach: syms h z2 Eqn = h^0.66+0.0805*h == 0.045; [h,prms,cndx] = solve(Eqn,h, 'ReturnConditions',1)...

5 dagar ago | 0

Answered

ode45: use span of length 2

The tspan argument can be anything you want it to be (within limits). To have ode45 to evaluate and output at only those valu...

ode45: use span of length 2

The tspan argument can be anything you want it to be (within limits). To have ode45 to evaluate and output at only those valu...

5 dagar ago | 0

Answered

confusing error message: 'Not enough input arguments'?

With a function in that format, it is necessary to pass a function handle to lsqnonlin (or any other function that takes a funct...

confusing error message: 'Not enough input arguments'?

With a function in that format, it is necessary to pass a function handle to lsqnonlin (or any other function that takes a funct...

5 dagar ago | 0

| accepted

Answered

How to plot only the real solutions of an implicit function ?

Nothing is being generated: R = 0.35; L = 0.25; H = 0.4; f1 = @(x,y) sqrt(L^2-(R*(cosd(x)-sind(x).*sind(y))-H).^2); figu...

How to plot only the real solutions of an implicit function ?

Nothing is being generated: R = 0.35; L = 0.25; H = 0.4; f1 = @(x,y) sqrt(L^2-(R*(cosd(x)-sind(x).*sind(y))-H).^2); figu...

6 dagar ago | 1

Answered

Plot and calculate Damping ratio

The findpeaks or islocalmax functions can return the information to do the calculations.

Plot and calculate Damping ratio

The findpeaks or islocalmax functions can return the information to do the calculations.

6 dagar ago | 1

Answered

Trouble taking the differential of a function

‘What am I missing here?’ Most likely, using the Symbolic Math Toolbox. Try this: syms E I L w_0 x func(x) = (((w_0*L)...

Trouble taking the differential of a function

‘What am I missing here?’ Most likely, using the Symbolic Math Toolbox. Try this: syms E I L w_0 x func(x) = (((w_0*L)...

6 dagar ago | 0