photo

Torsten


Last seen: Today |  Actif depuis 2013

Followers: 20   Following: 0

Statistiques

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

Afficher les badges

Feeds

Afficher par

Réponse apportée
how to use fmincon to optimize a functio with more than one control vector where the control vectors used in constrains also?
If all three vectors/matrices are to be treated as solution variables for "fmincon", use v = [p(:);X(:);L(:)] It's useful to t...

environ 7 heures il y a | 0

Réponse apportée
Double Integration of Numeric Data - Acceleration -> Displacement
Assuming x(0) = v(0) = 0 in both cases, I get load("A2.mat") load("A9.mat") a2 = Datashort1; a9 = Datashort21; figure plot...

3 jours il y a | 0

Réponse apportée
Find real valued parameters of a complex equation/model by using optimization or curve fitting.
I don't know if this helps to get better results, but of course you have to fit real and imaginary part of mu_eff simultaneously...

4 jours il y a | 0

| A accepté

Réponse apportée
Storing while loop results in an array and plotting
Why do you think the results you get from the second part of your code are wrong ? g = 10; y(1) = 100; v(1) = 0; t(1) = 0; ...

6 jours il y a | 0

| A accepté

Réponse apportée
how to solve this error? Objective function is returning undefined values at initial point. FSOLVE cannot continue.
In the computation of c0,...,c10, you divide by x(3),...,x(13), thus by 0 because x0 = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0] . After ...

6 jours il y a | 0

Réponse apportée
4th order Runge - Kutta vs. 4th order Predictor - Corrector
%--- Exercise 1 - A %--- Solving the diff. eqn. N'(t) = N(t) - c*N^2(t) with c = 0.5 using %--- Euler's method, Runge - Kutta ...

9 jours il y a | 0

| A accepté

Réponse apportée
4th order Runge - Kutta vs. 4th order Predictor - Corrector
For k>4, the variable RKN should no longer appear in the equations PredN(k) = PredN(k-1) + h/24*( 55*f(t(k-1),RKN(k-1)...

9 jours il y a | 0

Réponse apportée
Runge - kutta 4th order method for two different steps
y(j+1) = y(j) + h(i)*(K1 + K4 + 2*(K2 + K3))/6; instead of y(j+1) = y(j) + (K1 + K4 + 2*(K2 + K3))/6;

11 jours il y a | 0

| A accepté

Réponse apportée
Error using surf (line 71) Z must be a matrix, not a scalar or vector.
Works for me (assuming that FL_hat is 1001x1001 in your case): alpha = linspace(0, 2*pi, 51); beta = linspace(0, 2*pi,21); FL...

12 jours il y a | 0

Réponse apportée
How to avoid the warnings on my code
%--- Solving the diff. eqn. y'(t) = y(t) - c*y^2(t) with c = 0.5 using %--- Euler's method, Runge - Kutta method, and predictor...

12 jours il y a | 1

| A accepté

Réponse apportée
I am new to using Matlab, and wanted to ask if there is a way to group datapoints into bins/columns by numerical range? Thanks
https://uk.mathworks.com/help/matlab/ref/double.discretize.html https://uk.mathworks.com/help/matlab/ref/matlab.graphics.chart....

15 jours il y a | 0

| A accepté

Réponse apportée
Plugging a new value to ksdensity
rng("default") r = randn(100,1); x_new = 1.1; f = ksdensity(r,x_new) or rng("default") r = randn(100,1); f = @(x)ksdensit...

15 jours il y a | 0

| A accepté

Réponse apportée
which software better to learn?
If you ask this in a MATLAB forum, you will get a different answer than in a python forum :-) If you start learning about optim...

15 jours il y a | 1

Réponse apportée
Hello everyone, I have a question about numerical integration. The formula is shown below, where the integration path C_beta is a closed curve on the complex plane.
However, as mentioned in my question, the array list_beta forms a closed curve on the complex plane. Isn't then dq = q' dbeta a...

17 jours il y a | 0

| A accepté

Réponse apportée
graph is not showing
P=1; % assuming and constant Ex=20*10^4; Ey=2*10^4; Ea=0.1*10^4; k=0.1; t=0.01; b=0.05; h=0.1; ExGxy=17; Ah=((Ex/Ey)^0....

21 jours il y a | 0

Réponse apportée
Diffusion equation in Cartesian and Cylindrical coordinates
If m = 1, the condition at r = 0 set by "pdepe" is pl = 0, ql = 1. All other conditions make no sense mathematically because the...

23 jours il y a | 0

Réponse apportée
"Too Many Input/Output Arguments" error for nrlmsise00 function
The results for the defined testcase seem to coincide with what is stated in the Documentation (see below). You will have to ad...

24 jours il y a | 1

| A accepté

Réponse apportée
Weighted regression without priority zone, but with finding the minimum possible value for all the differences
x = [1 2 3 4 5 10 20 50].'; y = [5 10 15 20 25 30 40 50 60 90]; z = [265.110074796797 195.008035083193 156.258229529605 129.74...

25 jours il y a | 1

| A accepté

Réponse apportée
Code Matlab for a function that depends in a real parameter
integrand = @(f,y)y.*sin(f(y)); F_alpha = @(f,alpha,x)alpha.*x.*integral(@(y)integrand(f,y),0,1); f = @(u)u.^2; alpha = 1; x...

25 jours il y a | 0

Réponse apportée
Weighted regression without priority zone, but with finding the minimum possible value for all the differences
Use "fminimax". It minimizes the maximum error and is the nonlinear equivalent to minimizing the Inf-norm in linear regression p...

26 jours il y a | 2

Réponse apportée
Why am I getting error using assignin and syms while solving pde using pdepe
Here is the code for the dimensional model. It should be no problem to rewrite it in the non-dimensional form. Note that the s...

26 jours il y a | 0

| A accepté

Réponse apportée
Why am I getting error using assignin and syms while solving pde using pdepe
tau is a vector of values in your code: tau = linspace(0, tau_end, 100); % Non-dimensional time You can't use it as a symbolic...

27 jours il y a | 0

Réponse apportée
Andrica's Conjecture
Use nthprime(1:n+1) to generate the vector of the first n+1 prime numbers. The rest should be simple.

27 jours il y a | 0

| A accepté

Réponse apportée
Forward, Reverse finite difference question
x = 0:0.1:1; f = @(x)x.^2; derf = (f(x(2:end))-f(x(1:end-1)))/0.1; hold on plot(x(2:end),derf); % backward plot(x(1:end-...

27 jours il y a | 1

| A accepté

Réponse apportée
Time interation error while trying to solve a partial differential equation which describes pore diffusion resistance combined with with surface kinetics.
One error, one problem: The error: qr = De; must be replaced by qr = 1; because the boundary condition is given by p+q*f =...

28 jours il y a | 0

| A accepté

Réponse apportée
Numerical differentiation and integration
t = [0,5:10:115]; s = [0, 4892.3, 14041.3, 22371.3, 29926.1, 36764.4, 42965.4, 48634.8, 53910.3, 58961.8, 63981.6, 69165.5, 74...

30 jours il y a | 0

| A accepté

Réponse apportée
error in runge kutta -Unable to perform assignment because the left and right sides have a different number of elements.
You always have to use the variables at the ith time step, not the complete array in the assignments. I didn't check your compl...

environ un mois il y a | 0

Réponse apportée
How to plot the temperature unit which is small circle raised at the left of C upper case?
Maybe by adding the interpreter ? In R2024b, your title is correctly plotted without any modifications. x=0:0.1:1; y=x.^2; pl...

environ un mois il y a | 0

| A accepté

Réponse apportée
The method of how to calculate the FWHM.
xmax = 0; ymax = 4; x = xmax:0.01:3; f = -x.^2 + ymax; FWHM = 2*interp1(f,x,ymax/2) Or directly use the inverse function of...

environ un mois il y a | 0

Réponse apportée
How to make my code run faster
Maybe your system does not have a solution. sol = fsolve(@fun,[1 1 0.1]); fun(sol) sol = 1./sol function equations = fun(par...

environ un mois il y a | 0

Charger plus