# Bruno Luong

Last seen: Today 2009 年からアクティブ

Followers: 1   Following: 0

All

バッジを表示

#### Feeds

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...

2ヶ月 前 | 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...

2ヶ月 前 | 1

| 採用済み

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...

2ヶ月 前 | 0

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

2ヶ月 前 | 0

| 採用済み

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...

2ヶ月 前 | 0

| 採用済み

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...

2ヶ月 前 | 0

| 採用済み

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...

2ヶ月 前 | 1

| 採用済み

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 ...

2ヶ月 前 | 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 ...

2ヶ月 前 | 0

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

2ヶ月 前 | 1

| 採用済み

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=...

2ヶ月 前 | 0

| 採用済み

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...

2ヶ月 前 | 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>...

2ヶ月 前 | 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...

2ヶ月 前 | 0

| 採用済み

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 =...

2ヶ月 前 | 0

| 採用済み

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...

2ヶ月 前 | 0

| 採用済み

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...

2ヶ月 前 | 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...

2ヶ月 前 | 1

| 採用済み

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...

2ヶ月 前 | 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...

2ヶ月 前 | 0

| 採用済み

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)...

2ヶ月 前 | 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

3ヶ月 前 | 1

| 採用済み

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...

3ヶ月 前 | 0

| 採用済み

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 ...

3ヶ月 前 | 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 ...

3ヶ月 前 | 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 ...

3ヶ月 前 | 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...

3ヶ月 前 | 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...

3ヶ月 前 | 0

| 採用済み

how to multiply every single element in a matrix to entire of another matrix?
x = [100,500,900,1300,1700; 120,600,1080,1560,2040; 140,700,1260,1820,2380; 160,800,1440,2080,2720; 180,900,...

3ヶ月 前 | 1

symmetric solutions of linear matrix equations
% Generate test matrices n = 5; A = randn(n); X=rand(n); X = X+X.'; B = A*X; % Add small noise to B B = B + 1e-1*randn(siz...

3ヶ月 前 | 0