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))

13 Tage 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...

14 Tage 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...

etwa ein Monat 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...

2 Monate 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 Monate 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 Monate 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...

4 Monate 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 Monate 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 Monate 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...

5 Monate 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...

5 Monate 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 Monate 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 Monate 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 Monate 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. ...

6 Monate 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 Monate 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 Monate 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 Monate 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 Monate 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 Monate 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 Monate 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 Monate 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 Monate 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 Monate 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...

8 Monate 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...

8 Monate 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 Monate 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 Monate 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 Monate ago | 0

| accepted