Bruno Luong

Last seen: Today Active since 2009

Followers: 1   Following: 0

All

Feeds

View by

find a zero of a two-variable function
If you are not unlucky Trajectory = @(x) sum(x.^2,2)-34^2-1; % Test function, use use rather own X0 = [9.609 , 32.288]; %ini...

29 dagen ago | 0

call function name in the same function
This returns the name of the function being executed s = dbstack; funname = s(1).name; In recent MATLAB versions you can do i...

29 dagen ago | 1

| accepted

pinv operation in matlab
NOTE: This code is useful fonly or users who run version greater than R2021a and before R2024a, where pagesvd is supported by pa...

29 dagen ago | 0

How to do 3D matrix multiplication
A = rand(56,32,10); B = rand(32,32,10); C = pagemtimes(A,B); size(C)

29 dagen ago | 0

| accepted

Create an array of distinct objects of handle class
Why not simply calling the constructor as many time as you need arrDir = arrayfun(@(~)uilabel, 1:2); arrDir(1).Text = "textDir...

30 dagen ago | 0

| accepted

How to flip the x axis in a 3d plot?
Do set(gca,'XDir','reverse') one the surfc command is carried out s = peaks'; subplot(1,2,1) surfc(s) xlabel('x') subplot(1...

30 dagen ago | 0

| accepted

Extracting multiple signal segments without a for loop
As you like but this code is perhaps slower than the for-loop signal = rand(1000,1); segmentOnsetIndices = [50 120 550 600 750...

ongeveer een maand ago | 1

| accepted

How to have individual numbers combine into a larger number
One way: UserIput = [1 2 3 10 11]; %% Coding X = polyval(UserIput, 16); X % Combine number clear UserIput % forget the ...

ongeveer een maand ago | 0

Problem with solve() and "Empty sym: 0-by-1".
Your equations look illposed to me for three reasons* For a given j, LHS are known (4 x 1) vector Ax*Ax*B(:,j), let's call ...

ongeveer een maand ago | 0

Efficient Vectorization of For Loop
Not tested but the sign reading tell me M = 2*B' * (A.*C);

ongeveer een maand ago | 1

| accepted

uicontrol doesn't work in Matlab 2023a
The button position you give is way off push the button outside the universe. This happens with all MATLAB. figure;pauseButton=...

ongeveer een maand ago | 0

| accepted

Finding the index k of square nilpotent matrix A
This modified version of binary search only use matrix multiplication and keep track of A^2^p, p = 0,1,2 ...The number of matrix...

ongeveer een maand ago | 0

Finding the index k of square nilpotent matrix A
I modify your original code and now it seems working: function k = nilIndex(A) [n,n]=size(A); l=0; k=n; B =A^l; while k-l>...

ongeveer een maand ago | 0

When Should Generalized eig() Solver be Used for Ordinary eig() Problems?
Actually I think eig with 2 inputs is more approprite for multiple order eigen value and the degenerate case to non trivial jord...

ongeveer een maand ago | 0

| accepted

uniform knot vector for splines
Your knot sequence seems NOT to be suitable for interpolation. The interpolation matrix is singular. x = linspace(1,6,6); xi =...

ongeveer een maand ago | 0

| accepted

Strange behavior with ODE45 dependent on tspan
I would divide the time interval in three parts, as the typical time scales are different during the pulse and outside the pulse...

ongeveer een maand ago | 0

| accepted

Matrix multiplication not returning on certain matrix content (Intel)
The matrix contains many non normalized floating point numbers and operations on these numbers are very slow; >> load('C:\Users...

ongeveer een maand ago | 0

Using fmincon to solve objective function in integral form.
trapz(tspan/length(t), 2*x.^2 + 2*xd.^2) This looks wrong to me, you seem to assume t is equi distance (even in that case you s...

ongeveer een maand ago | 1

| accepted

Lognormal distribution parameters mu and sigma
Is mu=log(mean) and sigma=log(Standard deviation)? No: https://en.wikipedia.org/wiki/Log-normal_distribution If X is log-norma...

ongeveer een maand ago | 0

The best N = 200; % 750 gradX_1 = zeros(N,N); w = rand(N,N); AXB = rand(N,N)+1i*rand(N); A = rand(N,N)+1i*rand(N); B = ran...

ongeveer een maand ago | 0

| accepted

I propose this, and time testing for N = 200; N = 200; % 750 gradX_1 = zeros(N,N); w = rand(N,N); AXB = rand(N,N)+1i*rand(N)...

ongeveer een maand ago | 0

fmincon solution does not differ from initial guess if I provide gradient
To get more info turn the option 'CheckGradients' to 'on', and check the exitflag (third output) of fmincon

ongeveer een maand ago | 1

| accepted

Function matchpairs does not work in MATLAB R2023a when I use a sparse cost matrix. Please help me out ASAP!
Your sparse matrix filled default cost wit 0 not 99. Then your expectation on that the result is identocal for different cost ma...

ongeveer een maand ago | 0

| accepted

error using spline: The first input must contain unique values.
Your spline command wouldn't do what you want, since for every x you have made corresponding two y values. It is NOT a function ...

ongeveer een maand ago | 0

error using spline: The first input must contain unique values.
You did not post your spline command, my guess is you do spline(x,y) Try to remove the last (or the first value of th, then do ...

ongeveer een maand ago | 0

How to fit a sine curve with only the maximum and minimum values
piecewise (co)sine interpolation (not fitting). I don't think you have extra DOF to play with if you restrict data to be local ...

ongeveer een maand ago | 2

symmetric solutions of linear matrix equations
This is a method that use the small "original" linear system. I use pcg since it a linear solver that can accept function handle...

ongeveer een maand ago | 1

symmetric solutions of linear matrix equations
If you have tthe optimization toolbox % Generate test matrices n = 3; A = randn(n); X=rand(n); X = X+X.'; B = X*A; % Add s...

ongeveer een maand ago | 0

| accepted