Community Profile

# Torsten

Last seen: Today Active since 2013

#### Content Feed

View by

T is not defined, dont know how to fix
After deleting the @ before ode23t, the code with @locO works for me. %ODEs n0=[0 0 0]; % Initial conditions options = odes...

ungefär 9 timmar ago | 0

How to optimise the computation of multiple integral for a specific integrand?
Precalculate q over the hypercube on a predefined regular mesh and use "trapz" dimensionwise. Look at the example "Multiple nu...

ungefär 15 timmar ago | 1

substitution in symbolic.
syms s w assume(w,'real') f = s + s^2; fs = subs(f,s,1i*w); fs_real = real(fs) fs_imag = imag(fs)

ungefär 16 timmar ago | 1

Fit chi distribution to measurement data
Square your measurement data, and you can fit them against a chi-squared distribution.

ungefär 16 timmar ago | 0

| accepted

Calculating and plotting conditional distribution.
You mean this lambda = 0.3; p_bar = 1; p_tilda = (1 - lambda)/(1 + lambda); S = 0:0.01:p_tilda; n = 1000000; icount = 0; ...

en dag ago | 1

| accepted

odeEvent doesn't terminate the integration process?
The functions you can prescribe as "value" should be continous. Using logicals leads to step functions that are not appropriate ...

en dag ago | 0

Giving a rectangular pulse in ode45
I can't evaluate the equations and initial conditions you use. This is what the solver returns - slight over-and undershooting. ...

en dag ago | 0

Giving a rectangular pulse in ode45
time=linspace(0,20,300); is wrong. You must concatenate the t-vectors returned from ode45 in the three calls in the same way ...

en dag ago | 0

Calculating roots of an equation in Matlab.
Change syms theta a to syms theta

en dag ago | 0

how can i create Bernoulli random variable of a given probability
p = 0.2; n = 20000; A = binornd(ones(n,1),p) p_num = sum(A)/n B = 1 - A one_minus_p_num = sum(B)/n

en dag ago | 0

| accepted

Optimization: Reuse calculations from objective function in nonlinear constraints function
https://de.mathworks.com/help/optim/ug/objective-and-nonlinear-constraints-in-the-same-function.html

en dag ago | 2

| accepted

How apply the conformal mapping? (circle to rectangular)
https://de.mathworks.com/matlabcentral/fileexchange/67117-rectangular-conformal-map

2 dagar ago | 1

I'm having trouble plotting the temperature field for a thin plate with border temperatures known.
Write a function function T = Temp(x,y) T = ...; end in which you calculate the value of the infinite sum for x-y coordin...

2 dagar ago | 1

simultaneously, fitting two functions with two databases
I think you will have to work on the representation of f1 and f2 without the imaginary unit: c10 = 3.0; c210 = 1.0; c220 = -1...

2 dagar ago | 0

| accepted

Intersection of two curves
Do you see a point of intersection ? I don't. N = [1 2 3 4 5 6]; P = [1.48982012272225e-09 1.20413388047633e-08 4.174352140114...

2 dagar ago | 0

Error using SumOperator with Optimization Variables
nTasks = 3; nRess = 2; maxPT = 50; x = optimvar('x',nTasks,nRess,maxPT,'Type', 'integer','LowerBound',0,'UpperBound',1); Con...

3 dagar ago | 0

| accepted

plot based on if and else condition
MATLAB does not accept "double inequalities". Thus instead of if a < b < c you will have to use if a < b && b < c Further, ...

3 dagar ago | 0

Is ode45 always calculate the entire function even if she has been asked only for specific values?
I mean, will it do the full calculation in the background and return only 3 points? Yes, it will do the full calculation in the...

3 dagar ago | 1

| accepted

Constrained dynamic optimization with fmincon
J is minimized in the objective function. The constraint T(x(t_final)) = 0 is set in the function "nonlcon" of fmincon. u (and...

3 dagar ago | 0

| accepted

Where/How exactly do I create an ODE Event Location Function?
if option == 1 % Choosing Straight Run Vs = 10 * (1852 / 3600); % [m/s]-S...

3 dagar ago | 1

How to solve a system of equations using the iteration by assuming the values of the unknowns and check them later. If the error is large return to change the assumed values.
Unknowns are epsilon_ri and zi. Equations are sigma_Ti = 0 and psi_i = psi_i_given. Two equations in two unknowns. Use MATLAB...

3 dagar ago | 0

| accepted

Using fmincon but converging to an infeasible point because the size of the current step is less than the value of the step size tolerance but constraints are not satisfied to
Works for me. T = 1000; %K R = 8.314; %kJ/kmol.K P = 300; %kpa in reformer Po = 100; %kpa standard conditions mair = ...

4 dagar ago | 0

Create the program to return the orthogonal basis and orthonormal basis
As far as I can see, you only return the orthonormal basis from the Gram Schmidt process in the matrix Q. And saving the norms ...

4 dagar ago | 1

| accepted

Computing temperature of a fluid inside a cylinder using PDEPE
For comparison: R0 = .1; R1 = .15; % radii in meters t0 = 0; t1 = 1e4; % time limits wallConductivity = 150; wallDensity =...

4 dagar ago | 0

Fitting experimental data with experimental references
w = [Yls,Yhs,Y3]\Yt if Yls, Yhs, Y3 and Yt are column vectors of the same size.

5 dagar ago | 0

| accepted

How to determine an y value at a specific x value?
sm1 = v1.*N.*L1.*R./(4.*W1.*(L1./C).^2) gives W1 = v1.*N.*L1.*R./(4.*sm1.*(L1./C).^2)

5 dagar ago | 0

isqnonlin: compute part of objective function outside of matlab
Then call "lsqnonlin" as p0 = ...; % initial guess for the vector of material parameters sol = lsqnonlin(@(p)fun(p,y),p0) a...

5 dagar ago | 0

| accepted

Using function handles in a ODE function
x, L, M, V, xlow, xhigh are not given. function y = acceleromet(x) wdwdx = @(t) -6*(x(1)*t.^3 + x(2)*t.^2 + x(3)*t + x(4))*(3*...

5 dagar ago | 0