Community Profile

# David Wilson

### Auckland University of Technology

Last seen: 1 giorno fa Attivo dal 2011

Followers: 0   Following: 0

All

#### Feeds

Visto da

Risposto
Runge-Kutta 4th order method
Wow, you haven't given us too much to go on, so that makes a real challenge. First up, your 2nd order ODE is needlessly comple...

quasi 3 anni fa | 0

Risposto
Close-loop TF Calculation
You need a few multiplies. s = tf('s') Gc= 32.8802*(s+3)/s/(s+7)/(s+8.52) Now you can perform step or bode or whatever.

circa 3 anni fa | 1

Risposto
How to calculate the area between two curves
Try the following : I've used one of the many interpolation routines, but I notice that you have not sorted your (x2, y2) data, ...

oltre 3 anni fa | 0

| accettato

Risposto
Range-Kutta solving error.
Seems fine to me: Let's try to check the analhytical solution, and also find the derivative of F(t). That's pretty easy: clc...

oltre 3 anni fa | 1

Risposto
Why when I try to LLSQ fit to a quadratic the line of best fit has so many lines
A simple fix is to replace your independent variable with a variable ordered from min to max % xx = x xx = linspace(50, 230)';...

oltre 3 anni fa | 2

| accettato

Risposto
Filling the area between curves (loglog)
I too had some difficulty to do this elgantly. I changed the name of your data file, and commented out the material that didn't ...

oltre 3 anni fa | 0

Risposto
How to solve the riccati equations in matlab
Well, if Google translate is accurate, you seem to want to solve for the matrix Q. I'm assuming that everything is a matrix w...

oltre 3 anni fa | 1

| accettato

Risposto
How to solve the riccati equations in matlab
Do you have the control toolbox? I suspect that you could use are.m and related friends.

oltre 3 anni fa | 0

Risposto
How to find time constants from given 2+ order transfer function?
If you just want to display the transfer function in time-constant form, you can do the following: s = tf('s'); k = 45; % an...

oltre 3 anni fa | 1

Risposto
i am trying to made a code to int the numbers of time i want
Remove the z=x^2 from the loop. syms x n y z a; n= 5 z=x^2; for y=1:n z=int(z) end

quasi 4 anni fa | 1

| accettato

Risposto
How to obtain regression polynomial equation with more than 2 independent variables with degree 5.
OK, the fitting is easy, the plotting less so. Since you didn't give us too many hints, not did you give us any data to work wit...

circa 4 anni fa | 2

| accettato

Risposto
Text keeps going behind topographic map
Your mountains are covering the text (as they should). So you need to write your text above the highest peak, so that it is seen...

circa 4 anni fa | 0

| accettato

Risposto
function to produce line graph
I would strongly suggest you don't try a polynomial fit here. From the context (which you neglected to tell us), I suspect that ...

circa 4 anni fa | 0

| accettato

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

oltre 4 anni fa | 0

| accettato

Risposto
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...

oltre 4 anni fa | 0

| accettato

Risposto
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...

oltre 4 anni fa | 1

Risposto
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...

oltre 4 anni fa | 0

Domanda

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

oltre 4 anni fa | 1 risposta | 0

### 1

risposta

Risposto
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...

oltre 4 anni fa | 1

Risposto
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...

oltre 4 anni fa | 1

| accettato

Risposto
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...

oltre 4 anni fa | 1

Risposto
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...

oltre 4 anni fa | 0

Risposto
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...

oltre 4 anni fa | 1

| accettato

Risposto
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...

oltre 4 anni fa | 1

Risposto
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...

oltre 4 anni fa | 0

Risposto
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...

oltre 4 anni fa | 0

Risposto
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...

oltre 4 anni fa | 0

| accettato

Risposto
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. ...

quasi 5 anni fa | 0

| accettato

Risposto
(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 =...

quasi 5 anni fa | 0

| accettato