Community Profile

photo

Torsten


Last seen: Today Attivo dal 2013

Statistics

  • Most Accepted 2022
  • Explorer
  • Master
  • 36 Month Streak
  • Revival Level 2
  • Knowledgeable Level 5
  • First Answer

Visualizza badge

Content Feed

Visto da

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
ode45 with a array of vector
a1 = ...; a2 = ...; a3 = ...; b1 = ...; b2 = ...; b3 = ...; c1 = ...; c2 = ...; c3 = ...; vp = ...; tValues = linspace...

4 giorni fa | 0

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Risposto
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

Carica altro