Answered

problem in creating function file using ode45

Are you looking for something like this? (I might have misinterpreted your equations!) O=1; a=pi/2; g=9.81; L=100; tspan ...

problem in creating function file using ode45

Are you looking for something like this? (I might have misinterpreted your equations!) O=1; a=pi/2; g=9.81; L=100; tspan ...

1 year ago | 0

Answered

how can eliminate negative data in erf function?

Should term2(i) = (sqrt(lemdaD)./(2.*e1(i))).*exp(-(v1(i).*sqrt_avgsnr - v2(i).* mu).^2./(2.*e1(i).^2)).*.... erf((e4(...

how can eliminate negative data in erf function?

Should term2(i) = (sqrt(lemdaD)./(2.*e1(i))).*exp(-(v1(i).*sqrt_avgsnr - v2(i).* mu).^2./(2.*e1(i).^2)).*.... erf((e4(...

1 year ago | 0

| accepted

Answered

Why is my plot not showing for improved Euler method?

You need to keep track of each step. Your original simply overwrites the values of t and y at each step. dy=@(t,y) -(0.1/10^(1...

Why is my plot not showing for improved Euler method?

You need to keep track of each step. Your original simply overwrites the values of t and y at each step. dy=@(t,y) -(0.1/10^(1...

1 year ago | 1

Answered

Why I couldn't plot this graph f=(@(x,y) (x.^2)+(x*y)+(y.^2)<=68200.^2/3);

This will give you the outline (I've assumed you want 68200^(2/3) rather than 68200^2/3) f=@(x,y) x.^2+x.*y+y.^2-68200.^(2/3); ...

Why I couldn't plot this graph f=(@(x,y) (x.^2)+(x*y)+(y.^2)<=68200.^2/3);

This will give you the outline (I've assumed you want 68200^(2/3) rather than 68200^2/3) f=@(x,y) x.^2+x.*y+y.^2-68200.^(2/3); ...

1 year ago | 0

Answered

I want to solve a simple "y=mx+c" with initial condition to find value of "m" and "c" but i am getting error by using following code. Any idea to tackle this kind of problem?

These are simple linear equations in m and c so set: X = [x(1) 1; x(2) 1]; Y = [y(1); y(2)]; and solve using X\Y to get m a...

I want to solve a simple "y=mx+c" with initial condition to find value of "m" and "c" but i am getting error by using following code. Any idea to tackle this kind of problem?

These are simple linear equations in m and c so set: X = [x(1) 1; x(2) 1]; Y = [y(1); y(2)]; and solve using X\Y to get m a...

1 year ago | 0

Answered

I have been trying to code a solar pv module for 1 diode model, the following is my code, can someone tell me the mistake ive made and the correction to it pls.

Like this? Note: I've made kref positive to avoid the I's being complex. This might be completely incorrect!! Pmpp = 50; ...

I have been trying to code a solar pv module for 1 diode model, the following is my code, can someone tell me the mistake ive made and the correction to it pls.

Like this? Note: I've made kref positive to avoid the I's being complex. This might be completely incorrect!! Pmpp = 50; ...

1 year ago | 1

| accepted

Answered

Reorganization of data in matrix

Like this? M = [1 1 0.270588872500870 1 1 0.191034400222846 1 1 0.316821350897659 1 1 0.391794699189989 1 2...

Reorganization of data in matrix

Like this? M = [1 1 0.270588872500870 1 1 0.191034400222846 1 1 0.316821350897659 1 1 0.391794699189989 1 2...

1 year ago | 0

Answered

LogLog Plot is Blank

Like this? (I've assumed you want log base 10 everywhere): %Paschen Curve Plotting for H2 p = 1.2; %Torr d = 20; %cm GamSE ...

LogLog Plot is Blank

Like this? (I've assumed you want log base 10 everywhere): %Paschen Curve Plotting for H2 p = 1.2; %Torr d = 20; %cm GamSE ...

1 year ago | 1

Answered

How to combine multiple numbers into 1 number?

Do you mean R1 = 2; R2 = 4; Zr1 = [0 0 0]; R = R1*10^4+R2*10^3+Zr1(1)*10^2+Zr1(2)*10+Zr1(3) Alternatively: R1 = 2; R2 = 4; Z...

How to combine multiple numbers into 1 number?

Do you mean R1 = 2; R2 = 4; Zr1 = [0 0 0]; R = R1*10^4+R2*10^3+Zr1(1)*10^2+Zr1(2)*10+Zr1(3) Alternatively: R1 = 2; R2 = 4; Z...

1 year ago | 0

Answered

Runge Kutta Fourth order and Interpolation

Shouldn't f1 = @(z,r,ed,p11,m11) -(ed+p11)*z*m11.^3 *((r+1)./(((ed+p11)*(r.^2 +z.^2).^6) -(r/p11))) ./(2*pi*(r.^2 + z.^2)) ; ...

Runge Kutta Fourth order and Interpolation

Shouldn't f1 = @(z,r,ed,p11,m11) -(ed+p11)*z*m11.^3 *((r+1)./(((ed+p11)*(r.^2 +z.^2).^6) -(r/p11))) ./(2*pi*(r.^2 + z.^2)) ; ...

1 year ago | 0

Answered

Newton Raphson method for a system of non-linear equations?

Here's a simple example for you to follow: % Functions f = @(x,y) x^2 + y - 5.94; % i.e. the equation is x^2+y=5.94 g = @(x,...

Newton Raphson method for a system of non-linear equations?

Here's a simple example for you to follow: % Functions f = @(x,y) x^2 + y - 5.94; % i.e. the equation is x^2+y=5.94 g = @(x,...

1 year ago | 0

| accepted

Answered

No feasible solution in optimisation in linprog

I don't know if this will fix your problem, but I think the first number in the second row of A1 should be 8 not 2.

No feasible solution in optimisation in linprog

I don't know if this will fix your problem, but I think the first number in the second row of A1 should be 8 not 2.

1 year ago | 0

Answered

Simple relaxation method Matlab

For the simple relaxation technique to work the A matrix needs to be diagonally dominant. Your A is not diagonally dominant - s...

Simple relaxation method Matlab

For the simple relaxation technique to work the A matrix needs to be diagonally dominant. Your A is not diagonally dominant - s...

1 year ago | 0

| accepted

Answered

Solve a system of 3 second order equations

Write your equations as 6 first order equations like so: where I've assumed you have pre-defined omegaout and thetaout as fun...

Solve a system of 3 second order equations

Write your equations as 6 first order equations like so: where I've assumed you have pre-defined omegaout and thetaout as fun...

1 year ago | 0

Answered

Unable to solve nonlinear equation using fsolve as the message shows No solution found

fminsearch makes a reasonable attempt: F1 = @(a,b) -9.135789053E+00-log(a)+166.2509975*a+5.84074229*b+... -166.2509975*a^2...

Unable to solve nonlinear equation using fsolve as the message shows No solution found

fminsearch makes a reasonable attempt: F1 = @(a,b) -9.135789053E+00-log(a)+166.2509975*a+5.84074229*b+... -166.2509975*a^2...

1 year ago | 2

Answered

How to solve and plot this set of differential equations?

If rho starts at 0, then drho/dtheta will be zero and rho will not change from zero according to your first ode. If rho and drh...

How to solve and plot this set of differential equations?

If rho starts at 0, then drho/dtheta will be zero and rho will not change from zero according to your first ode. If rho and drh...

1 year ago | 0

Answered

How can I add two differential equations to the system in a given time interval?

Like this? (though I'm not sure I've interpreted which constants apply during which interval correctly!) Tspan = 0:30; C0 =...

How can I add two differential equations to the system in a given time interval?

Like this? (though I'm not sure I've interpreted which constants apply during which interval correctly!) Tspan = 0:30; C0 =...

1 year ago | 1

| accepted

Answered

How do I plot Rosette-like diagram in this case?

Something like this? data=[10 10 30 100; 20 50 70 150; 30 45 60 120]; theta0 = deg2rad(data(:,2)); thetaf = deg2r...

How do I plot Rosette-like diagram in this case?

Something like this? data=[10 10 30 100; 20 50 70 150; 30 45 60 120]; theta0 = deg2rad(data(:,2)); thetaf = deg2r...

1 year ago | 1

| accepted

Answered

Iteration of multiple nonlinear functions

Here's a quick and dirty way. I'll leave you to modify the following to record all the temperatures at every step. %Given valu...

Iteration of multiple nonlinear functions

Here's a quick and dirty way. I'll leave you to modify the following to record all the temperatures at every step. %Given valu...

1 year ago | 0

Answered

using num2str for subplot titles

Try title(['Q1S', num2str(k)]) i.e. enclose the two terms within square brackets.

using num2str for subplot titles

Try title(['Q1S', num2str(k)]) i.e. enclose the two terms within square brackets.

1 year ago | 0

| accepted

Answered

control and simulation of petroleum distillation column

The following works, but probably doesn't contain all the right data or equations! % CALCULATE MATERIAL BALANCES FOR THE DISTIL...

control and simulation of petroleum distillation column

The following works, but probably doesn't contain all the right data or equations! % CALCULATE MATERIAL BALANCES FOR THE DISTIL...

1 year ago | 0

| accepted

Answered

Error with plot command

You haven't specified function atten_func_dist, so we can't run your code to test it. You specify fre as a constant, then us...

Error with plot command

You haven't specified function atten_func_dist, so we can't run your code to test it. You specify fre as a constant, then us...

1 year ago | 2

| accepted

Answered

How to combine multiple curve fits in one plot

Use hold on after plotting the first figure

How to combine multiple curve fits in one plot

Use hold on after plotting the first figure

1 year ago | 0

Answered

How to represent [airfoil] coordinates as a polynomial

How about using a spline fit? After loading your data and calling the first column x and the second column y try xfit = linspa...

How to represent [airfoil] coordinates as a polynomial

How about using a spline fit? After loading your data and calling the first column x and the second column y try xfit = linspa...

1 year ago | 0

Answered

1D melting problem with Neumann's analytical solution

eps is a built-in constant so better not to use it as a variable. Your png files show it is xi not epsilon. Try the following ...

1D melting problem with Neumann's analytical solution

eps is a built-in constant so better not to use it as a variable. Your png files show it is xi not epsilon. Try the following ...

1 year ago | 0

| accepted

Answered

Overlay plots of different angles of incidence.

Remove the word figure and add hold on after the plot grid on [Fr,In] = meshgrid(Frequency,Incident); plot3(Fr,In,Reflectivit...

Overlay plots of different angles of incidence.

Remove the word figure and add hold on after the plot grid on [Fr,In] = meshgrid(Frequency,Incident); plot3(Fr,In,Reflectivit...

1 year ago | 0

Answered

How can write in MATLAB the following expression?

Try ket0 = [1; 0; 0] ket1 = [0; 1; 0]

How can write in MATLAB the following expression?

Try ket0 = [1; 0; 0] ket1 = [0; 1; 0]

1 year ago | 0

| accepted

Answered

I don't understand why does it give error?

Put the functions at the end. Tidy up one or two minor errors to get: %########### BF-PSO.m file ############# %####### Butter...

I don't understand why does it give error?

Put the functions at the end. Tidy up one or two minor errors to get: %########### BF-PSO.m file ############# %####### Butter...

1 year ago | 0

Answered

How to count steps represented by the acceleration data given in text file?

Like this (you can plot the other two - look up help on subplot): M = csvread('walking.txt'); x = M(:,1); y = M(:,2); z = M(:,...

How to count steps represented by the acceleration data given in text file?

Like this (you can plot the other two - look up help on subplot): M = csvread('walking.txt'); x = M(:,1); y = M(:,2); z = M(:,...

1 year ago | 0

| accepted