Hello,
this is more of a theory question but I would greatly appreciate if you could cite good sources with the theory behind my following question:
Why do we have to define (either numerically or analytically) the jacobians for the NMPC in Matlab?
I can understand that the Jacobian for the A, B and C matrices are for the linearization of the prediction model, but I don’t really understand why we need Jacobians for the Cost function and the Constraints as well.Hello,
this is more of a theory question but I would greatly appreciate if you could cite good sources with the theory behind my following question:
Why do we have to define (either numerically or analytically) the jacobians for the NMPC in Matlab?
I can understand that the Jacobian for the A, B and C matrices are for the linearization of the prediction model, but I don’t really understand why we need Jacobians for the Cost function and the Constraints as well. Hello,
this is more of a theory question but I would greatly appreciate if you could cite good sources with the theory behind my following question:
Why do we have to define (either numerically or analytically) the jacobians for the NMPC in Matlab?
I can understand that the Jacobian for the A, B and C matrices are for the linearization of the prediction model, but I don’t really understand why we need Jacobians for the Cost function and the Constraints as well. nmpc MATLAB Answers — New Questions

​

i plot the real one but when i change the plot for abs i can get the abs why i can’t get that of this function
0.2e1 * (0.4e1 * exp((-8 * t + 2 * x)) + 0.9e1 * exp((-27 * t + 3 * x)) + 0.25e2 * exp((-35 * t + 5 * x + 2))) / (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) – 0.2e1 * (0.2e1 * exp((-8 * t + 2 * x)) + 0.3e1 * exp((-27 * t + 3 * x)) + 0.5e1 * exp((-35 * t + 5 * x + 2))) ^ 2 / (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) ^ 2
and my code for ploting is
abs(0.2e1 .* (0.4e1 * exp((-8 * t + 2 * x)) + 0.9e1 .* exp((-27 * t + 3 * x)) + 0.25e2 .* exp((-35 * t + 5 * x + 2))) ./ (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) – 0.2e1 * (0.2e1 .* exp((-8 * t + 2 * x)) + 0.3e1 .* exp((-27 * t + 3 * x)) + 0.5e1 .* exp((-35 * t + 5 * x + 2))) .^ 2 ./ (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 .* x)) + exp((-35 * t + 5 * x + 2))) .^ 2).^2i plot the real one but when i change the plot for abs i can get the abs why i can’t get that of this function
0.2e1 * (0.4e1 * exp((-8 * t + 2 * x)) + 0.9e1 * exp((-27 * t + 3 * x)) + 0.25e2 * exp((-35 * t + 5 * x + 2))) / (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) – 0.2e1 * (0.2e1 * exp((-8 * t + 2 * x)) + 0.3e1 * exp((-27 * t + 3 * x)) + 0.5e1 * exp((-35 * t + 5 * x + 2))) ^ 2 / (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) ^ 2
and my code for ploting is
abs(0.2e1 .* (0.4e1 * exp((-8 * t + 2 * x)) + 0.9e1 .* exp((-27 * t + 3 * x)) + 0.25e2 .* exp((-35 * t + 5 * x + 2))) ./ (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) – 0.2e1 * (0.2e1 .* exp((-8 * t + 2 * x)) + 0.3e1 .* exp((-27 * t + 3 * x)) + 0.5e1 .* exp((-35 * t + 5 * x + 2))) .^ 2 ./ (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 .* x)) + exp((-35 * t + 5 * x + 2))) .^ 2).^2 i plot the real one but when i change the plot for abs i can get the abs why i can’t get that of this function
0.2e1 * (0.4e1 * exp((-8 * t + 2 * x)) + 0.9e1 * exp((-27 * t + 3 * x)) + 0.25e2 * exp((-35 * t + 5 * x + 2))) / (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) – 0.2e1 * (0.2e1 * exp((-8 * t + 2 * x)) + 0.3e1 * exp((-27 * t + 3 * x)) + 0.5e1 * exp((-35 * t + 5 * x + 2))) ^ 2 / (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) ^ 2
and my code for ploting is
abs(0.2e1 .* (0.4e1 * exp((-8 * t + 2 * x)) + 0.9e1 .* exp((-27 * t + 3 * x)) + 0.25e2 .* exp((-35 * t + 5 * x + 2))) ./ (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 * x)) + exp((-35 * t + 5 * x + 2))) – 0.2e1 * (0.2e1 .* exp((-8 * t + 2 * x)) + 0.3e1 .* exp((-27 * t + 3 * x)) + 0.5e1 .* exp((-35 * t + 5 * x + 2))) .^ 2 ./ (exp(0.2e1) + exp((-8 * t + 2 * x)) + exp((-27 * t + 3 .* x)) + exp((-35 * t + 5 * x + 2))) .^ 2).^2 plotting, 3d plots MATLAB Answers — New Questions

​

I have a text document, called transition.txt, that has this inside it.
L,R,D’,B,U,D,F’,R,B’,L’,B2,R,F2,D2,F,L’,U2,L,U2,L,B2,L’,B2,L2,U2,L2,R2,U2,B2,U2
I need to send it to my arduino. Not sure what the code should do on the Matlab side. Also, how would I make sure the information got to the arduino in the way I wanted it to?

Here is my arduino code.
const size_t maxCommandNumber = 50;
const size_t maxCommandSize = 3 * maxCommandNumber; // if commands fit on 2 characters and a comma, this is enough for maxCommandNumber commands
char receivedChars[maxCommandSize + 1]; // +1 for trailing null char required by cStrings
const char endMarker = ‘n’;

void func_L(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Lp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_L2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_R(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Rp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_R2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_D(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Dp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_D2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_B(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Bp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_B2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_U(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Up(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_U2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_F(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Fp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_F2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_unknown(const char* s) {Serial.print(F("Unknown commmand: ")); Serial.println(s);}

struct {
const char* label;
void (*func)(const char*);
} keys[] = {
{"L", func_L},
{"L’", func_Lp},
{"L2", func_L2},
{"R", func_R},
{"R’", func_Rp},
{"R2", func_R2},
{"D", func_D},
{"D’", func_Dp},
{"D2", func_D2},
{"B", func_B},
{"B’", func_Bp},
{"B2", func_B2},
{"U", func_U},
{"U’", func_Up},
{"U2", func_U2},
{"F", func_F},
{"F’", func_Fp},
{"F2", func_F2},
};
const size_t nbKeys = sizeof keys / sizeof keys[0];

bool findAndExecute(const char* s)
{
bool success = false;
for (size_t i = 0; i < nbKeys; i++)
if (!strcmp(s, keys[i].label)) {
keys[i].func(s);
success = true;
break;
}
if (!success) func_unknown(s);
return success;
}

bool commandAvailable()
{
static size_t commandIndex = 0;
static bool awaitEndMarkerError = false;

int incomingByte = Serial.read();
bool commandComplete = false;
if (incomingByte != -1) { // read returns -1 if there is nothing to read (faster than using available)
if (awaitEndMarkerError) {
if (incomingByte == endMarker) {
awaitEndMarkerError = false;
commandIndex = 0;
receivedChars[0] = ”; // discard what was there
commandComplete = true; // if we want to still notify the parser we got a command (here will be empty)
}
} else {
char incomingCharacter = incomingByte; // grab LSB as a char
if (incomingByte != endMarker) {
if (commandIndex < maxCommandSize) {
if (!isspace(incomingCharacter)) { // discard all spaces. could also add a test to only keep A-Z and ‘ and 2, whatever is legit in the command language
receivedChars[commandIndex++] = incomingCharacter;
receivedChars[commandIndex] = ”; // always maintain a correct cString
}
} else {
Serial.println(F("ERROR: Command line too big. Discarded.nt-> increase maxCommandNumber or send shorter command lines."));
awaitEndMarkerError = true;
}
} else {
commandIndex = 0;
commandComplete = true;
}
}
}
return commandComplete;
}

bool handleCommand()
{
bool success = true;
const char* separators = ",";
Serial.print(F("nparsing : ")); Serial.println(receivedChars);
char* ptr = strtok(receivedChars, separators);
if (*ptr == ”) Serial.println(F("Error: empty command"));
else
while (ptr) {
success &= findAndExecute(ptr);
ptr = strtok(NULL, separators);
}
return success;
}

void setup() {
Serial.begin(115200);
}

void loop() {
if (commandAvailable())
if (! handleCommand())
Serial.println(F("Some commands were not recognized"));
}I have a text document, called transition.txt, that has this inside it.
L,R,D’,B,U,D,F’,R,B’,L’,B2,R,F2,D2,F,L’,U2,L,U2,L,B2,L’,B2,L2,U2,L2,R2,U2,B2,U2
I need to send it to my arduino. Not sure what the code should do on the Matlab side. Also, how would I make sure the information got to the arduino in the way I wanted it to?

Here is my arduino code.
const size_t maxCommandNumber = 50;
const size_t maxCommandSize = 3 * maxCommandNumber; // if commands fit on 2 characters and a comma, this is enough for maxCommandNumber commands
char receivedChars[maxCommandSize + 1]; // +1 for trailing null char required by cStrings
const char endMarker = ‘n’;

void func_L(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Lp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_L2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_R(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Rp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_R2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_D(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Dp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_D2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_B(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Bp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_B2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_U(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Up(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_U2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_F(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Fp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_F2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_unknown(const char* s) {Serial.print(F("Unknown commmand: ")); Serial.println(s);}

struct {
const char* label;
void (*func)(const char*);
} keys[] = {
{"L", func_L},
{"L’", func_Lp},
{"L2", func_L2},
{"R", func_R},
{"R’", func_Rp},
{"R2", func_R2},
{"D", func_D},
{"D’", func_Dp},
{"D2", func_D2},
{"B", func_B},
{"B’", func_Bp},
{"B2", func_B2},
{"U", func_U},
{"U’", func_Up},
{"U2", func_U2},
{"F", func_F},
{"F’", func_Fp},
{"F2", func_F2},
};
const size_t nbKeys = sizeof keys / sizeof keys[0];

bool findAndExecute(const char* s)
{
bool success = false;
for (size_t i = 0; i < nbKeys; i++)
if (!strcmp(s, keys[i].label)) {
keys[i].func(s);
success = true;
break;
}
if (!success) func_unknown(s);
return success;
}

bool commandAvailable()
{
static size_t commandIndex = 0;
static bool awaitEndMarkerError = false;

int incomingByte = Serial.read();
bool commandComplete = false;
if (incomingByte != -1) { // read returns -1 if there is nothing to read (faster than using available)
if (awaitEndMarkerError) {
if (incomingByte == endMarker) {
awaitEndMarkerError = false;
commandIndex = 0;
receivedChars[0] = ”; // discard what was there
commandComplete = true; // if we want to still notify the parser we got a command (here will be empty)
}
} else {
char incomingCharacter = incomingByte; // grab LSB as a char
if (incomingByte != endMarker) {
if (commandIndex < maxCommandSize) {
if (!isspace(incomingCharacter)) { // discard all spaces. could also add a test to only keep A-Z and ‘ and 2, whatever is legit in the command language
receivedChars[commandIndex++] = incomingCharacter;
receivedChars[commandIndex] = ”; // always maintain a correct cString
}
} else {
Serial.println(F("ERROR: Command line too big. Discarded.nt-> increase maxCommandNumber or send shorter command lines."));
awaitEndMarkerError = true;
}
} else {
commandIndex = 0;
commandComplete = true;
}
}
}
return commandComplete;
}

bool handleCommand()
{
bool success = true;
const char* separators = ",";
Serial.print(F("nparsing : ")); Serial.println(receivedChars);
char* ptr = strtok(receivedChars, separators);
if (*ptr == ”) Serial.println(F("Error: empty command"));
else
while (ptr) {
success &= findAndExecute(ptr);
ptr = strtok(NULL, separators);
}
return success;
}

void setup() {
Serial.begin(115200);
}

void loop() {
if (commandAvailable())
if (! handleCommand())
Serial.println(F("Some commands were not recognized"));
} I have a text document, called transition.txt, that has this inside it.
L,R,D’,B,U,D,F’,R,B’,L’,B2,R,F2,D2,F,L’,U2,L,U2,L,B2,L’,B2,L2,U2,L2,R2,U2,B2,U2
I need to send it to my arduino. Not sure what the code should do on the Matlab side. Also, how would I make sure the information got to the arduino in the way I wanted it to?

Here is my arduino code.
const size_t maxCommandNumber = 50;
const size_t maxCommandSize = 3 * maxCommandNumber; // if commands fit on 2 characters and a comma, this is enough for maxCommandNumber commands
char receivedChars[maxCommandSize + 1]; // +1 for trailing null char required by cStrings
const char endMarker = ‘n’;

void func_L(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Lp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_L2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_R(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Rp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_R2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_D(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Dp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_D2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_B(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Bp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_B2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_U(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Up(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_U2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_F(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_Fp(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_F2(const char* s) {Serial.print(F("in Function ")); Serial.println(s);}
void func_unknown(const char* s) {Serial.print(F("Unknown commmand: ")); Serial.println(s);}

struct {
const char* label;
void (*func)(const char*);
} keys[] = {
{"L", func_L},
{"L’", func_Lp},
{"L2", func_L2},
{"R", func_R},
{"R’", func_Rp},
{"R2", func_R2},
{"D", func_D},
{"D’", func_Dp},
{"D2", func_D2},
{"B", func_B},
{"B’", func_Bp},
{"B2", func_B2},
{"U", func_U},
{"U’", func_Up},
{"U2", func_U2},
{"F", func_F},
{"F’", func_Fp},
{"F2", func_F2},
};
const size_t nbKeys = sizeof keys / sizeof keys[0];

bool findAndExecute(const char* s)
{
bool success = false;
for (size_t i = 0; i < nbKeys; i++)
if (!strcmp(s, keys[i].label)) {
keys[i].func(s);
success = true;
break;
}
if (!success) func_unknown(s);
return success;
}

bool commandAvailable()
{
static size_t commandIndex = 0;
static bool awaitEndMarkerError = false;

int incomingByte = Serial.read();
bool commandComplete = false;
if (incomingByte != -1) { // read returns -1 if there is nothing to read (faster than using available)
if (awaitEndMarkerError) {
if (incomingByte == endMarker) {
awaitEndMarkerError = false;
commandIndex = 0;
receivedChars[0] = ”; // discard what was there
commandComplete = true; // if we want to still notify the parser we got a command (here will be empty)
}
} else {
char incomingCharacter = incomingByte; // grab LSB as a char
if (incomingByte != endMarker) {
if (commandIndex < maxCommandSize) {
if (!isspace(incomingCharacter)) { // discard all spaces. could also add a test to only keep A-Z and ‘ and 2, whatever is legit in the command language
receivedChars[commandIndex++] = incomingCharacter;
receivedChars[commandIndex] = ”; // always maintain a correct cString
}
} else {
Serial.println(F("ERROR: Command line too big. Discarded.nt-> increase maxCommandNumber or send shorter command lines."));
awaitEndMarkerError = true;
}
} else {
commandIndex = 0;
commandComplete = true;
}
}
}
return commandComplete;
}

bool handleCommand()
{
bool success = true;
const char* separators = ",";
Serial.print(F("nparsing : ")); Serial.println(receivedChars);
char* ptr = strtok(receivedChars, separators);
if (*ptr == ”) Serial.println(F("Error: empty command"));
else
while (ptr) {
success &= findAndExecute(ptr);
ptr = strtok(NULL, separators);
}
return success;
}

void setup() {
Serial.begin(115200);
}

void loop() {
if (commandAvailable())
if (! handleCommand())
Serial.println(F("Some commands were not recognized"));
} arduino, serial, send MATLAB Answers — New Questions

​

For example if I were to type this in matlab:
>> 1/[1, 2 ,3]
There would be an error:
>> Error using /
Matrix dimensions must agree.

But if I were to type this in Matlab
>> 1/[1; 2; 3]
The answer would be:
[0 0 0.3333]

Why is Matlab doing this ?For example if I were to type this in matlab:
>> 1/[1, 2 ,3]
There would be an error:
>> Error using /
Matrix dimensions must agree.

But if I were to type this in Matlab
>> 1/[1; 2; 3]
The answer would be:
[0 0 0.3333]

Why is Matlab doing this ? For example if I were to type this in matlab:
>> 1/[1, 2 ,3]
There would be an error:
>> Error using /
Matrix dimensions must agree.

But if I were to type this in Matlab
>> 1/[1; 2; 3]
The answer would be:
[0 0 0.3333]

Why is Matlab doing this ? matlab, matrix MATLAB Answers — New Questions

​

lambda = [1.064e-6];
R = [30];
w=[0.001];
q = (1./R- i* lambda./pi./w.^2).^(-1);
a=1;
p=1;
m=1;
probe_r=linspace(0,0.003,100);
probe_theta=linspace(0,0.003,100);
rseed=[0*max(w):max(w)/30:3*max(w)];
thetaseed=[0:360]*pi/180;
[r,theta]=meshgrid(rseed,thetaseed);
E=LaguerreGaussianE([p,m,q,lambda,a],r,theta);
V=interp2(r,theta,E,probe_r,probe_theta);

column_names = {‘r’, ‘theta’, ‘V’};
fid = fopen(‘fidtext.txt’,’wt’);
fprintf(fid, ‘%s ‘, column_names{:});
fprintf(fid, ‘n’);
block_of_data = [probe_r, probe_theta, V];
fmt = repmat(‘%15g ‘, 1, 3);
fmt(end:end+1) = ‘n’;
fprintf(fid, fmt, block_of_data.’);
fclose(fid);

With the current code I have I get a .txt file of only the real numbers from my function V along with the values of probe_r and probe_theta. How do I alter this to produce both the real and complex numbers as a 3 column .txt file of r, theta and V as I am unable to see a formatSpec to include complex numbers.lambda = [1.064e-6];
R = [30];
w=[0.001];
q = (1./R- i* lambda./pi./w.^2).^(-1);
a=1;
p=1;
m=1;
probe_r=linspace(0,0.003,100);
probe_theta=linspace(0,0.003,100);
rseed=[0*max(w):max(w)/30:3*max(w)];
thetaseed=[0:360]*pi/180;
[r,theta]=meshgrid(rseed,thetaseed);
E=LaguerreGaussianE([p,m,q,lambda,a],r,theta);
V=interp2(r,theta,E,probe_r,probe_theta);

column_names = {‘r’, ‘theta’, ‘V’};
fid = fopen(‘fidtext.txt’,’wt’);
fprintf(fid, ‘%s ‘, column_names{:});
fprintf(fid, ‘n’);
block_of_data = [probe_r, probe_theta, V];
fmt = repmat(‘%15g ‘, 1, 3);
fmt(end:end+1) = ‘n’;
fprintf(fid, fmt, block_of_data.’);
fclose(fid);

With the current code I have I get a .txt file of only the real numbers from my function V along with the values of probe_r and probe_theta. How do I alter this to produce both the real and complex numbers as a 3 column .txt file of r, theta and V as I am unable to see a formatSpec to include complex numbers. lambda = [1.064e-6];
R = [30];
w=[0.001];
q = (1./R- i* lambda./pi./w.^2).^(-1);
a=1;
p=1;
m=1;
probe_r=linspace(0,0.003,100);
probe_theta=linspace(0,0.003,100);
rseed=[0*max(w):max(w)/30:3*max(w)];
thetaseed=[0:360]*pi/180;
[r,theta]=meshgrid(rseed,thetaseed);
E=LaguerreGaussianE([p,m,q,lambda,a],r,theta);
V=interp2(r,theta,E,probe_r,probe_theta);

column_names = {‘r’, ‘theta’, ‘V’};
fid = fopen(‘fidtext.txt’,’wt’);
fprintf(fid, ‘%s ‘, column_names{:});
fprintf(fid, ‘n’);
block_of_data = [probe_r, probe_theta, V];
fmt = repmat(‘%15g ‘, 1, 3);
fmt(end:end+1) = ‘n’;
fprintf(fid, fmt, block_of_data.’);
fclose(fid);

With the current code I have I get a .txt file of only the real numbers from my function V along with the values of probe_r and probe_theta. How do I alter this to produce both the real and complex numbers as a 3 column .txt file of r, theta and V as I am unable to see a formatSpec to include complex numbers. fprintf, complex numbers MATLAB Answers — New Questions

​

I got the following error when I tried a simulink model (named UPFC compensated line) of a power system using only the Simscape powersystem specialised technology components.
Complex signal mismatch. Input port 1 of ‘UPFC_compensated_line/Unified Power Flow Controller (Phasor Type)/Power Components Modeling (Shunt & Series Converter)/Power Components modeling (Shunt Converter)/Demux1’ expects a signal of numeric type real. However, it is driven by a signal of numeric type complex
Component: Simulink | Category: Block error
Complex signal mismatch. Output port 1 of ‘UPFC_compensated_line/Unified Power Flow Controller (Phasor Type)/Power Components Modeling (Shunt & Series Converter)/Power Components modeling (Shunt Converter)/dq-axis model of 3-ph RL branch’ is a signal of numeric type complex. However, it is driving a signal of numeric type real

What might be the reason? Please help to resolve the problemI got the following error when I tried a simulink model (named UPFC compensated line) of a power system using only the Simscape powersystem specialised technology components.
Complex signal mismatch. Input port 1 of ‘UPFC_compensated_line/Unified Power Flow Controller (Phasor Type)/Power Components Modeling (Shunt & Series Converter)/Power Components modeling (Shunt Converter)/Demux1’ expects a signal of numeric type real. However, it is driven by a signal of numeric type complex
Component: Simulink | Category: Block error
Complex signal mismatch. Output port 1 of ‘UPFC_compensated_line/Unified Power Flow Controller (Phasor Type)/Power Components Modeling (Shunt & Series Converter)/Power Components modeling (Shunt Converter)/dq-axis model of 3-ph RL branch’ is a signal of numeric type complex. However, it is driving a signal of numeric type real

What might be the reason? Please help to resolve the problem I got the following error when I tried a simulink model (named UPFC compensated line) of a power system using only the Simscape powersystem specialised technology components.
Complex signal mismatch. Input port 1 of ‘UPFC_compensated_line/Unified Power Flow Controller (Phasor Type)/Power Components Modeling (Shunt & Series Converter)/Power Components modeling (Shunt Converter)/Demux1’ expects a signal of numeric type real. However, it is driven by a signal of numeric type complex
Component: Simulink | Category: Block error
Complex signal mismatch. Output port 1 of ‘UPFC_compensated_line/Unified Power Flow Controller (Phasor Type)/Power Components Modeling (Shunt & Series Converter)/Power Components modeling (Shunt Converter)/dq-axis model of 3-ph RL branch’ is a signal of numeric type complex. However, it is driving a signal of numeric type real

What might be the reason? Please help to resolve the problem upfc, signal mismatch MATLAB Answers — New Questions

​

Hello,
I’m using the envelope function with ‘peak’ option, [yupper,ylower] = envelope(x,np,’peak’).
The description says, ‘The envelopes are determined using spline interpolation over local maxima separated by at least np samples.’
I’m wondering that what order is used in the spline interpolation and also how the local maxima are calculated, separated by ‘at least’ np samples.
If I can, I would like to know the algorithm too.
Thank you.Hello,
I’m using the envelope function with ‘peak’ option, [yupper,ylower] = envelope(x,np,’peak’).
The description says, ‘The envelopes are determined using spline interpolation over local maxima separated by at least np samples.’
I’m wondering that what order is used in the spline interpolation and also how the local maxima are calculated, separated by ‘at least’ np samples.
If I can, I would like to know the algorithm too.
Thank you. Hello,
I’m using the envelope function with ‘peak’ option, [yupper,ylower] = envelope(x,np,’peak’).
The description says, ‘The envelopes are determined using spline interpolation over local maxima separated by at least np samples.’
I’m wondering that what order is used in the spline interpolation and also how the local maxima are calculated, separated by ‘at least’ np samples.
If I can, I would like to know the algorithm too.
Thank you. envelope, algorithm MATLAB Answers — New Questions

​

In a 3D comsol model, how to output the solution on a specific line (in reference or current configuration)? In .mph comsol models, it’s like the
comp1.at3
command that is equivalent.In a 3D comsol model, how to output the solution on a specific line (in reference or current configuration)? In .mph comsol models, it’s like the
comp1.at3
command that is equivalent. In a 3D comsol model, how to output the solution on a specific line (in reference or current configuration)? In .mph comsol models, it’s like the
comp1.at3
command that is equivalent. comsol MATLAB Answers — New Questions

​

I have simulated FOC with BLDC using mathematical model but Iam getting negative speed. Can anyone help me identify what is the issue with the mode. I have used the same with bldc motor block(PMSM block with trapizoidal back emf) that is working but this one is going to negative speed.I have simulated FOC with BLDC using mathematical model but Iam getting negative speed. Can anyone help me identify what is the issue with the mode. I have used the same with bldc motor block(PMSM block with trapizoidal back emf) that is working but this one is going to negative speed. I have simulated FOC with BLDC using mathematical model but Iam getting negative speed. Can anyone help me identify what is the issue with the mode. I have used the same with bldc motor block(PMSM block with trapizoidal back emf) that is working but this one is going to negative speed. foc MATLAB Answers — New Questions

​

I would like to do image processing to recognize a needle in an ultrasound image. Such as the one below
This image is scaled down to 20%

Currently I am thinking about blurring, created a region of interest, thresholding the image and then identifiying the most continous blob of white pixels.
Does anyone have a way to create a blur that utilizes the brightest pixels to change the pixels around it?
Other feedback and methods and input would also be appreciated.
Also @Image Analyst might be interetstedI would like to do image processing to recognize a needle in an ultrasound image. Such as the one below
This image is scaled down to 20%

Currently I am thinking about blurring, created a region of interest, thresholding the image and then identifiying the most continous blob of white pixels.
Does anyone have a way to create a blur that utilizes the brightest pixels to change the pixels around it?
Other feedback and methods and input would also be appreciated.
Also @Image Analyst might be interetsted I would like to do image processing to recognize a needle in an ultrasound image. Such as the one below
This image is scaled down to 20%

Currently I am thinking about blurring, created a region of interest, thresholding the image and then identifiying the most continous blob of white pixels.
Does anyone have a way to create a blur that utilizes the brightest pixels to change the pixels around it?
Other feedback and methods and input would also be appreciated.
Also @Image Analyst might be interetsted image analysis, digital image processing, ultrasound MATLAB Answers — New Questions

​