Community Profile

David Wilson

Auckland University of Technology

96 total contributions since 2011

View details...

Contributions in
View by

creating a transfer function from a determinant
How about: [num,den] = numden(G_vel) G_vel_tf = tf(sym2poly(num), sym2poly(den))

7 jours ago | 0

| accepted

Problem Calculating Nonlinear Indefinite Integral
Since you neglected to tell us values for a and b, I'll just set them to a=2 and b=1. You have also two definitions for g(x), so...

8 jours ago | 0

| accepted

Plot electric potential distribution
Well I kind of get what you want, but you do seem to have a singularity at the origin when r=0. p = 1; sig = 1; % constants th...

environ un mois ago | 1

fitting scatter data into multiple cosine functions
Yes, but it is always tricky to fit such sinusoids without good starting estimates. I'm going to use lsqcurvefit from the optimi...

environ 2 mois ago | 0

Question

Is the Moving Variance block in Simulink vectorised?
I'm using the Moving Variance block in Simulink and it seems to be giving unexplained results for vectorised (muxed) inputs. My ...

2 mois ago | 0 answers | 0

0

2D Convolution - Sobel Filter. What is wrong?
Here's my (old) code for a sobel filter: img = imread('Davis_Hall.jpg'); % 'Davis_Hall.jpg':color image X =double(rgb2gray(i...

3 mois ago | 0

Related to Matlab Step Comment
You are possibly confusing the closed loop response with the open loop. The closed loop gain is 1/4, so the error is 1-1/4 = 3/4...

3 mois ago | 1

| accepted

how to solve matrix differential equation
Here's my attempt: M = [ 97227/21875, 0, 82269/109375, 0, 67311/43750, 0, -97227/218750, 0; ... 0, 97227/21875, 0, 822...

4 mois ago | 1

From Fir2 output to transfer function
tr_fun should be a discrete transfer function, (not continuous) say tr_fun=tf(BF,1,1); % Generate a DISCRETE Transfer Function...

4 mois ago | 0

How to identify corners of data plotted to estimate data between corners
How about using interp1 with the 'pchip' option? You did speculate the interpolation is overly complex, but then again it is onl...

4 mois ago | 1

| accepted

How I can find co ordinates of intersection of two curves ?
Oh, looks like this is a common question: See (from 2016)https://au.mathworks.com/matlabcentral/answers/318475-how-to-find-the-i...

4 mois ago | 1

Iterative Method of solving multiple variables without using the curvefit tool
Here's my attempt at fitting uusing lsqcurvefit from the optimisation toolbox. First I generate some synthetic data to subsequ...

5 mois ago | 0

how can i calculate the angle of each gradient and x axis?
You could always differentiate to get the slope and then take the arc-tan. syms x real s = 1/(1 + exp(-2*(x-4))) dsdx = dif...

5 mois ago | 0

Extracting the daily time series from a text file to a single column
A quick hack is: fname = 'D59.txt' Flow = []; fid=fopen(fname); Yr = 2012; R = []; while 1 tline = fgetl(fi...

5 mois ago | 0

| accepted

How to Rearrange data?
Try something like the following: First read in the data file (which I've called numData.txt). I've used your test data above. ...

5 mois ago | 0

| accepted

(Food for Thought) Peak Detection on ECG signal
OK, here's my (ugly) solution (withou using findpeaks). I've re-named your data set, and did some pre-processig. load ecg y =...

6 mois ago | 0

| accepted

Simulating a FIR filter
Use the filter command. In this case your denominator, (A), coefficients are just 1. x = [ 0 0 0 1 1 1 1 1 ]' % input b = [...

6 mois ago | 1

| accepted

Obtain input knowing output and transfer function
Since you have a simple transfer function, here's what I would do. (1) First generate some data that we will subsequently try ...

6 mois ago | 0

curve fitting of non linear data for strain-time of creep experiment
Here's my attempt: load data1 t = time; e = strain; plot(t,e) %% Now try curve fit creep = @(p, t) p(1)*(1-exp(-p...

6 mois ago | 0

| accepted

Trouble with time delay differential equation
I've only made one small change to your function. You have tried to index a non-integer number back, but the "Z" variable holds ...

6 mois ago | 0

How can I simultaneously solve for these equations?
If you have the symbolic toolbox, a straight forward approach works, although it takes some time ... % ODE propblem syms B...

7 mois ago | 0

Error when using ode15
It is good practice to try and avoid globals, but rather pass them as auxilary variables into the function. Fao = 1487.351; % ...

7 mois ago | 0

| accepted

Plotting Step Response of a Differential Equation
You are close to a viable solution. Solve for an analytical solution, then make it a Matlab (anonymous) function with matlabFunc...

7 mois ago | 0

Runge-Kutta 4th order method
First up, you will need a much smaller step size to get an accurate solution using this explicit RK4 (with no error control). I ...

7 mois ago | 0

| accepted

How to draw Fig. 1 from the attached pdf with this code
I didn't bother draw the other 3 lines, but you just ned to make the necessary changes to gamma for that. If you run somethi...

7 mois ago | 0

| accepted

fsurf not drawing cant explain
Do you really need to use the symbolic toolbox? Why not just plot a numerical surface? But you are right, your numbers are widel...

7 mois ago | 0

Fitting Monod Equation with ODE45 to data using lsqcurvefit function
Here's my solution. First set up the measured data, plot it, and look to see that it is all OK. t = [0 3 5 8 9.5 11.5 14 16 ...

8 mois ago | 0

| accepted

Solve integrals with Matlab
Do you have the symbolic toolbox? If so, start with the int command. You can do both indefinite and definite integration. sy...

8 mois ago | 1

| accepted

MATLAB Help - Rolling Dice Simulation / Central Limit Theorem
A quick hack using randi and cumsum is N = 100 x = randi(6,N,6); cs = cumsum(x')'; for i=1:6 subplot(3,2,i) histogram...

8 mois ago | 0

| accepted