Beantwortet
How do I edit a loop to be able to store the variables and outputs used in each loop in a way which can be graphed?
Be careful here. In correlations for friction factors, log is usually meant as log to the basis of 10, not the natural logarithm...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Question solve(eqn)
Replace solve (f1, f2, f3, f4, f5, f6) by [a0,a1,a2,a3,a4,a5] = solve (f1, f2, f3, f4, f5, f6) But note that extrapolating t...

mehr als ein Jahr vor | 0

Beantwortet
find is not working in a loop, it skips numbers
Never test two double precision numbers for equality. Because of floating point arithmetic, they will almost always be classifie...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Can someone rectify this code?
My guess is that you work with an older MATLAB version where script and function parts cannot be mixed. Thus save function F =...

mehr als ein Jahr vor | 0

Beantwortet
elseif that goes through previous statements but with different parameters
if e(t)==1 nn1 = nn11; nn2 = nn21; elseif e(t)==2 nn1 = nn12; nn2 = nn22; elseif e(t)==3 nn1 = nn13; nn2 = nn23...

mehr als ein Jahr vor | 0

Beantwortet
syms to another .m file
Call the other .m-file with the syms-equation (or function) as input argument. syms x f = x.^2; result = fun(f) function res...

mehr als ein Jahr vor | 0

Beantwortet
Clarification regarding Sensitivity Analysis of a coupled ODE system using the odesensitivity function
Start with the simple example y'(t) = p(1)*y(t), y(0) = 1. The solution is y(t) = exp(p(1)*t). The (local) sensitivity is ...

mehr als ein Jahr vor | 1

| akzeptiert

Beantwortet
How to make ode45 store/return only results of every n-th time step to reduce memory use?
You can specify the times when the solution should be saved in the input variable "tspan" to "ode45". And you should test wheth...

mehr als ein Jahr vor | 1

| akzeptiert

Beantwortet
Energy Harvesting of Cantilever Beam with MSMA alloy1
I inserted "warning('off')" at the start of your code because your matrix M in your function "system_eqs" is singular and MATLAB...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Sampling from inverse gamma distribution by specifing mean and standard deviation
It works for mu = 0.5, sig = 0.5, e.g. . It seems that your distribution parameter "a" is too close to the critical value of a =...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
How to code an equation with that integrates a vector.
If there is no closed-form expression for u (e.g. if u depends on x), you can solve the vector-differential equation dx/dt = u ...

mehr als ein Jahr vor | 0

Beantwortet
The surface area won't round to one decimal place and I don't know how to get rid of the 'e's.
https://uk.mathworks.com/help/matlab/ref/fprintf.html https://uk.mathworks.com/help/matlab/ref/format.html

mehr als ein Jahr vor | 1

Beantwortet
Need to correct the following script in part of Infinitely Long Horizontal Cylinder (to be of perfect circular base)
%% 2. Infinitely Long Horizontal Cylinder subplot(1,3,2); theta = linspace(0, 2*pi, 30); x_v = linspace(-10, 0, 20); [Theta,...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Grey box model of 1d heat diffusion in a rod (official mathworks example)
I think it will be simpler to couple "pdepe" for the solution of your PDE and "lsqnonlin" for optimizing your model parameters. ...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
trouble using pdepe to solve system of pdes
Your initial condition for u(2) at x = 0 is not consistent with your boundary condition. Use function u0 = icfun(x) global ...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
How plot 2d of more than 3 function and substitute the number in each function separate parameter?
y=@(x,t,m)real(((-0.385e1 .* tanh(-x - 0.310e1) - 0.385e1 * coth(-x - 0.310e1)) .^ (0.1e1 ./ m ./ 0.2e1)) .* exp(i * (0.431e1 * ...

mehr als ein Jahr vor | 1

| akzeptiert

Beantwortet
Help with the Deflection of a Plate in MATLAB
You are trying to solve the heat conduction equation with a homogeneous heat sink -q/D and boundary temperature 0. Is this the s...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Time dependent variable in plot.
The simplest method is to define x, compute y in a loop and plot: x = 0:0.1:15; for i = 1:numel(x) if x(i) >= 0 & x(i) < ...

mehr als ein Jahr vor | 1

| akzeptiert

Beantwortet
Solving Advection Equation PDE with 2nd order Accuracy
Usual second-order schemes for the advection equation will always produce oscillations. Section 3.1.3 of the enclosed article su...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
How can I use fmincon with general (user defined) nonlinear constraints?
g = @(x)deal([gl(x(1:nx),x(nx+1:nx+ny));gf(x(1:nx),x(nx+1:nx+ny))],[])

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Assign fitlme output to a variable
See https://uk.mathworks.com/help/stats/linearmixedmodel.html The functions that can be applied to a fitlme object are: a...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Avoiding -0.0000 as an output
The results for d2,...,d4 aren't exactly 0 because of floating point errors in their computation. So without artificially manipu...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
solve symbolic eigenvalue problem
If you want an analytical expression for the eigenvalues dependent on the parameter a, your matrix A must have dimension <= 4. ...

mehr als ein Jahr vor | 0

Beantwortet
problem to solve incorrrect
%% R-RRR clear all; clc; close all % Input data AB=0.325; %(m) BC=0.938; %(m) CD=0.675; %(m) CE=0.6; %(m) xD=0.8; %(m) y...

mehr als ein Jahr vor | 0

Beantwortet
Monte Carlo integration (hit or miss) to find the area of a circle of radius R
You compare the area of the square with side length R and corner points (0,0), (R,0), (R,R) and (0,R) with the area of the quart...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Finding parameters by fitting data to a system of ODEs with lsqnonlin
syms T1(t) T0(t) s d g eqn1 = diff(T0,t) == s - (d+g) *T0 ; eqn2 = diff(T1,t) == 2*d*T0 + (d-g)*T1; sol = dsolve([eqn1,eqn2])...

mehr als ein Jahr vor | 1

Beantwortet
Can pdepe solve a system of two second-order equations?
Goto https://uk.mathworks.com/matlabcentral/fileexchange/97437-pde1dm download the code, goto folder "documents" and study ex...

mehr als ein Jahr vor | 0

| akzeptiert

Beantwortet
Error in flood_algorithm_rtc_pv (line 22) best_fitness = objective_function(best_solution, V_data, I_data); i don not understand this error
V_data and I_data must be of the same size. This is not the case in your code.

mehr als ein Jahr vor | 1

| akzeptiert

Beantwortet
Issues solving a system of 2, second order, differential equations with ode45
Replace Cd = 0,45 ;% drag coef aircraft [-] rho_Air = 1,2 ...

mehr als ein Jahr vor | 2

| akzeptiert

Beantwortet
simulating rolling 1 and 2 dice
You get a discrete pdf and cdf for the one and two dice-roll experiments you describe. So it doesn't make sense to smoothe resul...

mehr als ein Jahr vor | 0