Answered

ODE15s with non-constant Jacobian

JacFun = J(u_int, deltaZ, Nz, cFeed); options = odeset( 'Jacobian' , @(t,x)JacFun(t, x, u_int, deltaZ, Nz, cFeed), 'RelTol',1e-...

ODE15s with non-constant Jacobian

JacFun = J(u_int, deltaZ, Nz, cFeed); options = odeset( 'Jacobian' , @(t,x)JacFun(t, x, u_int, deltaZ, Nz, cFeed), 'RelTol',1e-...

9 giorni ago | 1

| accepted

Answered

Integration of the expected value

b = 1; N = 2000000; for i = 1:N p(i) = b*rand; pe(i) = p(i) + (b-p(i))*rand; p2(i) = p(i) + pe(i); end ecdf(p...

Integration of the expected value

b = 1; N = 2000000; for i = 1:N p(i) = b*rand; pe(i) = p(i) + (b-p(i))*rand; p2(i) = p(i) + pe(i); end ecdf(p...

9 giorni ago | 1

| accepted

Answered

Why do I get the wrong value when I increase N?

a = 0; % Lower Boundary b = 10; % Upper Boundary fun = @(x)sin(pi*x/10); % Function zeta = [-sqrt(3/5) 0 sqrt(3/5)]; % absci...

Why do I get the wrong value when I increase N?

a = 0; % Lower Boundary b = 10; % Upper Boundary fun = @(x)sin(pi*x/10); % Function zeta = [-sqrt(3/5) 0 sqrt(3/5)]; % absci...

9 giorni ago | 0

Answered

Trouble with least squares fit to 3D ODE model with 2D observed data

[~,Sv] = ode45(@(t,y) growth(t,y),t,y0); instead of Sv = ode45(@(t,y) growth(t,y),t,y0);

Trouble with least squares fit to 3D ODE model with 2D observed data

[~,Sv] = ode45(@(t,y) growth(t,y),t,y0); instead of Sv = ode45(@(t,y) growth(t,y),t,y0);

9 giorni ago | 0

Answered

How can I change the definition of x to find the value of y(3.5112)( i mean y(x1)) for this question ?

I don't know if it helps: %%main func clear all; close all; clf; clc; % Define the function x = linspace(-10,10,100); ...

How can I change the definition of x to find the value of y(3.5112)( i mean y(x1)) for this question ?

I don't know if it helps: %%main func clear all; close all; clf; clc; % Define the function x = linspace(-10,10,100); ...

9 giorni ago | 0

Answered

Trouble defining a set of parameters while using ode45 to solve a system of differential equations

This page should help you: https://de.mathworks.com/help/matlab/math/ode-event-location.html

Trouble defining a set of parameters while using ode45 to solve a system of differential equations

This page should help you: https://de.mathworks.com/help/matlab/math/ode-event-location.html

10 giorni ago | 0

Answered

Finding the turning point on a curve

v = [-0.000499915992890149 -0.00136716055731299 -0.00312202879147695 -0.00415143850349407 -0.003723995799058...

Finding the turning point on a curve

v = [-0.000499915992890149 -0.00136716055731299 -0.00312202879147695 -0.00415143850349407 -0.003723995799058...

10 giorni ago | 0

Answered

how to substitute a symbolic differential equation?

syms x(t) y(t) lambda h = 2*diff(x,t,2) + (diff(y,t))^2; xtt = diff(x,t,2); yt = diff(y,t); h_subs = subs(h,[xtt,yt],[lambda...

how to substitute a symbolic differential equation?

syms x(t) y(t) lambda h = 2*diff(x,t,2) + (diff(y,t))^2; xtt = diff(x,t,2); yt = diff(y,t); h_subs = subs(h,[xtt,yt],[lambda...

10 giorni ago | 0

| accepted

Answered

How to pre-allocate a double for following script ?

What about syncat = zeros(100000,2) ?

How to pre-allocate a double for following script ?

What about syncat = zeros(100000,2) ?

10 giorni ago | 1

| accepted

Answered

Getting a system of equations and the using ODE45

fun = @(x) Xi.'*[1;x(1);x(2);x(3);x(1)^2;x(1)*x(2);x(1)*x(3);x(2)^2;x(2)*x(3);x(3)^2;x(1)^3;x(1)^2*x(2);x(1)^2*x(3);x(1)*x(2)^2;...

Getting a system of equations and the using ODE45

fun = @(x) Xi.'*[1;x(1);x(2);x(3);x(1)^2;x(1)*x(2);x(1)*x(3);x(2)^2;x(2)*x(3);x(3)^2;x(1)^3;x(1)^2*x(2);x(1)^2*x(3);x(1)*x(2)^2;...

11 giorni ago | 1

| accepted

Answered

i need to find the optimize values of R for a values in an interval but R values are in form of the root of the equation; how to solve this sir? for simple equation it was ok.

I think this is the safest way: syms R a p = -282091600255531819/9878934086330590709531986654911361482361187039433493704830155...

i need to find the optimize values of R for a values in an interval but R values are in form of the root of the equation; how to solve this sir? for simple equation it was ok.

I think this is the safest way: syms R a p = -282091600255531819/9878934086330590709531986654911361482361187039433493704830155...

11 giorni ago | 0

| accepted

Answered

Index in position 3 exceeds array bounds (must not exceed 1)

"orph" has dimension 1 x 20000 because you overwrite the image by the command orph = [i,j]. Then you try to access orph(:,:,2)...

Index in position 3 exceeds array bounds (must not exceed 1)

"orph" has dimension 1 x 20000 because you overwrite the image by the command orph = [i,j]. Then you try to access orph(:,:,2)...

11 giorni ago | 1

| accepted

Answered

i need to find the optimize values of R for a values in an interval but R values are in form of the root of the equation; how to solve this sir? for simple equation it was ok.

Try "fmincon" for the problem min: R s.c. f(R,a) = 0 (f is the complete expression between "t=solve(" and "==0,R)) a_left <...

i need to find the optimize values of R for a values in an interval but R values are in form of the root of the equation; how to solve this sir? for simple equation it was ok.

Try "fmincon" for the problem min: R s.c. f(R,a) = 0 (f is the complete expression between "t=solve(" and "==0,R)) a_left <...

11 giorni ago | 0

Answered

Error using cat Dimensions of arrays being concatenated are not consistent.

Xtrain = cat(2,Xtrain,smat.'); would work, I guess.

Error using cat Dimensions of arrays being concatenated are not consistent.

Xtrain = cat(2,Xtrain,smat.'); would work, I guess.

12 giorni ago | 0

Answered

How to find L from the given code?

lags = 0.6; k1 = 3; k_1 = 1; k2 = 2.5; k3 = 1; k_3 = 1; k4 = 2; k5 = 1; E1 = 1; E2 = 2; K1 = (k_1+k2)/k1; K2 = (k_3+k...

How to find L from the given code?

lags = 0.6; k1 = 3; k_1 = 1; k2 = 2.5; k3 = 1; k_3 = 1; k4 = 2; k5 = 1; E1 = 1; E2 = 2; K1 = (k_1+k2)/k1; K2 = (k_3+k...

12 giorni ago | 1

| accepted

Answered

i have to plot the percentage of heads by the times token is trowed times, in a game of head or tail, but a need help for doing that. What can I use in my plot?

N = 500; x=rand(1,N); for i = 1:length(x) if x(i) < 0.5 x(i) = 1; % cara else x(i) = 0; % coroa ...

i have to plot the percentage of heads by the times token is trowed times, in a game of head or tail, but a need help for doing that. What can I use in my plot?

N = 500; x=rand(1,N); for i = 1:length(x) if x(i) < 0.5 x(i) = 1; % cara else x(i) = 0; % coroa ...

12 giorni ago | 0

| accepted

Answered

ode15s solver produces wrong solution

t0 = 0; % initial time tf = 1; % final time nt = 20; % number of output times xl = 0; % left point xr = 1; % right point nx...

ode15s solver produces wrong solution

t0 = 0; % initial time tf = 1; % final time nt = 20; % number of output times xl = 0; % left point xr = 1; % right point nx...

12 giorni ago | 0

| accepted

Answered

How to store n iterations of randomly generated matrices

n = 10 ; bounds = [1 100] ; %interval over which the random number generator will produce values randA = zeros(3,3); randA_al...

How to store n iterations of randomly generated matrices

n = 10 ; bounds = [1 100] ; %interval over which the random number generator will produce values randA = zeros(3,3); randA_al...

12 giorni ago | 0

Answered

Numerical answer to Partial derivative of Hopfield-style Energy in symbolic expression

p is (4x1), not (16x1). And keep r_mat as matrix, then E = -0.5*p.'*r_mat*p - c.'*p ( r_mat can be assumed symmetric ). By...

Numerical answer to Partial derivative of Hopfield-style Energy in symbolic expression

p is (4x1), not (16x1). And keep r_mat as matrix, then E = -0.5*p.'*r_mat*p - c.'*p ( r_mat can be assumed symmetric ). By...

14 giorni ago | 0

| accepted

Answered

Error from ode15s

psi = 0.63; psip = 0.27; L = 150e-6; Ts = 60; k = 0.00158; q0 = 4.62e3; C0 =1.03; gama = 5021.8; Jw = 8.5e-6; a = psi*...

Error from ode15s

psi = 0.63; psip = 0.27; L = 150e-6; Ts = 60; k = 0.00158; q0 = 4.62e3; C0 =1.03; gama = 5021.8; Jw = 8.5e-6; a = psi*...

14 giorni ago | 0

Answered

Array indices must be positive integers or logical values.

e = exp(-SNR_index.^2./2) SNR_index is undefined. compute_the_Q = integral(@(SNR_index) (1/sqrt(2*pi)).*e, SNR_index,Inf,'Arra...

Array indices must be positive integers or logical values.

e = exp(-SNR_index.^2./2) SNR_index is undefined. compute_the_Q = integral(@(SNR_index) (1/sqrt(2*pi)).*e, SNR_index,Inf,'Arra...

15 giorni ago | 0

Answered

How can I simulate Buffon's needle on MATLAB?

https://simonensemble.github.io/2018-04/buffon

How can I simulate Buffon's needle on MATLAB?

https://simonensemble.github.io/2018-04/buffon

15 giorni ago | 0

| accepted

Answered

Solve Three non-linear equations for 3 variables system

You might want to try a numerical solver, but the two I tested also didn't succeed. fun = @(A,B,C)[72000-(B*4000)/(log(A*4000)-...

Solve Three non-linear equations for 3 variables system

You might want to try a numerical solver, but the two I tested also didn't succeed. fun = @(A,B,C)[72000-(B*4000)/(log(A*4000)-...

15 giorni ago | 0

Answered

Coupled ODE using ODE15s with multiple outputs

When communicating with ODE15S, the solution must be treated as one vector of size 200, not as two vectors of size 100. E.g. f...

Coupled ODE using ODE15s with multiple outputs

When communicating with ODE15S, the solution must be treated as one vector of size 200, not as two vectors of size 100. E.g. f...

16 giorni ago | 0

| accepted

Answered

Is it possible to (non-linear) minimize x^y by choosing both x and y?

fun = @(x) x(1)^sqrt(x(2)); lb = [1.01,1.01]; ub = [100,2]; x0 = [1.01,1.01]; sol = fmincon(fun,x0,[],[],[],[],lb,ub,@nonlco...

Is it possible to (non-linear) minimize x^y by choosing both x and y?

fun = @(x) x(1)^sqrt(x(2)); lb = [1.01,1.01]; ub = [100,2]; x0 = [1.01,1.01]; sol = fmincon(fun,x0,[],[],[],[],lb,ub,@nonlco...

16 giorni ago | 0

Answered

What is wrong in script?

I would have progrmmed it this way, but the result seems to be the same. My guess is that there is a mathematical trick to extr...

What is wrong in script?

I would have progrmmed it this way, but the result seems to be the same. My guess is that there is a mathematical trick to extr...

17 giorni ago | 0

Answered

How to generate plot

a = 5; n = 30; P0 = 100; c = 3; s = cos((2*(1:n)-1)*pi/(2*n)); x = a*s; P = P0*exp(-0.5*(x/c).^2).*(1-x/c).^2; K1 = sqrt(...

How to generate plot

a = 5; n = 30; P0 = 100; c = 3; s = cos((2*(1:n)-1)*pi/(2*n)); x = a*s; P = P0*exp(-0.5*(x/c).^2).*(1-x/c).^2; K1 = sqrt(...

18 giorni ago | 1

Answered

What is an efficient way of finding the maximum value under a certain limit?

A = [1,2,3;4,5,6;7,8,9]; A(A>=8) = -Inf [m,i] = max(A,[],'All') [row,col] = ind2sub(size(A),i)

What is an efficient way of finding the maximum value under a certain limit?

A = [1,2,3;4,5,6;7,8,9]; A(A>=8) = -Inf [m,i] = max(A,[],'All') [row,col] = ind2sub(size(A),i)

18 giorni ago | 1

| accepted

Answered

How I can give condition & plot the solution of this differential equation. . . . . . . Please Guide

dsolve doesn't succeed. Thus use a numerical solver (bvp4c) to solve your equation. syms A theta(z) dtheta=diff(theta,z) d...

How I can give condition & plot the solution of this differential equation. . . . . . . Please Guide

dsolve doesn't succeed. Thus use a numerical solver (bvp4c) to solve your equation. syms A theta(z) dtheta=diff(theta,z) d...

20 giorni ago | 0

Answered

diffusion from point source boundary conditions

if i use the "s" variable for the solver, does this input it at r=0, or for every value of r? For every value of r (at least if...

diffusion from point source boundary conditions

if i use the "s" variable for the solver, does this input it at r=0, or for every value of r? For every value of r (at least if...

21 giorni ago | 1

| accepted