Risposto
Split non-bijective function into two bijective regions
x = [4;2;3;6;6.5;4.5;1;1.5;3.5;2.5;5;6.2;2.8;6.3]; y = [5.5;6.5;1;4;1.5;1;2;7;6;1.8;5;1;6.2;3]; % TSP to order the data, See...

oltre 2 anni fa | 1

| accettato

Risposto
Why is there a difference in computing least squares although same symbolic operator \ is used?
"I'm guessing for my case only the documentation of page * matters? " Yes For numerical case and vector of RHS b x = A\b is ...

oltre 2 anni fa | 0

| accettato

Risposto
Does parfor inside lsqnonlin with UseParallel=True work properly?
If your function use parfor internally then you should NOT set UseParallel to true.

oltre 2 anni fa | 0

Risposto
Help finding datapoints in x,y falling outside the boundaries of a polygon. Incongruencies between inpolygon and inpoly2, time issues, and possible solutions
The reason that inpoly2/insidepoly fail because the polygons p1 and p3 have holes, therefore px.Vertices have NaN values. Here ...

oltre 2 anni fa | 1

| accettato

Risposto
generate a plane that has the coordinates of a 3D matrix as its boundary
Why don't you do the very obvious? load('C.mat') patch(C(:,1),C(:,2),C(:,3),'k') hold on plot3(C(:,1),C(:,2),C(:,3),'g.','Ma...

oltre 2 anni fa | 0

| accettato

Risposto
can i velocize it? (if possible to vectorize)
On my computer the acceleratoion of Method V2 is more impressive load('matlab_DD.mat') load('matlab_minTrades.mat') load('m...

oltre 2 anni fa | 0

| accettato

Risposto
clc command window restore
From MATLAB windows go to the menu Layout -> Command History -> Dock. You then have a tab "Command History" appears on the right...

oltre 2 anni fa | 0

Risposto
How can I log commands being executed by matlab? Like diary() but for all function calls.
Check out profile see if it meets your need.

oltre 2 anni fa | 1

Risposto
calculate the normal for a plane passing through more than three points
load('V1_in.mat') [~,~,V]=svd( V_1-mean(V_1)); N=V(:,3)

oltre 2 anni fa | 0

| accettato

Risposto
Changing default encoding of MATLAB
I see someone succeed in older MATLAB by editing the file [matlabroot '\bin\lcdata.xml'] Ref in this thread https://www.math...

oltre 2 anni fa | 0

| accettato

Risposto
Dynamic Assign Field Access in MATLAB Structures
clear s=struct('a',1,'b',2,'c',3) % this is a dirty command that pops the field values to the workspace cellfun(@(f) assignin...

oltre 2 anni fa | 0

Risposto
Skipping iteration on child throwing a warning
Modification of script_Z for k = 1:200 lastwarn('',''); % clean the warning state x = input_postprocessed_1/input_pos...

oltre 2 anni fa | 1

| accettato

Risposto
(Matlab Coder)Domain error. To compute complex results from real x, use 'acosd(complex(x))'
"Obviously, the a,b I gave are equal and the val value should be equal to 1" wrong. Before claming anything, you should first c...

oltre 2 anni fa | 0

| accettato

Risposto
Compute group statistics when the group definition is an OR of values found in several columns of a table
I submit a method without loop. I expected to be faster but to my surprise it is slower than v2. But at least it is scalable. %...

oltre 2 anni fa | 0

| accettato

Risposto
Why I am getting "Matrix is singular to working precision"?
The warning message is there "Matrix is singular to working precision" because your equations are dependent % 2*w+3*x+y-11*z =...

oltre 2 anni fa | 0

Risposto
how to effectively find a starting point for fzero?
The problem is not finding the right starting point, the problem is that your function func implementation returns NaN for large...

oltre 2 anni fa | 0

Risposto
How to estimate 95% probability contour in bivariate kernel density function MATLAB
I'm not sure if there is a standard definition of area of 95 confidence in 2D pdf since it like asking what if the shape of 95% ...

oltre 2 anni fa | 0

Risposto
How to keep the values for particular dimension of a matrix and rest of the dimension will be zeros?
A single line solution A(~(ismember(1:size(A,2),350:400) & ismember((1:size(A,1))',20:35))) = 0

oltre 2 anni fa | 0

Risposto
Finding value in Matrix to corresponding minvalue Position in different Matrix
A=randi(10,5) B=randi(10,5) [Amin,jmin]=min(A,[],2) % Amin is A(sub2ind(size(A),(1:size(A,1))',jmin)) % Bmin = B(sub2ind(si...

oltre 2 anni fa | 0

| accettato

Risposto
generate new coordinates (starting from initial coordinates) that are arranged outwards (by a distance H) and on the same plane
load('trace.mat') mt = mean(trace); dt = trace-mt; [~,~,V]=svd(dt); Q=V(:,[1:2]); t2D = dt*Q; [~,is]=sort(atan2(t2D(:,...

oltre 2 anni fa | 0

| accettato

Risposto
What will be a suitable option to circumvent a non-invertible matrix?
When you encount warning message like this, the first thing to ask is not "how do I change linear inversion algorithm?", but "wh...

oltre 2 anni fa | 0

Risposto
Determine the intersection coordinates between a 3D object and an infinite plane
Mesh your plane then use this https://www.mathworks.com/matlabcentral/fileexchange/49160-fast-mesh-mesh-intersection-using-ray-...

oltre 2 anni fa | 0

| accettato

Risposto
How to tune matrices to specific output?
Use lsqnonlin (or such) if you have optimization toolbox.

oltre 2 anni fa | 1

| accettato

Risposto
Converged to an infeasible point
I interpret fmincon message for you : Your system is not really controllable and the solution of the final step is driven by the...

oltre 2 anni fa | 0

Risposto
I will find an xy dataset satisfying an implicit equation.
f = @(x,y)(x.^2 + y.^2).^3 + (15.*x + 3.3.*y).*(x.^2 + y.^2).^2 + (62.5.*x.^2 + 20.*x.*y - 8.37.*y.^2).*(x.^2 + y.^2) - 17.*x.^3...

oltre 2 anni fa | 0

| accettato

Risposto
Problems with savepath in R2023b
On Windows to fix it you need to change the Security setup of the file pathdef.m Right Click on file, Properties, Tab Security, ...

oltre 2 anni fa | 2

Risposto
generate triangular mesh starting from the external nodes (rx3) of a circle
load V; n=size(V,1); % 360 stlwrite(triangulation([ones(n-2,1) (2:n-1)' (3:n)'],V),'mycircle.stl')

oltre 2 anni fa | 1

Risposto
LU factorization with decreasing elements on the main diagonal of U
Obviously not A=[1 10 9; 5 1 9; 2 8 1] [L,U,P]=lu(A) But you can fix the progression of abs(diag(U)) in any arbitra...

oltre 2 anni fa | 1

Risposto
How to solve out of memory error when constructing sparse matrix
You might try the syntax nz = ... ; < 222263280 S = sparse(i,j,v,m,n,nz); with nz is the estimation of non-zero elements of y...

oltre 2 anni fa | 0

| accettato

Risposto
ascribe zeros to an array
To pad a vector with zeros you only need to set the last-element of the extended vector to 0 f=rand(1,5) % 1000 in your case ...

oltre 2 anni fa | 0

Carica altro