Community Profile

photo

Torsten


Last seen: Today Con actividad desde 2013

Statistics

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

Ver insignias

Content Feed

Ver por

Respondida
custom equation fit in curveFitter is not working
pc2 = [ 0 0.0001 0.0011 0.0216 0.2847 0.8708 1.0000]; lam_2 = [ 0 30 55 110 2...

31 minutos hace | 0

Respondida
numerical integration inf to inf
Set a=-1+eps; b=1-eps; instead of a=-1; b=1; For comparison: syms x e = 0.01; f = 1/sqrt((1-x^2)+(e/2)*(1-x^4)); valu...

alrededor de 2 horas hace | 0

Respondida
Bessel Equation VPASolve for the order
nu = 0:0.1:7; B = besselk(nu,3)-5; plot(nu,B) fun = @(nu)besselk(nu,3)-5; nu = fzero(fun,[6 7])

alrededor de 3 horas hace | 0

Respondida
How to plot non-quadratic functions?
I suggest you evaluate the three profit functions for a rectangular region [alpha_min,alpha_max] x [ F_min,F_max] and introduce ...

alrededor de 3 horas hace | 0

Respondida
Problem with undefined function
An alternative without the toolbox: logncdf_user = @(x,mu,sigma) 0.5*(erf((log(x)-mu)/(sqrt(2)*sigma))+1); x=0:0.1:10; hold o...

1 día hace | 0

Respondida
A system of nonlinear equations with three variables
The equations for x and y (eqnx and eqny) are polynomials of degree 3 in x resp. y. They might have several real solutions x an...

1 día hace | 0

| aceptada

Respondida
Way to solve AX=XB
dim = 4; X = sym('X',[dim dim]); A = rand(dim); B = A.'; [M, ~] = equationsToMatrix(A*X==X*B) if rank(M) < size(A,1)^2 N...

2 días hace | 0

Respondida
Trouble using lsqnonlin to fit the IV curve of a solar cell
You did not include global V Voc Jsc T k q J in the script part. clc clear global V Voc Jsc T k q J %reads input file ...

3 días hace | 0

| aceptada

Respondida
How do I represent this sym function as a written formula
syms s f = vpa(ilaplace(1/((s+7)*(s^3+4*s^2+6*s+5))),4)

3 días hace | 0

Respondida
Invalid use of operator. error, dont know what to change to fix it.
Elementwise operations are .* ./ and .^ .) does not exist. f = (x.^2).*exp(-x)

3 días hace | 0

| aceptada

Respondida
I want to check if something algebraically simplifies to 0 when i sub in values
Euler equations of gas dynamics ? Eigenvalues and Eigenvectors of the Jacobian are well-studied. Why reinvent the wheel ? syms...

3 días hace | 0

Respondida
How do I solve a linear optimization problem with a constrained output?
Define the X to be solved for as X = (u,x).Then you can constrain x according to your needs and continue using "linprog".

3 días hace | 0

Respondida
ODE45 - must return a column vector.
B = readtable('t13.xlsx'); t13 = table2array(B); c = readtable('ambient.xlsx'); ambient = table2array(c); d = readtable('sol...

3 días hace | 0

Respondida
How to solve optimization problems when the objective function includes a symbolic function
The usual way to deal with this problem is to define Func2Var1 and Func2Var2 as additional optimization variables and to specify...

4 días hace | 1

Respondida
Plot a graph from arrays
plot(hv, cell2mat(Xm))

4 días hace | 0

Respondida
Function with many outputs
X=ga(@cost,numel(X)); function J = cost(X); [Y_1,Y_2] = MyFun(X); J = Y_1.^2+Y_2.^2; end

4 días hace | 0

| aceptada

Respondida
Trial and error problem
The results are not convincing, but that's the way to go with three data points. fun1 = @(x)[0.0701-x(1)*exp(x(2)/(293.65-x(3))...

4 días hace | 0

Respondida
How to plot a graph correctly?
Try "fimplicit".

4 días hace | 0

Respondida
Why cant I get a stable Cranck Nicolson discretization?
We use a Dirichet boundary condition so that x=0 for t = 0 and a = 0. And which conditions at a = a_end ? Try "pdepe". If th...

4 días hace | 0

Respondida
Optimising Limits of Chain integrals
fun1 =@(g1) 208.4*exp(-0.02704.*g1)-293.9*exp(-0.1286.*g1); fun2 =@(g2) 108.3*exp(-0.01137.*g2)-622.8*exp(-0.1578.*g2); A = ...

4 días hace | 1

| aceptada

Respondida
Is it possible to add more constraints than the default ones in LSQLIN function?
Use "fmincon" with the constraint (C*x).*d >=0 and define this constraint in function "nonlcon".

4 días hace | 2

Respondida
Hi, i am trying plot my thesis graph using dde23 but i keep on getting this error, derivative and history function have different lengths.. HELP PLEASE!
z(:,1) is a 3x1 vector with z(1,1) = y1(t-0.02), z(2,1) = y2(t-0.02) and z(3,1) = y3(t-0.02), You have to decide which of the t...

4 días hace | 0

| aceptada

Respondida
How to solve error using integral (line 83) first input argument must be a function handle?
T = 1; m = 64; x = 1:0.1:2; y = 2:0.2:10; [X,Y] = meshgrid(x,y); g = @(u)sin(2*pi*u); phi = @(t,x,y) ((sign(x/y+g(t/T)).*(...

4 días hace | 0

Respondida
Double monod maximum rate estimation
% Initial guess for rmax rmax0 = ...; % Time is a (nx1) vector, % Sdata is a (nx2) matrix with measurements of c_CH4 (1st c...

5 días hace | 0

| aceptada

Respondida
How to calculate the symbolic integration of rational functions?
syms x f = x^4+x^3+x^2+x+1; fp = partfrac(1/f,x,'FactorMode','real') int(fp,x,0,1) vpaintegral(1/f,x,0,1) f = matlabFunctio...

5 días hace | 3

Respondida
I couldn't understand why this program took a long time to run, can anyone help me to solve the problem so that the code can run fast
The reason is this line: U = sort(rand(1,20).*(0.005)+rand(1,20).*(0.003),'ascend'); MATLAB integrators are not suited to cope...

5 días hace | 0

Respondida
Trying to use Newton's method until solution converges
% Variables Eps = 0.011; E = 73100; H = 662; n = 0.07; maxIter = 100; % # of iterations Tol = 1e-8; i = 1; error = 1.0; ...

5 días hace | 0

| aceptada

Respondida
Solve System of Equations to calculate the inverse kinematics
X = cos(q1)*(cos(q2 + q3) + cos(q2)); Y = sin(q1)*(cos(q2 + q3) + cos(q2)); Z = sin(q2 + q3) + sin(q2) + 1/5; Dividing equati...

6 días hace | 1

Respondida
I have a problem with solve command
syms x y [x y] = solve([3*x-y==2,x+y==1])

6 días hace | 0

Respondida
Write a row and column vector as matrix index
a = [4 5 3]; % Generate index vectors A and B [A,B] = ndgrid(a,a); A = A(:) B = B(:) % Use A and B as index vectors for a m...

6 días hace | 1

| aceptada

Cargar más