
Torsten
Statistics
0 Domande
5.037 Risposte
RANK
19
of 275.924
REPUTAZIONE
12.254
CONTRIBUTI
0 Domande
5.037 Risposte
ACCETTAZIONE DELLE RISPOSTE
0.00%
VOTI RICEVUTI
1.365
RANK
of 125.707
CONTRIBUTI
0 Problemi
0 Soluzioni
PUNTEGGIO
0
NUMERO DI BADGE
0
CONTRIBUTI
0 Post
CONTRIBUTI
0 Pubblico Canali
VALUTAZIONE MEDIA
CONTRIBUTI
0 Punti principali
NUMERO MEDIO DI LIKE
Content Feed
Nested Function is not working correctly with ode45
You can also get an analytical solution for v since your piecewise equation is easily integrated. EGA215_1234567_A1_Q2 functio...
circa 13 ore fa | 0
Constrained Optimization Problem With Obj Function Which Is Numerical Solution To A System Of ODE's
You never use the 53 solution variables of u and v in your code. And you solve the following ODEs for u and v du/dt = u, u(0)...
circa 16 ore fa | 0
Each element of INTCON must be an integer in the range [1 numVars], ([1 1]). What is the meaning??
f = c' * P_max + c0; % objective function is a scalar value. Thus intlinprog assumes you want to solve a problem with only one ...
circa 20 ore fa | 1
Error using symfun/subsref (line 203) Invalid argument at position 2.
N(i) = symfun(((1/(2*A))*(a(i)+b(i)*x+c(i)*y)),[x y]); instead of N = symfun(((1/(2*A))*(a(i)+b(i)*x+c(i)*y)),[x y]); And N...
circa 20 ore fa | 0
I have these 3 columns and I want to divide p/z and put the p/z values in a new columns
p=[2080,1885,1620,1205,888,645]; z=[0.759,0.767,0.787,0.828,0.866,0.900]; gp=[0,6.873,14.002,23.687,31.009,36.207]; new_row =...
1 giorno fa | 0
How to obtained matris vector solution using fminunc function?
You will have to use fminunc three times: c_1=[3;5;8]; c_2=[2;6;8]; x1_sol = zeros(size(c_1)); x2_sol = zeros(size(c_1)); f...
2 giorni fa | 0
| accettato
Applying updated initial conditions at every points of tspan and solve system of ODE using ode45
Then why do you solve for v1 if you want to prescribe it ? Simply set v1 = u0x + vp*p1 in your code and determine displacement...
3 giorni fa | 0
Using lsqnonlin for deconvolution
lsqnonlin tries to adjust x, and you reset it to 0: function F=insulin(x) ... x=zeros(89,1); This destroys the complete fitt...
3 giorni fa | 1
The following error was reported evaluating the function in FunctionLine update: Unable to convert expression into double array.
Your calculations partially produce complex and NaN values. And you cannot prevent this by trying to force them to be real. pv1...
3 giorni fa | 0
Solve of implicity equations
You should first check whether a zero of your function really exists. It doesn't seem to be the case: alpha = 0.14; A = 3300; n...
3 giorni fa | 0
why do I get this error in the simulation
Before calling "interp1", use % Sort array oxygen_energy_tot in ascending order [oxygen_energy_tot_sorted,idx_sort] = sort(oxy...
4 giorni fa | 0
| accettato
ode45 with a array of vector
a1 = ...; a2 = ...; a3 = ...; b1 = ...; b2 = ...; b3 = ...; c1 = ...; c2 = ...; c3 = ...; vp = ...; tValues = linspace...
4 giorni fa | 0
function handle array Problem
result = integral(@fun,0,1000,'ArrayValued',1) function value = fun(a) m = [1,2]; U(3) = a + 2; V(3) = a - 2; for...
4 giorni fa | 0
| accettato
i want to solve a set of homogeneous linear equation
Or use the following code to produce an optimal solution in the least-squares sense: A =1.0e+04*[6.6064,-3.5642,0,0;-3.5642,6.6...
4 giorni fa | 0
is there any way to get a 2nd order polynominal through 2 points(each end point)
rng("default") syms a b c x y n = 30; % number of data points xdata = sym('xdata',[1 n]); % x values ydata = sym('ydat...
4 giorni fa | 0
Error using plot Vectors must be the same length
As you can see above, S and f have a different number of elements. So they cannot be plotted against each other.
5 giorni fa | 1
incorrect dimensions for matrix multiplication in piecewise function in matlab
tau = 1.0; V_D = 0.72; vi = @(t) sin(2*pi*t); f = @(t) max(vi(t)-V_D,0).*exp(-t/tau); t = 0:0.01:2*pi; plot(t,f(t))
5 giorni fa | 0
How do I store an ODE in a vector?
Either you use the code below and work with the time derivatives after the integration to decide when the curves have flattened ...
6 giorni fa | 0
How to implement the Gaussian radial basis function in MATLAB ?
Before you start the optimization, check whether K in your case is positive-definite. Lt = 10; p = rand(Lt,4); sigma = 1.0; ...
6 giorni fa | 0
| accettato
Interface boundary condition in transient heat problem
Here is a sample code for what I wrote above: % Set parameters x1start = 0.0; x1end = 0.25; x2start = x1end; x2en...
7 giorni fa | 0
| accettato
How to solve below PDE based model using finite difference scheme.
Take a grid (ai,xj) over your rectangular region with 0 = a1 < a2 < ... < an = 60 0 = x1 < x2 < ... < xm = 40 Approximate the...
7 giorni fa | 2
Error: Array indices must be positive integers or logical values. implementing fletcher reeves method
[xmin,valormin]=fletcherreeves2([5;5],@(x) x(1)^2-x(1)*x(2)+4*x(2)^2+3*x(1)+2,@(x) [2*x(1)-x(2)+3;-x(1)+8*x(2)],10^-(6)) func...
8 giorni fa | 0
Why do I get 'Array indices must be positive or logical values' in symsum function?
h=0.05; U=0.02; delta_t=1; v=1.846*10^-5; d=1.117; t=50; y=0.2; P=pi; n = 1:300; Fn = 1./n.*exp(-(n.^2*P^2*v*t)/(d*h^2)...
8 giorni fa | 0
| accettato
How to get numerical values of nonlinear implicit function?
k = -1.5; F = @(t,y) log(abs(y)) - y.^2/2 - t - k; fimplicit(F,[0,4],'*') t = 0:0.01:0.99; y0 = 2.0; y = zeros(size(t)); f...
8 giorni fa | 0
| accettato
How to write the equation in matlab I= A A* T² exp (qob/kT) exp (q(V - Rs I /kT)) (exp (-1/n) ? Why is the curve not shown in the drawing?
Maybe you mean I = A*a*T^2*exp(q*phi_b/(K*T))*exp(q*(V-I1*Rs/(K*T)))*exp(-1/n) ?
8 giorni fa | 0
How to solve three coupled partial differential equations using matlab code?
Use MATLAB's "pdepe".
8 giorni fa | 0
Solving economic dispatch problem
The message from "quadprog" says that no feasible solution can be found. So you will have to reconsider your constraints - it se...
8 giorni fa | 0
| accettato
it says it cant plot because my vectors aren't the same size but they are
Check whether size(yvar) equals size(xvar). I think this cannot be true because the following code works: yvec=[6.67 17.33 42...
8 giorni fa | 0
I'm trying to run a code to plot conversion vs catalyst mass in a packed bed reactor, but something is wrong.
You tell the integrator that you have four additional inputs to your function [W,F]=ode45(@(t,F)ODEfun(t,F,T,P_total,F0_total,R...
9 giorni fa | 0
Bisection function to determine the value of (k) not working.
Please fill in the missing values and try again. If you don't succeed, it usually helps to plot f in the range 1e-12:1-1e-12 fo...
9 giorni fa | 0
| accettato