Community Profile # Torsten

Last seen: Today 2013 以来アクティブ

Professional Interests:

#### Statistics

• • • • • バッジを表示

#### Content Feed

where is the error in my code?
V1 = @(r,w) -acosh(10*(w/(1600*r + 21))^(1/2))/20000000000; V2 = @(r,w) acosh(10*(w/(1600*r + 21))^(1/2))/20000000000; % Defin...

I have three coupled differential equation and need help to solve them
You have a linear system of equations in dx/dt, dy/dt and dz/dt. Solve it for these variables. Then you can write your system e...

Unable to convert the following expression into double array
clc; clear all; close all; al = 0.5; n = 5; sig2_fsnrn = 1; sig2_ksln = 1; sig2_glnrn = 10; sig2_e_fsnrn = 0; sig2...

1日 前 | 0

Can't get the solution on this nonlinear system of equations using vpasolve
For that your two equations can have a solution, one can show that abs((PX^2+PY^2-l1^2-L2^2)/(2*l1*L2)) must be <= 1. This i...

1日 前 | 0

| 採用済み

How to load intermediate variable into workspace while using ode45 ?
tspan = 1:200; x0 = zeros(38,1); [t,x] = ode45(@system, tspan,x0); ts=40; ssActive = [1,2,3,4,1,3,4,2,4,3,3,4,3,3,1,2,3,4,1,...

2日 前 | 0

[f,x]=ecdf(y) definationf of x
https://de.mathworks.com/help/stats/ecdf.html There are six literature reference at the bottom of the page. Usually, the minim...

3日 前 | 0

trying to use the fmincon function
Works for me after changing "fun" in "@fun": % script calling fmincon start clc clear all x0=[-1...

3日 前 | 0

How to generate all possible vectors where each element can take one of two values
Generate all possible combinations of 0 and 1: m = 3; justRows = dec2bin(0:2^m-1)-'0' Now replace the 0's by a(i) and the 1 b...

4日 前 | 1

Local minima of an interpolant fit
Pass the evaluating function for the interpolant as objective function to fminsearch or another optimizer (fmincon,...). E.g. ...

4日 前 | 0

How can plot 1D manifold in 3d space
syms x y z p=0.255*z.^2./(z.^2+(0.055)^2); m1=0.15^2./(p.^2+0.15^2); m2=z.^4./(z.^4+0.1^4); m3=p.^2/(p.^2+0.15^2); g = 0.05...

4日 前 | 0

| 採用済み

How to create a function that can work out an input provided the output and and all other variables are known?
Use "fzero" or "fsolve" to solve G(x) := f(L1,l1,h1) - f(random L2, x, random h2) = 0 for x.

4日 前 | 0

| 採用済み

Workaround for the expinv function in statistics toolbox
p = 0.3; mu = 1/4; x = -mu.*log(1-p) x2 = expinv(p,mu)

4日 前 | 0

| 採用済み

Creating a matrix for a function
A = [ones(numel(t),1),sin(t),cos(t),sin(2*t),cos(2*t)]

4日 前 | 0

| 採用済み

slightly off euler method code
f = @(t,y) exp(t) - t; fSol = @(t) exp(t) - t.^2/2; [t,y]= eulerMethod(f,0.2,2,0,1); plot(t,y) hold on plot(t,fSol...

5日 前 | 0

How to numerically integrate the planar equations of motion?
syms V(t) a1 a2 sol = dsolve(diff(V,t)==a1*V^2+a2) Now substitute a1 = -rho/(2*beta) and a2 = -g0*sin(gamma)

5日 前 | 0

Why am I not getting desired output?
h = 0.05; t = 1:h:500; y0 = [2.5, 3.6, 9.5, 8.8, 2, 3, 7, 7.5]; % t(1) =0; y = zeros(length(t),8); y(1,:) = y0; for i = 1:...

5日 前 | 0

| 採用済み

Why does matlab not evaluate equations properly?
syms x f = log(0.7*x+0.5); x_0 = 4.2857; t0 = taylor(f,x,x_0,'Order',2); h = matlabFunction(t0); h(6.42855) log(0.7*6.4285...

5日 前 | 0

solve a forced vibration equation for x and w
https://de.mathworks.com/help/symbolic/eig.html First multiply by M^(-1) = [1/(3*m) 0; 0 1/m]

5日 前 | 0

recur.m function
Works for me. a = 1:10; %N=10 b = 1:5; %M=4 n = 11:25; x = 11:25; x0 = 3:6; y0 = 1:10; y = recur(a,b,n,x,x0,y0) f...

6日 前 | 0

dsolve gives wrong equations?
There are often several ways to write a function. But it seems to me that your S2 solutions are not independent. syms y(x) C1 ...

7日 前 | 0

| 採用済み

Solving a System of Trig Equations
Use the symbolic toolbox to solve the last equation for cos(theta1). It's a polynomial equation of degree 4 which will have 4 r...

7日 前 | 0

Optimization of Logistic Function Variable K for multiple Inputs Simultaneously
I think you will have to work on the approximation function "fun" in order to get better results. promAngles = [146.7589 115.7...

7日 前 | 0

| 採用済み

How to stop ODE fuction in a certain value
Better you use the Event facility of the ODE integrators. Note that your results become complex-valued because y(4) becomes gre...

8日 前 | 0

How to solve and store a system of non-linear equations with an array of outputs?
Taking cotd on both sides of the second equation gives cotd(dRF)-cotd(dLF)=Ap*w/L; cotd(dRF)+cotd(dLF)=2*cotd(davg(i)); Thus ...

9日 前 | 1

| 採用済み

Optimize transition rate matrix
Instead of a 4x4 matrix use a (16x1) vector. Then it is easy to implement your constraints in Aeq and beq. If somewhere in your...

9日 前 | 0

| 採用済み

Suggestions on solving this complex-vector-nonlinear system of equations in MATLAB
This is a simple linear system of equations in the unknowns W and Z: alpha2 = 205*(pi/180); alpha3 = 55*(pi/180); beta2 = 39*...

10日 前 | 0

| 採用済み

Intersection of two curves in matlab
clc %crate=0:1:10; G=-I_BT; h=4000;%no. of cycles T=30+273;%temp in kelvin R=8.314;%universal gas const capa=200;%battery ...

10日 前 | 0

Polyfit and Polyval on Reduced Matrices
You are aware that you fit parameters a and b of the function f(x) = a + 10*b*log10(x) to your data vector y ? Did you try f...

10日 前 | 0

Attempting to Plot In A Specific Range
x1 = (linspace(-3.5,3.5,numel(y))).'; plot(x1,y,'-')

10日 前 | 0

Itterate over array in steps of 120.
r=readmatrix('data_jong_edit.txt'); r(isnan(r(:,1)),:)=[]; empty = nan(16800,2); r = [r empty]; u=unique(r(:,1)); for k=1...

10日 前 | 0

| 採用済み