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

How to make this signal linear

If you want to eliminate the low-frequency parabolic(?) trend, the easiest way is to use a digital filter: signal_2 = load('si...

How to make this signal linear

If you want to eliminate the low-frequency parabolic(?) trend, the easiest way is to use a digital filter: signal_2 = load('si...

alrededor de 19 horas ago | 0

Answered

How to solve 3 non-linear equations with 3 unknowns?

The ‘x0’ in the solve call is not appropriate. The solve function wants name-value pairs after the first two arguments, and ‘x0...

How to solve 3 non-linear equations with 3 unknowns?

The ‘x0’ in the solve call is not appropriate. The solve function wants name-value pairs after the first two arguments, and ‘x0...

alrededor de 22 horas ago | 0

Answered

I'm trying to understand your tables command bugs

Both of these work for me: A = randi([0 20],5, 5); Atable1 = array2table(A, 'RowNames',{'a' 'b' 'c' 'd' 'e'},'VariableNames',...

I'm trying to understand your tables command bugs

Both of these work for me: A = randi([0 20],5, 5); Atable1 = array2table(A, 'RowNames',{'a' 'b' 'c' 'd' 'e'},'VariableNames',...

1 día ago | 0

Answered

Designing IIR Butterworth filter (pass band)

Yes. Just define the order as 5. There are different design functions in the Signal Processing Toolbox (butter that I have a...

Designing IIR Butterworth filter (pass band)

Yes. Just define the order as 5. There are different design functions in the Signal Processing Toolbox (butter that I have a...

1 día ago | 0

Answered

Play 2 wav files after each other

Use the playblocking function to get them to play the way you want them. [y1, Fs1] = audioread('file1.wav'); %this file i...

Play 2 wav files after each other

Use the playblocking function to get them to play the way you want them. [y1, Fs1] = audioread('file1.wav'); %this file i...

2 días ago | 1

Answered

I have a function which is ranged from 0 microseconds to 50 microseconds. I just want to plot its exact copies from 0to50 microsecond,50 to 100microseconds,100 to 150 and so on

Try this: Ts=1/200e6; t=Ts:Ts:0.00005; y=[ones(1,2000) zeros(1,8000)]; x=sin(1256e6*t + 2198e10*t.^2); z=x.^y; NR = ...

I have a function which is ranged from 0 microseconds to 50 microseconds. I just want to plot its exact copies from 0to50 microsecond,50 to 100microseconds,100 to 150 and so on

Try this: Ts=1/200e6; t=Ts:Ts:0.00005; y=[ones(1,2000) zeros(1,8000)]; x=sin(1256e6*t + 2198e10*t.^2); z=x.^y; NR = ...

2 días ago | 0

Answered

How to condense a plot with a shorter x-axis?

Divide the x-axis vector by max(x), assuming ‘x’ is the x-axis vector. That will re-scale it without losing any information.

How to condense a plot with a shorter x-axis?

Divide the x-axis vector by max(x), assuming ‘x’ is the x-axis vector. That will re-scale it without losing any information.

2 días ago | 0

Answered

How to avoid a for loop in functions?

Use ‘logical indexing’ and anonymous functions to eliminate the for loop and the if block: PsolWB1 = @(press) (973.-(70400./(1...

How to avoid a for loop in functions?

Use ‘logical indexing’ and anonymous functions to eliminate the for loop and the if block: PsolWB1 = @(press) (973.-(70400./(1...

3 días ago | 0

| accepted

Answered

fsolve for non linear equation

One problem is that you called fsolve inside the function. That is not appropriate. The other is that you use global variables...

fsolve for non linear equation

One problem is that you called fsolve inside the function. That is not appropriate. The other is that you use global variables...

3 días ago | 0

Answered

power law fit to find exponent

I would do this (no Toolboxes required): x = linspace(0, 5, 50); % Create Data y = x.^pi ...

power law fit to find exponent

I would do this (no Toolboxes required): x = linspace(0, 5, 50); % Create Data y = x.^pi ...

3 días ago | 0

Answered

Definite integral without symbolic variables

Create ‘fun’ as an anonymous function, then use the integral function: fun = @(z) (G1 + G2 - 2*cos(k*L1/2) * G0) .* sin(k*(L2/...

Definite integral without symbolic variables

Create ‘fun’ as an anonymous function, then use the integral function: fun = @(z) (G1 + G2 - 2*cos(k*L1/2) * G0) .* sin(k*(L2/...

3 días ago | 0

Answered

kinetic energy of the wind (u and v vectors)

Apparently, ‘h’, ‘u’, and ‘v’ are vectors of the same dimensions. One approach is to code the function as an anonymous functi...

kinetic energy of the wind (u and v vectors)

Apparently, ‘h’, ‘u’, and ‘v’ are vectors of the same dimensions. One approach is to code the function as an anonymous functi...

3 días ago | 0

| accepted

Answered

how to plot euler formula in matlab?

The plot function plots lines between points. You have only one point, so if you change the plot call to: y = exp(100*i*pi) ...

how to plot euler formula in matlab?

The plot function plots lines between points. You have only one point, so if you change the plot call to: y = exp(100*i*pi) ...

3 días ago | 0

Answered

Please help me understand the use of dot operator

The dot operator, used with multiplication, division, and exponentiation, creates element-wise oiperations. See Array vs. Matri...

Please help me understand the use of dot operator

The dot operator, used with multiplication, division, and exponentiation, creates element-wise oiperations. See Array vs. Matri...

3 días ago | 1

| accepted

Answered

How to sort the rows of an array according to another vector?

Try this: D=[1 1 0 1 0 1; 4 6 7 8 9 9; 1 1 1 2 3 4]; b=[2 44 0]; [bb,idx] = sort(b,'descend'); Out = D(idx,:) producing: ...

How to sort the rows of an array according to another vector?

Try this: D=[1 1 0 1 0 1; 4 6 7 8 9 9; 1 1 1 2 3 4]; b=[2 44 0]; [bb,idx] = sort(b,'descend'); Out = D(idx,:) producing: ...

3 días ago | 0

| accepted

Answered

How to change the values in a matrix

One approach: Before = [ 1 9 9 9 3 8 9 9 9 4 8 9 9 9 ...

How to change the values in a matrix

One approach: Before = [ 1 9 9 9 3 8 9 9 9 4 8 9 9 9 ...

3 días ago | 0

| accepted

Answered

Inverse of algebraic expression in Matlab

Possibly: rfcn = @(E) fsolve(@(r) 1-1/2*((1+r).*log2(1+r)+(1-r).*log2(1-r)) - E, 10); % Anonymous Function: r(E) N = 200;...

Inverse of algebraic expression in Matlab

Possibly: rfcn = @(E) fsolve(@(r) 1-1/2*((1+r).*log2(1+r)+(1-r).*log2(1-r)) - E, 10); % Anonymous Function: r(E) N = 200;...

4 días ago | 0

| accepted

Answered

Designing an Arbitrary FIR Filter

This is probably as close as you can get: Frequency = [100 125 160 200 250 315 400 500 630 800 1000 1250 1600 200...

Designing an Arbitrary FIR Filter

This is probably as close as you can get: Frequency = [100 125 160 200 250 315 400 500 630 800 1000 1250 1600 200...

4 días ago | 1

| accepted

Answered

Plot the Intensity and phase of variation of electric field.

The Fourier transform is easy enough to calculate and plot: t = dlmread('t_spectrum_11480.txt','',1,0); E = dlmread('y_spectr...

Plot the Intensity and phase of variation of electric field.

The Fourier transform is easy enough to calculate and plot: t = dlmread('t_spectrum_11480.txt','',1,0); E = dlmread('y_spectr...

4 días ago | 0

| accepted

Answered

How to fix axes in dynamic plot

Try this: x=@(t)(sin(4*t)).^2; y=@(t)cos(5*t); z=pi; figure(1) for t=0:0.002:z plot(x(t),y(t),'.','MarkerSize',4,'c...

How to fix axes in dynamic plot

Try this: x=@(t)(sin(4*t)).^2; y=@(t)cos(5*t); z=pi; figure(1) for t=0:0.002:z plot(x(t),y(t),'.','MarkerSize',4,'c...

5 días ago | 0

| accepted

Answered

How can plot this curve in matlab ?

This appears to approximate it: k = linspace(0,1); Finesse = 5./k; figure plot(k, Finesse)

How can plot this curve in matlab ?

This appears to approximate it: k = linspace(0,1); Finesse = 5./k; figure plot(k, Finesse)

5 días ago | 0

| accepted

Answered

Extraction of FFT from each segment for computation

If this relates to your earlier Question, the fft results are in each column of the ’P1’ matrix. Y = fft(X-mean(X)); ...

Extraction of FFT from each segment for computation

If this relates to your earlier Question, the fft results are in each column of the ’P1’ matrix. Y = fft(X-mean(X)); ...

5 días ago | 0

Answered

How to integrate a data set?

I would use the trapz function, or cumtrapz (linked to in the trapz documentation), depending on the result you want. For exa...

How to integrate a data set?

I would use the trapz function, or cumtrapz (linked to in the trapz documentation), depending on the result you want. For exa...

5 días ago | 0

| accepted

Answered

Equations of 3d Curves

If you want to differentiate them using the Symbolic Math Toolbox, start with: syms t x(t) = 5*t*sin(4*t) y(t) = 12*t + 200...

Equations of 3d Curves

If you want to differentiate them using the Symbolic Math Toolbox, start with: syms t x(t) = 5*t*sin(4*t) y(t) = 12*t + 200...

5 días ago | 0

Answered

Why won't my x-axis tick label format change?

Add these two lines to your code: xt_sigfig = round(xt, 3, 'significant'); ax.XTickLabel = xt_sigfig; so the complet...

Why won't my x-axis tick label format change?

Add these two lines to your code: xt_sigfig = round(xt, 3, 'significant'); ax.XTickLabel = xt_sigfig; so the complet...

5 días ago | 0

| accepted

Answered

how characters connected in quotes

See the char documentation for information on character arrays, and string for information on string arrays. Those should pro...

how characters connected in quotes

See the char documentation for information on character arrays, and string for information on string arrays. Those should pro...

6 días ago | 0

| accepted

Answered

Design IIR filter from absorption data

Use the freqz function: figure freqz(b, a, 2^14, Fs) set(subplot(2,1,1), 'XLim',[0 1E+3]) % Optional set(sub...

Design IIR filter from absorption data

Use the freqz function: figure freqz(b, a, 2^14, Fs) set(subplot(2,1,1), 'XLim',[0 1E+3]) % Optional set(sub...

6 días ago | 0

Answered

Split every current signal into N segments of Length L

Assuming I understand waht you want to do, I would use the Signal Processing Toolbox buffer function to do that. You would ha...

Split every current signal into N segments of Length L

Assuming I understand waht you want to do, I would use the Signal Processing Toolbox buffer function to do that. You would ha...

6 días ago | 0

| accepted

Answered

problem in z transform

Use the simplify function: ZT = simplify(ZT, 'Steps',500); to get: ZT = 2/(z - 2)

problem in z transform

Use the simplify function: ZT = simplify(ZT, 'Steps',500); to get: ZT = 2/(z - 2)

6 días ago | 0

| accepted

Answered

How to use fzero to solve for a variable in an exponent?

Take the fzero call out of the function it calls: X1 = 0.5; X2 = 0.5; A12 = 1.6798; A21 = 0.9227; AW = 7.96681; BW = 1668...

How to use fzero to solve for a variable in an exponent?

Take the fzero call out of the function it calls: X1 = 0.5; X2 = 0.5; A12 = 1.6798; A21 = 0.9227; AW = 7.96681; BW = 1668...

6 días ago | 0