Community Profile

# Torsten

Last seen: Today Active since 2013

#### Content Feed

View by

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...

15 hours ago | 0

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 ...

1 day ago | 0

| accepted

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)

1 day ago | 0

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)

1 day ago | 0

| accepted

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...

1 day ago | 0

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".

1 day ago | 0

ODE45 - must return a column vector.

2 days ago | 0

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...

2 days ago | 1

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

3 days ago | 0

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

3 days ago | 0

| accepted

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))...

3 days ago | 0

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

3 days ago | 0

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...

3 days ago | 0

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 = ...

3 days ago | 1

| accepted

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".

3 days ago | 2

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...

3 days ago | 0

| accepted

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)).*(...

3 days ago | 0

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...

3 days ago | 0

| accepted

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...

4 days ago | 3

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...

4 days ago | 0

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; ...

4 days ago | 0

| accepted

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...

5 days ago | 1

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

5 days ago | 0

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...

5 days ago | 1

| accepted

3D Contour plot. Filled slices/surfaces are black
The mesh lines are black and cover the colors. To remove them, use gr = slice(Xg, Yg, Zg, Vg,[],[],zslice); set(gr,'Edgecolor...

7 days ago | 0

Hoe can I convert my script to a function?
If the coefficients don't change in the course of the simulation, you should read "Coefficient.mat" only once and pass the matri...

7 days ago | 0

| accepted

a = b<0
"a" is a logical variable. It is set to "true" (=1) if b<0, else it is set to "false" (=0). The next operation c = 3*a works a...

7 days ago | 2

How to compare previous results in iteration?
You should never check whether two variables are equal. Always check whether they are approximately equal: difference = 1.0; o...

7 days ago | 1