Beantwortet
Estimating the parameter of a complex equation using maximum likelihood estimation (MLE) method
And this is a probability distribution for arbitrary values of q1, q2, A, lambda, beta1 and beta2 ? I can't believe. If you wer...

fast 3 Jahre vor | 0

Beantwortet
Output issues with deformation calculation
Use double(DEFAB) etc. to see the result in the usual notation.

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
contour plot with a circular boundary
If you want to plot it within a circular contour, you have a 3d-plot. For a 3d-plot, you can only plot cuts in length direction ...

fast 3 Jahre vor | 0

Beantwortet
How do I get largest distance between two rows in an nx3 matrix?
max(pdist(X),[],"all") where X is your matrix.

fast 3 Jahre vor | 0

Beantwortet
interpolation of values gives Inf
InterpCorr.(TitleStr)(1) is equal to 0, so you divide by 0 every time. And my guess is that your line of code should be some ki...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
How to find the value of x if y=0 from two arrays of data?
Use interp1 if y_temperature_soil is monotonic with depth: y0 = interp1(y_temperature_soil,x_depth,0)

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
How do you solve a coupled ODE when one of the ODE results in a vector of length 3 and the other results in a scalar of length 1?
All solution components have to be aggregated in one big vector z, and also the derivatives have to be supplied in this vector f...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
How to get best fitting Parameters for a function for all datapoints
You cannot estimate lh and lo separately in your approximating model, only their sum (lh+lo). But your model seems to return c...

fast 3 Jahre vor | 1

Beantwortet
Linear Interpolation code without using interp1
How can I write a code that works faster? Allocate ynew = zeros(size(n2)) at the beginning of your code. Check whether n is mo...

fast 3 Jahre vor | 1

Beantwortet
creating smaller matrix from a large matrix
You mean rng("default") A = rand(16); for i = 1:4 for j = 1:4 B = A((i-1)*4+1:i*4,(j-1)*4+1:j*4); A_c...

fast 3 Jahre vor | 0

Beantwortet
indexing in the nested loop
sx = 0; sy = 0; for s = T-1:-1:0 sx = sx + x^(T-1-s); sy = sy + y^(T-1-s); ratio(s+1) = sy/(sx+K^(1/gamma)*x^(T...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
nested loop shows only the result of last iteration
... W_t(n+1)=(1+r*(1-tau)).^(-(T-n)).*C_0*epsilon.^T.*K^(1/gamma); for m=(n+1:1:T-1) W_t(n+1)=W_t(n+1)+((C_0.*epsilon.^m...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
I'm not sure as to why MATLAB is not outputting the figure, as well as the vector [T_Euler]. When I run the code it is fine, but when publishing I don't see either. Any fixes?
%% problem description % Determine the temperature vs time for hot iron with convection and % radiation. HW3_prob4_driver()...

fast 3 Jahre vor | 0

Beantwortet
I don't know why this code doesn't stops :(
res(3) == 0 and res(4) == 0 are two nonlinear equations in the unknowns yp(3) and yp(4). ode15i seems to have problems following...

fast 3 Jahre vor | 0

Beantwortet
Delay differential equations where every variable can take on multiple delays
I don't understand your mathematical formulation. Should the summation be over j instead of i ? And shouldn't there be a bracket...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
Solving Multi-variable Partial Differential Equation
The only thing you can get with the "boundary condition" you gave is the solution on the characteristic curve of the PDE through...

fast 3 Jahre vor | 1

| akzeptiert

Beantwortet
Curve fitting with more than 3 variables
Do you really think that there is a causal correlation between the housenumber and the energy from the grid ? So if the house nu...

fast 3 Jahre vor | 1

Beantwortet
i have some error in my MATLAB code .
B*delta_e(t) is undefined since B and delta_e(t) are both 4x1. Maybe you mean B.*delta_e(t) instead.

fast 3 Jahre vor | 0

Beantwortet
How to draw ellipse around my data
https://uk.mathworks.com/help/stats/fitgmdist.html After fitting, you can plot a 90% ellipse around your data. What does it me...

fast 3 Jahre vor | 0

Beantwortet
Event function with DAE system
options = odeset('Mass', M, 'MassSingular', 'yes', 'RelTol', 1e-4, 'AbsTol', 1e-6,'Events',@event); function [position,isterm...

fast 3 Jahre vor | 1

| akzeptiert

Beantwortet
I get error in unrecognized function or variable 'T'.
I changed the bold parts below: T = t; for j = 1:length(t) %with drugs T(j+1)=T(j)+dt; M13(j+1) = M13(j)+...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
Error in ode15s: incompatible sizes arrays, while no errors in function
if Re(i,j)/Re_check(i,j)<0.95 || Re(i,j)/Re_check(i,j)>1.05 fprintf('Re difference too big \n') return end if Re(i...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
How to draw 3d plot with system of equations
Use https://uk.mathworks.com/matlabcentral/fileexchange/9261-plot-2d-3d-region %plotregion([1,2;2,1;1 -1;1 -4;-4,1],[-1;4;-1;-...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
interpolation of coordinates in space using interp3
Given two points P1 = [1 3 4]; and P2 = [3 6 -1]; you can connect them by a line s = @(t) (1-t)*P1 + t*P2; and choose poi...

fast 3 Jahre vor | 0

Beantwortet
How can I plot an ellipsoid using the ellipsoid function with a given equation in the form of ax^2+by^2+cz^2 = 1.
f = @(x,y,z) x.^2 + y.^2 + z.^2 - 1; fimplicit3(f)

fast 3 Jahre vor | 0

Beantwortet
incorrect dimensions in data set
The code works technically. But after all the changes I made, it's not clear if it does what you want. This line needs correcti...

fast 3 Jahre vor | 0

| akzeptiert

Beantwortet
Index in position 1 is invalid. Array indices must be positive integers or logical values. "Ideas of fractal fractional using exponential kernel"
sum(((u(i+1-k,j)-u(i-k,j))/deleta) and k = 1:j For i <= j, the first index in the matrix u becomes <=0.

fast 3 Jahre vor | 0

Beantwortet
Unable to perform assignment because the size of the left side is 5-by-75 and the size of the right side is 5-by-87.
Problem1() function Problem1 %% Parameters values delta1=0.1; Ha=0.1; d1=0.1; d2=0.1; beta=0.1; Pr=1; alpha= pi/4; S1=0.1; ...

fast 3 Jahre vor | 1

| akzeptiert

Beantwortet
Create an array with only the increasing values of a pressure time series.
I prefer a clean for-loop: a_monotone(1) = a(1); for i = 2:numel(a) if a(i) > a_monotone(end) a_monotone(end+1) = a(i)...

fast 3 Jahre vor | 1

Beantwortet
calculate the normal for a plane passing through more than three points
If the points in V1_in.mat all lie in a common plane, you can arbitrarily pick three points and define these points as P. If th...

fast 3 Jahre vor | 1