Hi all, I am trying to install matlab on my ubuntu 18.04. I have matlab license from the institution and do the sudo installation accordingly.

The problem is that whenever I get to the download stage, it stuck/freeze at less than 2 percent and pop-up error telling me the download has to restart.

So apparently googling "matlab download freeze" tells me that this is a common occurrence on windows due to security settings/firewall etc.

So I turned off my ubuntu UFW, but the issue persisted even after rebooting.

May I know anyone has similar experience and how did you overcome this? Thanks!Hi all, I am trying to install matlab on my ubuntu 18.04. I have matlab license from the institution and do the sudo installation accordingly.

The problem is that whenever I get to the download stage, it stuck/freeze at less than 2 percent and pop-up error telling me the download has to restart.

So apparently googling "matlab download freeze" tells me that this is a common occurrence on windows due to security settings/firewall etc.

So I turned off my ubuntu UFW, but the issue persisted even after rebooting.

May I know anyone has similar experience and how did you overcome this? Thanks! Hi all, I am trying to install matlab on my ubuntu 18.04. I have matlab license from the institution and do the sudo installation accordingly.

The problem is that whenever I get to the download stage, it stuck/freeze at less than 2 percent and pop-up error telling me the download has to restart.

So apparently googling "matlab download freeze" tells me that this is a common occurrence on windows due to security settings/firewall etc.

So I turned off my ubuntu UFW, but the issue persisted even after rebooting.

May I know anyone has similar experience and how did you overcome this? Thanks! ubuntu, download, installation MATLAB Answers — New Questions

​

s= [0,5*10^-2,1*10^-1,1.5*10^-1,2*10^-1,2.5*10^-1,3*10^-1,3.5*10^-1,4*10^-1,4.5*10^-1,5*10^-1,5.5*10^-1,6*10^-1,6.5*10^-1,7*10^-1,7.5*10^-1,8*10^-1,8.5*10^-1,9*10^-1,9.5*10^-1,9.9*10^-1,9.99*10^-1,1,1,1,1,1];
omega= [3.74*10^-1, 3.8*10^-1,3.87*10^-1,3.94*10^-1,4.02*10^-1,4.11*10^-1,4.2*10^-1,4.29*10^-1,4.4*10^-1,4.51*10^-1,4.64*10^-1,4.78*10^-1,4.94*10^-1,5.12*10^-1,5.33*10^-1,5.57*10^-1,5.86*10^-1,6.23*10^-1,6.72*10^-1,7.46*10^-1,8.71*10^-1,9.56*10^-1,9.68*10^-1,9.72*10^-1,9.75*10^-1,9.8*10^-1,9.86*10^-1];
Z1= 1+((1-(s.^2)).^(1/3)).*(((1+s).^(1/3))+(1-s).^(1/3));
Z2= sqrt(3*(s.^2)+(Z1.^2));
Ri= 3+Z2-sqrt((3-Z1).*(3+Z1+2.*Z2));R= 3+Z2+(sqrt((3-Z1).*(3+Z1+2.*Z2)));
Wi= (s+((Ri).^(3/2))).^-1;
WR= (s+((R).^(3/2))).^-1;
Rlr= 2*(1+cos((2/3)*acos(-s)));
Wlr= (s+((Rlr).^2)).^-1;
plot(s,Wi,’r’);
hold on
plot(s,WR,’y’),plot(s,Wlr,’b’),plot(s,omega,’g’);
hold off
xlabel(‘spin’),ylabel(‘angular velocity’),legend(‘r-prograde Isco angular velocity’,’y-retrograde ISCO angular velocity’,’b-LR angular velocity’,’g-given angular velocity’);
y= diff(Wi);
s=diff(s);
c=diff(omega);
plot(s,y,’r’);
how can I plot this ? because it doesn’t give me to plot y over ss= [0,5*10^-2,1*10^-1,1.5*10^-1,2*10^-1,2.5*10^-1,3*10^-1,3.5*10^-1,4*10^-1,4.5*10^-1,5*10^-1,5.5*10^-1,6*10^-1,6.5*10^-1,7*10^-1,7.5*10^-1,8*10^-1,8.5*10^-1,9*10^-1,9.5*10^-1,9.9*10^-1,9.99*10^-1,1,1,1,1,1];
omega= [3.74*10^-1, 3.8*10^-1,3.87*10^-1,3.94*10^-1,4.02*10^-1,4.11*10^-1,4.2*10^-1,4.29*10^-1,4.4*10^-1,4.51*10^-1,4.64*10^-1,4.78*10^-1,4.94*10^-1,5.12*10^-1,5.33*10^-1,5.57*10^-1,5.86*10^-1,6.23*10^-1,6.72*10^-1,7.46*10^-1,8.71*10^-1,9.56*10^-1,9.68*10^-1,9.72*10^-1,9.75*10^-1,9.8*10^-1,9.86*10^-1];
Z1= 1+((1-(s.^2)).^(1/3)).*(((1+s).^(1/3))+(1-s).^(1/3));
Z2= sqrt(3*(s.^2)+(Z1.^2));
Ri= 3+Z2-sqrt((3-Z1).*(3+Z1+2.*Z2));R= 3+Z2+(sqrt((3-Z1).*(3+Z1+2.*Z2)));
Wi= (s+((Ri).^(3/2))).^-1;
WR= (s+((R).^(3/2))).^-1;
Rlr= 2*(1+cos((2/3)*acos(-s)));
Wlr= (s+((Rlr).^2)).^-1;
plot(s,Wi,’r’);
hold on
plot(s,WR,’y’),plot(s,Wlr,’b’),plot(s,omega,’g’);
hold off
xlabel(‘spin’),ylabel(‘angular velocity’),legend(‘r-prograde Isco angular velocity’,’y-retrograde ISCO angular velocity’,’b-LR angular velocity’,’g-given angular velocity’);
y= diff(Wi);
s=diff(s);
c=diff(omega);
plot(s,y,’r’);
how can I plot this ? because it doesn’t give me to plot y over s s= [0,5*10^-2,1*10^-1,1.5*10^-1,2*10^-1,2.5*10^-1,3*10^-1,3.5*10^-1,4*10^-1,4.5*10^-1,5*10^-1,5.5*10^-1,6*10^-1,6.5*10^-1,7*10^-1,7.5*10^-1,8*10^-1,8.5*10^-1,9*10^-1,9.5*10^-1,9.9*10^-1,9.99*10^-1,1,1,1,1,1];
omega= [3.74*10^-1, 3.8*10^-1,3.87*10^-1,3.94*10^-1,4.02*10^-1,4.11*10^-1,4.2*10^-1,4.29*10^-1,4.4*10^-1,4.51*10^-1,4.64*10^-1,4.78*10^-1,4.94*10^-1,5.12*10^-1,5.33*10^-1,5.57*10^-1,5.86*10^-1,6.23*10^-1,6.72*10^-1,7.46*10^-1,8.71*10^-1,9.56*10^-1,9.68*10^-1,9.72*10^-1,9.75*10^-1,9.8*10^-1,9.86*10^-1];
Z1= 1+((1-(s.^2)).^(1/3)).*(((1+s).^(1/3))+(1-s).^(1/3));
Z2= sqrt(3*(s.^2)+(Z1.^2));
Ri= 3+Z2-sqrt((3-Z1).*(3+Z1+2.*Z2));R= 3+Z2+(sqrt((3-Z1).*(3+Z1+2.*Z2)));
Wi= (s+((Ri).^(3/2))).^-1;
WR= (s+((R).^(3/2))).^-1;
Rlr= 2*(1+cos((2/3)*acos(-s)));
Wlr= (s+((Rlr).^2)).^-1;
plot(s,Wi,’r’);
hold on
plot(s,WR,’y’),plot(s,Wlr,’b’),plot(s,omega,’g’);
hold off
xlabel(‘spin’),ylabel(‘angular velocity’),legend(‘r-prograde Isco angular velocity’,’y-retrograde ISCO angular velocity’,’b-LR angular velocity’,’g-given angular velocity’);
y= diff(Wi);
s=diff(s);
c=diff(omega);
plot(s,y,’r’);
how can I plot this ? because it doesn’t give me to plot y over s derivative, diff MATLAB Answers — New Questions

​

There is a machine, in my institution, with an installed and activated Matlab (R2019b). The installation/activation was likely done using an online account. Is it possible to figure out the account associated to that license? ThanksThere is a machine, in my institution, with an installed and activated Matlab (R2019b). The installation/activation was likely done using an online account. Is it possible to figure out the account associated to that license? Thanks There is a machine, in my institution, with an installed and activated Matlab (R2019b). The installation/activation was likely done using an online account. Is it possible to figure out the account associated to that license? Thanks licensing, online account, activation MATLAB Answers — New Questions

​

I am stimulating the branch line coupler in the first ppt you have given, and I found that while it could stimulate S11 and S31 in different wavelength, I found it could not stimulate S21 and S41 correctly, instead they overlap with S31 and S21 respectively.
This is my code:
Z=50;
Z1=Z/sqrt(2);
omega=0:0.01:2;
syms S11 S21 S31 S41;
for i=1:1:201
resOdd=1;
Odd1=[1,0;-1i/Z,1];
Odd2=[cos(pi/2/omega(i)),1i*Z1*sin(pi/2/omega(i));1i*sin(pi/2/omega(i))/Z1,cos(pi/2/omega(i))];
resOdd=resOdd*Odd1*Odd2*Odd1;
AOdd=resOdd(1,1);
BOdd=resOdd(1,2);
COdd=resOdd(2,1);
DOdd=resOdd(2,2);
S11Odd=(-DOdd+AOdd-Z*COdd+BOdd/Z)/(DOdd+AOdd+Z*COdd+BOdd/Z);
S21Odd=2*Z/(Z*DOdd+Z*AOdd+Z^2*COdd+BOdd);
resEven=1;
Even1=[1,0;1i/Z,1];
Even2=[cos(pi/2/omega(i)),1i*Z1*sin(pi/2/omega(i));1i*sin(pi/2/omega(i))/Z1,cos(pi/2/omega(i))];
resEven=resEven*Even1*Even2*Even1;
AEven=resEven(1,1);
BEven=resEven(1,2);
CEven=resEven(2,1);
DEven=resEven(2,2);
S11Even=(-DEven+AEven-Z*CEven+BEven/Z)/(DEven+AEven+Z*CEven+BEven/Z);
S21Even=2*Z/(Z*DEven+Z*AEven+Z^2*CEven+BEven);
S11a=abs((S11Even+S11Odd)/2);
S21a=abs((S21Even+S21Odd)/2);
S31a=abs((S21Even-S21Odd)/2);
S41a=abs((S11Even-S11Odd)/2);
S11(i)=20*log10(S11a);
S21(i)=20*log10(S21a);
S31(i)=20*log10(S31a);
S41(i)=20*log10(S41a);
end
plot(omega,S11)
hold on
plot(omega,S21)
hold on
plot(omega,S31)
hold on
plot(omega,S41)
hold off
P1 is my stimulation result:
and P2&P3 is the stimulation result from other research of branch line coupler:

I have checked both my code and the scale conversion between frequency and wavelength, from which I still could not understand how the results are not matched.I am stimulating the branch line coupler in the first ppt you have given, and I found that while it could stimulate S11 and S31 in different wavelength, I found it could not stimulate S21 and S41 correctly, instead they overlap with S31 and S21 respectively.
This is my code:
Z=50;
Z1=Z/sqrt(2);
omega=0:0.01:2;
syms S11 S21 S31 S41;
for i=1:1:201
resOdd=1;
Odd1=[1,0;-1i/Z,1];
Odd2=[cos(pi/2/omega(i)),1i*Z1*sin(pi/2/omega(i));1i*sin(pi/2/omega(i))/Z1,cos(pi/2/omega(i))];
resOdd=resOdd*Odd1*Odd2*Odd1;
AOdd=resOdd(1,1);
BOdd=resOdd(1,2);
COdd=resOdd(2,1);
DOdd=resOdd(2,2);
S11Odd=(-DOdd+AOdd-Z*COdd+BOdd/Z)/(DOdd+AOdd+Z*COdd+BOdd/Z);
S21Odd=2*Z/(Z*DOdd+Z*AOdd+Z^2*COdd+BOdd);
resEven=1;
Even1=[1,0;1i/Z,1];
Even2=[cos(pi/2/omega(i)),1i*Z1*sin(pi/2/omega(i));1i*sin(pi/2/omega(i))/Z1,cos(pi/2/omega(i))];
resEven=resEven*Even1*Even2*Even1;
AEven=resEven(1,1);
BEven=resEven(1,2);
CEven=resEven(2,1);
DEven=resEven(2,2);
S11Even=(-DEven+AEven-Z*CEven+BEven/Z)/(DEven+AEven+Z*CEven+BEven/Z);
S21Even=2*Z/(Z*DEven+Z*AEven+Z^2*CEven+BEven);
S11a=abs((S11Even+S11Odd)/2);
S21a=abs((S21Even+S21Odd)/2);
S31a=abs((S21Even-S21Odd)/2);
S41a=abs((S11Even-S11Odd)/2);
S11(i)=20*log10(S11a);
S21(i)=20*log10(S21a);
S31(i)=20*log10(S31a);
S41(i)=20*log10(S41a);
end
plot(omega,S11)
hold on
plot(omega,S21)
hold on
plot(omega,S31)
hold on
plot(omega,S41)
hold off
P1 is my stimulation result:
and P2&P3 is the stimulation result from other research of branch line coupler:

I have checked both my code and the scale conversion between frequency and wavelength, from which I still could not understand how the results are not matched. I am stimulating the branch line coupler in the first ppt you have given, and I found that while it could stimulate S11 and S31 in different wavelength, I found it could not stimulate S21 and S41 correctly, instead they overlap with S31 and S21 respectively.
This is my code:
Z=50;
Z1=Z/sqrt(2);
omega=0:0.01:2;
syms S11 S21 S31 S41;
for i=1:1:201
resOdd=1;
Odd1=[1,0;-1i/Z,1];
Odd2=[cos(pi/2/omega(i)),1i*Z1*sin(pi/2/omega(i));1i*sin(pi/2/omega(i))/Z1,cos(pi/2/omega(i))];
resOdd=resOdd*Odd1*Odd2*Odd1;
AOdd=resOdd(1,1);
BOdd=resOdd(1,2);
COdd=resOdd(2,1);
DOdd=resOdd(2,2);
S11Odd=(-DOdd+AOdd-Z*COdd+BOdd/Z)/(DOdd+AOdd+Z*COdd+BOdd/Z);
S21Odd=2*Z/(Z*DOdd+Z*AOdd+Z^2*COdd+BOdd);
resEven=1;
Even1=[1,0;1i/Z,1];
Even2=[cos(pi/2/omega(i)),1i*Z1*sin(pi/2/omega(i));1i*sin(pi/2/omega(i))/Z1,cos(pi/2/omega(i))];
resEven=resEven*Even1*Even2*Even1;
AEven=resEven(1,1);
BEven=resEven(1,2);
CEven=resEven(2,1);
DEven=resEven(2,2);
S11Even=(-DEven+AEven-Z*CEven+BEven/Z)/(DEven+AEven+Z*CEven+BEven/Z);
S21Even=2*Z/(Z*DEven+Z*AEven+Z^2*CEven+BEven);
S11a=abs((S11Even+S11Odd)/2);
S21a=abs((S21Even+S21Odd)/2);
S31a=abs((S21Even-S21Odd)/2);
S41a=abs((S11Even-S11Odd)/2);
S11(i)=20*log10(S11a);
S21(i)=20*log10(S21a);
S31(i)=20*log10(S31a);
S41(i)=20*log10(S41a);
end
plot(omega,S11)
hold on
plot(omega,S21)
hold on
plot(omega,S31)
hold on
plot(omega,S41)
hold off
P1 is my stimulation result:
and P2&P3 is the stimulation result from other research of branch line coupler:

I have checked both my code and the scale conversion between frequency and wavelength, from which I still could not understand how the results are not matched. branchline coupler, s-parameter MATLAB Answers — New Questions

​

This is not a question per se but a comment/complaint. Changing current directory is so disruptive. All I want to do is to just read the example. Anyone feel the same ?This is not a question per se but a comment/complaint. Changing current directory is so disruptive. All I want to do is to just read the example. Anyone feel the same ? This is not a question per se but a comment/complaint. Changing current directory is so disruptive. All I want to do is to just read the example. Anyone feel the same ? matlab help command window example MATLAB Answers — New Questions

​

Hi, please help me to devewlop the code of the following problem.
I am stating the problem in very simple way:
”I have 15 by 15 matrix. I need a closed path with shortest distance value starting from node i and passing through every node exactly once, where i= 1,2,…15.”
distanceMatrix = [
0 1214 1450 1130 1669 1964 1474 2286 827 176 1377 1735 1465 302 1471;
1214 0 2143 114 1802 1808 2169 1241 1158 1293 1074 690 2168 1029 2164;
1450 2143 0 2126 1317 1568 139 3318 989 1266 1679 2767 71.2 1360 30.8;
1130 114 2126 0 1852 1859 2196 1192 1185 1226 1125 641 2194 985 2190;
1669 1802 1317 1852 0 400 1451 3027 852 1591 754 2476 1383 1378 1317;
1964 1808 1568 1859 400 0 1700 3034 1159 1883 742 2483 1632 1670 1566;
1474 2169 139 2196 1451 1700 0 3349 1072 1289 1799 2798 69.4 1383 142;
2286 1241 3318 1192 3027 3034 3349 0 2341 2377 2300 560 3370 2135 3365;
827 1158 989 1185 852 1159 1072 2341 0 746 781 1792 1048 544 1016;
176 1293 1266 1226 1591 1883 1289 2377 746 0 1319 1825 1283 276 1290;
1377 1074 1679 1125 754 742 1799 2300 781 1319 0 1715 1695 1087 1656;
1735 690 2767 641 2476 2483 2798 560 1792 1825 1715 0 2819 1584 2815;
1465 2168 71.2 2194 1383 1632 69.4 3370 1048 1283 1695 2819 0 1389 87.4;
302 1029 1360 985 1378 1670 1383 2135 544 276 1087 1584 1389 0 1385;
1471 2164 30.8 2190 1317 1566 142 3365 1016 1290 1656 2815 87.4 1385 0
];Hi, please help me to devewlop the code of the following problem.
I am stating the problem in very simple way:
”I have 15 by 15 matrix. I need a closed path with shortest distance value starting from node i and passing through every node exactly once, where i= 1,2,…15.”
distanceMatrix = [
0 1214 1450 1130 1669 1964 1474 2286 827 176 1377 1735 1465 302 1471;
1214 0 2143 114 1802 1808 2169 1241 1158 1293 1074 690 2168 1029 2164;
1450 2143 0 2126 1317 1568 139 3318 989 1266 1679 2767 71.2 1360 30.8;
1130 114 2126 0 1852 1859 2196 1192 1185 1226 1125 641 2194 985 2190;
1669 1802 1317 1852 0 400 1451 3027 852 1591 754 2476 1383 1378 1317;
1964 1808 1568 1859 400 0 1700 3034 1159 1883 742 2483 1632 1670 1566;
1474 2169 139 2196 1451 1700 0 3349 1072 1289 1799 2798 69.4 1383 142;
2286 1241 3318 1192 3027 3034 3349 0 2341 2377 2300 560 3370 2135 3365;
827 1158 989 1185 852 1159 1072 2341 0 746 781 1792 1048 544 1016;
176 1293 1266 1226 1591 1883 1289 2377 746 0 1319 1825 1283 276 1290;
1377 1074 1679 1125 754 742 1799 2300 781 1319 0 1715 1695 1087 1656;
1735 690 2767 641 2476 2483 2798 560 1792 1825 1715 0 2819 1584 2815;
1465 2168 71.2 2194 1383 1632 69.4 3370 1048 1283 1695 2819 0 1389 87.4;
302 1029 1360 985 1378 1670 1383 2135 544 276 1087 1584 1389 0 1385;
1471 2164 30.8 2190 1317 1566 142 3365 1016 1290 1656 2815 87.4 1385 0
]; Hi, please help me to devewlop the code of the following problem.
I am stating the problem in very simple way:
”I have 15 by 15 matrix. I need a closed path with shortest distance value starting from node i and passing through every node exactly once, where i= 1,2,…15.”
distanceMatrix = [
0 1214 1450 1130 1669 1964 1474 2286 827 176 1377 1735 1465 302 1471;
1214 0 2143 114 1802 1808 2169 1241 1158 1293 1074 690 2168 1029 2164;
1450 2143 0 2126 1317 1568 139 3318 989 1266 1679 2767 71.2 1360 30.8;
1130 114 2126 0 1852 1859 2196 1192 1185 1226 1125 641 2194 985 2190;
1669 1802 1317 1852 0 400 1451 3027 852 1591 754 2476 1383 1378 1317;
1964 1808 1568 1859 400 0 1700 3034 1159 1883 742 2483 1632 1670 1566;
1474 2169 139 2196 1451 1700 0 3349 1072 1289 1799 2798 69.4 1383 142;
2286 1241 3318 1192 3027 3034 3349 0 2341 2377 2300 560 3370 2135 3365;
827 1158 989 1185 852 1159 1072 2341 0 746 781 1792 1048 544 1016;
176 1293 1266 1226 1591 1883 1289 2377 746 0 1319 1825 1283 276 1290;
1377 1074 1679 1125 754 742 1799 2300 781 1319 0 1715 1695 1087 1656;
1735 690 2767 641 2476 2483 2798 560 1792 1825 1715 0 2819 1584 2815;
1465 2168 71.2 2194 1383 1632 69.4 3370 1048 1283 1695 2819 0 1389 87.4;
302 1029 1360 985 1378 1670 1383 2135 544 276 1087 1584 1389 0 1385;
1471 2164 30.8 2190 1317 1566 142 3365 1016 1290 1656 2815 87.4 1385 0
]; matlab MATLAB Answers — New Questions

​

Hello all, I am currently coding a value/policy function iteration code to solve an economic problem. The problem that this is solving for is that there is a certain amount of S and R currently owned (these are the state variables). S and R generate income, which must be allocated to purchase more S and R for the next "round" (S’ and R’, equivalently. S’ and R’ are the control variables). This code solves for the optimal amount of income that should be spent on R’ at any given level of S and R. Right now my code below only works at specific parameter values and does not converge–the error term gets stuck, and I cannot seem to figure out why. My code is below. Any help on what could possibly prevent it from converging would be greatly appreciated. Thank you all in advance!
clc;
clear;
% Model Parameters and initalize matrices
numRisky = 100;
numSafe = 100;
v = (1:numSafe)’ + (1:numRisky);
vnew=ones(numSafe, numRisky);
policy = zeros(numSafe, numRisky);
delta=0.99;
ctol=0.001;
ntol=200;
count=0;
norm=1;
vx = ones(numSafe);
vy = ones(numRisky);
alpha = 1;
while norm>ctol && count<ntol
for R=1:numRisky
riskyReturns = 1:1:3;
numReturns = numel(riskyReturns);
pReturns = (1/numReturns)*ones(1, numReturns);
u = zeros(1, numReturns);
p = zeros(1, numReturns);
for S = 1:numSafe
safeReturn = 2;
for r = 1:numReturns
safeIncome = S*safeReturn;
riskyIncome = R*riskyReturns(r);
posI = safeIncome + riskyIncome;
U = alpha*log(S+R);
W = zeros(1, posI);

for i = 1:posI
if S + (posI-i) > numSafe && R + i <= numRisky
vBar = v(numSafe,R+i) + vx(R+i)*(S+posI-i-numSafe); %conduct linear approximation
W(i) = (1-delta)*U + delta*vBar;
elseif R + i > numRisky && S + (posI-i) <= numSafe
vBar = v(S+posI-i,numRisky) + vy(S+posI-i)*(R+i-numRisky);
W(i) = (1-delta)*U + delta*vBar; %conduct linear approximation
elseif R + i > numRisky && S + (posI-i) > numSafe
vBar = v(numSafe,numRisky) + vx(numRisky)*(S+posI-i-numSafe) + vy(numSafe)*(R+i-numRisky);
W(i) = (1-delta)*U + delta*vBar; %conduct linear approximation
else
W(i) = (1-delta)*U + delta*v(S + (posI – i), R + i); %no linear approximation needed here
end
end
[w, ind] = max(W);
u(r) = w;
p(r) = ind/posI;
end
vnew(S, R) = dot(u, pReturns);
policy(S, R) = dot(p, pReturns);
end
end
%calculate slopes for linear approximation if S’ and R’ are greater than numSafe and numRisky, respectively
for n = 1:numSafe
vx(n) = v(numSafe,n)-v(numSafe-1,n);
end
for m = 1:numRisky
vy(m) = v(m,numRisky)-v(m,numRisky-1);
end
%calculate convergence criteria and update value function
count = count + 1;
norm = max(abs(log(vnew(:)) – log(v(:))));
v = vnew;
disp(‘Iteration:’)
disp(count)
disp(‘Error Term:’)
disp(norm)
endHello all, I am currently coding a value/policy function iteration code to solve an economic problem. The problem that this is solving for is that there is a certain amount of S and R currently owned (these are the state variables). S and R generate income, which must be allocated to purchase more S and R for the next "round" (S’ and R’, equivalently. S’ and R’ are the control variables). This code solves for the optimal amount of income that should be spent on R’ at any given level of S and R. Right now my code below only works at specific parameter values and does not converge–the error term gets stuck, and I cannot seem to figure out why. My code is below. Any help on what could possibly prevent it from converging would be greatly appreciated. Thank you all in advance!
clc;
clear;
% Model Parameters and initalize matrices
numRisky = 100;
numSafe = 100;
v = (1:numSafe)’ + (1:numRisky);
vnew=ones(numSafe, numRisky);
policy = zeros(numSafe, numRisky);
delta=0.99;
ctol=0.001;
ntol=200;
count=0;
norm=1;
vx = ones(numSafe);
vy = ones(numRisky);
alpha = 1;
while norm>ctol && count<ntol
for R=1:numRisky
riskyReturns = 1:1:3;
numReturns = numel(riskyReturns);
pReturns = (1/numReturns)*ones(1, numReturns);
u = zeros(1, numReturns);
p = zeros(1, numReturns);
for S = 1:numSafe
safeReturn = 2;
for r = 1:numReturns
safeIncome = S*safeReturn;
riskyIncome = R*riskyReturns(r);
posI = safeIncome + riskyIncome;
U = alpha*log(S+R);
W = zeros(1, posI);

for i = 1:posI
if S + (posI-i) > numSafe && R + i <= numRisky
vBar = v(numSafe,R+i) + vx(R+i)*(S+posI-i-numSafe); %conduct linear approximation
W(i) = (1-delta)*U + delta*vBar;
elseif R + i > numRisky && S + (posI-i) <= numSafe
vBar = v(S+posI-i,numRisky) + vy(S+posI-i)*(R+i-numRisky);
W(i) = (1-delta)*U + delta*vBar; %conduct linear approximation
elseif R + i > numRisky && S + (posI-i) > numSafe
vBar = v(numSafe,numRisky) + vx(numRisky)*(S+posI-i-numSafe) + vy(numSafe)*(R+i-numRisky);
W(i) = (1-delta)*U + delta*vBar; %conduct linear approximation
else
W(i) = (1-delta)*U + delta*v(S + (posI – i), R + i); %no linear approximation needed here
end
end
[w, ind] = max(W);
u(r) = w;
p(r) = ind/posI;
end
vnew(S, R) = dot(u, pReturns);
policy(S, R) = dot(p, pReturns);
end
end
%calculate slopes for linear approximation if S’ and R’ are greater than numSafe and numRisky, respectively
for n = 1:numSafe
vx(n) = v(numSafe,n)-v(numSafe-1,n);
end
for m = 1:numRisky
vy(m) = v(m,numRisky)-v(m,numRisky-1);
end
%calculate convergence criteria and update value function
count = count + 1;
norm = max(abs(log(vnew(:)) – log(v(:))));
v = vnew;
disp(‘Iteration:’)
disp(count)
disp(‘Error Term:’)
disp(norm)
end Hello all, I am currently coding a value/policy function iteration code to solve an economic problem. The problem that this is solving for is that there is a certain amount of S and R currently owned (these are the state variables). S and R generate income, which must be allocated to purchase more S and R for the next "round" (S’ and R’, equivalently. S’ and R’ are the control variables). This code solves for the optimal amount of income that should be spent on R’ at any given level of S and R. Right now my code below only works at specific parameter values and does not converge–the error term gets stuck, and I cannot seem to figure out why. My code is below. Any help on what could possibly prevent it from converging would be greatly appreciated. Thank you all in advance!
clc;
clear;
% Model Parameters and initalize matrices
numRisky = 100;
numSafe = 100;
v = (1:numSafe)’ + (1:numRisky);
vnew=ones(numSafe, numRisky);
policy = zeros(numSafe, numRisky);
delta=0.99;
ctol=0.001;
ntol=200;
count=0;
norm=1;
vx = ones(numSafe);
vy = ones(numRisky);
alpha = 1;
while norm>ctol && count<ntol
for R=1:numRisky
riskyReturns = 1:1:3;
numReturns = numel(riskyReturns);
pReturns = (1/numReturns)*ones(1, numReturns);
u = zeros(1, numReturns);
p = zeros(1, numReturns);
for S = 1:numSafe
safeReturn = 2;
for r = 1:numReturns
safeIncome = S*safeReturn;
riskyIncome = R*riskyReturns(r);
posI = safeIncome + riskyIncome;
U = alpha*log(S+R);
W = zeros(1, posI);

for i = 1:posI
if S + (posI-i) > numSafe && R + i <= numRisky
vBar = v(numSafe,R+i) + vx(R+i)*(S+posI-i-numSafe); %conduct linear approximation
W(i) = (1-delta)*U + delta*vBar;
elseif R + i > numRisky && S + (posI-i) <= numSafe
vBar = v(S+posI-i,numRisky) + vy(S+posI-i)*(R+i-numRisky);
W(i) = (1-delta)*U + delta*vBar; %conduct linear approximation
elseif R + i > numRisky && S + (posI-i) > numSafe
vBar = v(numSafe,numRisky) + vx(numRisky)*(S+posI-i-numSafe) + vy(numSafe)*(R+i-numRisky);
W(i) = (1-delta)*U + delta*vBar; %conduct linear approximation
else
W(i) = (1-delta)*U + delta*v(S + (posI – i), R + i); %no linear approximation needed here
end
end
[w, ind] = max(W);
u(r) = w;
p(r) = ind/posI;
end
vnew(S, R) = dot(u, pReturns);
policy(S, R) = dot(p, pReturns);
end
end
%calculate slopes for linear approximation if S’ and R’ are greater than numSafe and numRisky, respectively
for n = 1:numSafe
vx(n) = v(numSafe,n)-v(numSafe-1,n);
end
for m = 1:numRisky
vy(m) = v(m,numRisky)-v(m,numRisky-1);
end
%calculate convergence criteria and update value function
count = count + 1;
norm = max(abs(log(vnew(:)) – log(v(:))));
v = vnew;
disp(‘Iteration:’)
disp(count)
disp(‘Error Term:’)
disp(norm)
end value function iteration, bellman, dynamic program, matlab, bellman equation MATLAB Answers — New Questions

​

Buen dia
Estamos implementando la conexion de un motor nema 17 con driver L293 en Matlab pero en el momento de correr el codigo es lento
Se realiza la prueba en Arduino y funciona bien pero en matlab resulta lento
el motor se encuentra probado, al igual que el driver
se tiene la secuencia de funcionamiento de las bobinas Probada
pero hay un conflicto con la velocidad al momento de subirlo en matlab (Envia las Señales con un intervalo de tiempo muy lento aprox 1 segundo entre pasos)
le estaria muy agradecido si fuera posible responder la duda.
Adjunto el codigo:

function Paso_a_Paso(Pos_Sig)
global ARD; global BASE;
if Pos_Sig>BASE
DIR=0;
Pasos=((Pos_Sig-BASE)*8);
end
if Pos_Sig<BASE
DIR=1;
Pasos=((BASE-Pos_Sig)*8);
end
steps={‘1010’, ‘0010’, ‘0110’, ‘0100’, ‘0101’, ‘0001’, ‘1001’ ,’1000′}; %SECUENCIA

if DIR==1
for(p=0:1:Pasos)
for (i = 1:1:8)
step = steps{i}; % step=1,2,3,4
writeDigitalPin(ARD,’D8′ ,str2double(step(1)));
writeDigitalPin(ARD,’D9′ ,str2double(step(2)));
writeDigitalPin(ARD,’D10′ , str2double(step(3)));
writeDigitalPin(ARD,’D11′ ,str2double(step(4)));
end
end
end
if DIR==0
for(p=0:1:Pasos)
for (i = 8:1:1)
step = steps{i}; % step=1,2,3,4
writeDigitalPin(ARD,’D8′ ,str2double(step(1)));
writeDigitalPin(ARD,’D9′ ,str2double(step(2)));
writeDigitalPin(ARD,’D10′ , str2double(step(3)));
writeDigitalPin(ARD,’D11′ ,str2double(step(4)));
end
end
end
BASE=Pos_Sig;
writeDigitalPin(ARD,’D8′ ,0);
writeDigitalPin(ARD,’D9′ ,0);
writeDigitalPin(ARD,’D10′ , 0);
writeDigitalPin(ARD,’D11′ ,0);
end

Gracias y buen diaBuen dia
Estamos implementando la conexion de un motor nema 17 con driver L293 en Matlab pero en el momento de correr el codigo es lento
Se realiza la prueba en Arduino y funciona bien pero en matlab resulta lento
el motor se encuentra probado, al igual que el driver
se tiene la secuencia de funcionamiento de las bobinas Probada
pero hay un conflicto con la velocidad al momento de subirlo en matlab (Envia las Señales con un intervalo de tiempo muy lento aprox 1 segundo entre pasos)
le estaria muy agradecido si fuera posible responder la duda.
Adjunto el codigo:

function Paso_a_Paso(Pos_Sig)
global ARD; global BASE;
if Pos_Sig>BASE
DIR=0;
Pasos=((Pos_Sig-BASE)*8);
end
if Pos_Sig<BASE
DIR=1;
Pasos=((BASE-Pos_Sig)*8);
end
steps={‘1010’, ‘0010’, ‘0110’, ‘0100’, ‘0101’, ‘0001’, ‘1001’ ,’1000′}; %SECUENCIA

if DIR==1
for(p=0:1:Pasos)
for (i = 1:1:8)
step = steps{i}; % step=1,2,3,4
writeDigitalPin(ARD,’D8′ ,str2double(step(1)));
writeDigitalPin(ARD,’D9′ ,str2double(step(2)));
writeDigitalPin(ARD,’D10′ , str2double(step(3)));
writeDigitalPin(ARD,’D11′ ,str2double(step(4)));
end
end
end
if DIR==0
for(p=0:1:Pasos)
for (i = 8:1:1)
step = steps{i}; % step=1,2,3,4
writeDigitalPin(ARD,’D8′ ,str2double(step(1)));
writeDigitalPin(ARD,’D9′ ,str2double(step(2)));
writeDigitalPin(ARD,’D10′ , str2double(step(3)));
writeDigitalPin(ARD,’D11′ ,str2double(step(4)));
end
end
end
BASE=Pos_Sig;
writeDigitalPin(ARD,’D8′ ,0);
writeDigitalPin(ARD,’D9′ ,0);
writeDigitalPin(ARD,’D10′ , 0);
writeDigitalPin(ARD,’D11′ ,0);
end

Gracias y buen dia Buen dia
Estamos implementando la conexion de un motor nema 17 con driver L293 en Matlab pero en el momento de correr el codigo es lento
Se realiza la prueba en Arduino y funciona bien pero en matlab resulta lento
el motor se encuentra probado, al igual que el driver
se tiene la secuencia de funcionamiento de las bobinas Probada
pero hay un conflicto con la velocidad al momento de subirlo en matlab (Envia las Señales con un intervalo de tiempo muy lento aprox 1 segundo entre pasos)
le estaria muy agradecido si fuera posible responder la duda.
Adjunto el codigo:

function Paso_a_Paso(Pos_Sig)
global ARD; global BASE;
if Pos_Sig>BASE
DIR=0;
Pasos=((Pos_Sig-BASE)*8);
end
if Pos_Sig<BASE
DIR=1;
Pasos=((BASE-Pos_Sig)*8);
end
steps={‘1010’, ‘0010’, ‘0110’, ‘0100’, ‘0101’, ‘0001’, ‘1001’ ,’1000′}; %SECUENCIA

if DIR==1
for(p=0:1:Pasos)
for (i = 1:1:8)
step = steps{i}; % step=1,2,3,4
writeDigitalPin(ARD,’D8′ ,str2double(step(1)));
writeDigitalPin(ARD,’D9′ ,str2double(step(2)));
writeDigitalPin(ARD,’D10′ , str2double(step(3)));
writeDigitalPin(ARD,’D11′ ,str2double(step(4)));
end
end
end
if DIR==0
for(p=0:1:Pasos)
for (i = 8:1:1)
step = steps{i}; % step=1,2,3,4
writeDigitalPin(ARD,’D8′ ,str2double(step(1)));
writeDigitalPin(ARD,’D9′ ,str2double(step(2)));
writeDigitalPin(ARD,’D10′ , str2double(step(3)));
writeDigitalPin(ARD,’D11′ ,str2double(step(4)));
end
end
end
BASE=Pos_Sig;
writeDigitalPin(ARD,’D8′ ,0);
writeDigitalPin(ARD,’D9′ ,0);
writeDigitalPin(ARD,’D10′ , 0);
writeDigitalPin(ARD,’D11′ ,0);
end

Gracias y buen dia ok MATLAB Answers — New Questions

​

I would have expected an error if the second argument is strictly required. Giving me a wrong map just pushes the error to runtime when I call it with an ‘double’ key. No idea what went into this design.I would have expected an error if the second argument is strictly required. Giving me a wrong map just pushes the error to runtime when I call it with an ‘double’ key. No idea what went into this design. I would have expected an error if the second argument is strictly required. Giving me a wrong map just pushes the error to runtime when I call it with an ‘double’ key. No idea what went into this design. containers map MATLAB Answers — New Questions

​

I am encountering several issues with deepNetworkDesigner.
(1) I have a large network whose layers do not fit entirely on one screen, or at least not at a magnification level that allows me to easily see, distinguish, and manipulate the layers. How do I connect one layer to a second layer that is off-screen? When I try to drag a connection wire toward the destination layer, the pane does not autoscroll to let me reach it.
(2) I have one layer A that fans out to many output layers C1,C2,…,CN. I now wish to replace A with a different layer. Is there a way to do that, without having to tediously and manually reconnect the N outputs? Similarly, what if I wish to insert a layer B, so that A now feeds only to B but it is B that fans out to C1,C2,C3,…CN. Again, can this be done without breaking all N connections and having to manually redraw them?
(3) Is there a way to undo a change to a network (something like Ctrl-Z)?
(4) Is there a way to reorient the pane so the network flows left to right, instead of top to bottom (something like a landscape mode, rather than portrait).I am encountering several issues with deepNetworkDesigner.
(1) I have a large network whose layers do not fit entirely on one screen, or at least not at a magnification level that allows me to easily see, distinguish, and manipulate the layers. How do I connect one layer to a second layer that is off-screen? When I try to drag a connection wire toward the destination layer, the pane does not autoscroll to let me reach it.
(2) I have one layer A that fans out to many output layers C1,C2,…,CN. I now wish to replace A with a different layer. Is there a way to do that, without having to tediously and manually reconnect the N outputs? Similarly, what if I wish to insert a layer B, so that A now feeds only to B but it is B that fans out to C1,C2,C3,…CN. Again, can this be done without breaking all N connections and having to manually redraw them?
(3) Is there a way to undo a change to a network (something like Ctrl-Z)?
(4) Is there a way to reorient the pane so the network flows left to right, instead of top to bottom (something like a landscape mode, rather than portrait). I am encountering several issues with deepNetworkDesigner.
(1) I have a large network whose layers do not fit entirely on one screen, or at least not at a magnification level that allows me to easily see, distinguish, and manipulate the layers. How do I connect one layer to a second layer that is off-screen? When I try to drag a connection wire toward the destination layer, the pane does not autoscroll to let me reach it.
(2) I have one layer A that fans out to many output layers C1,C2,…,CN. I now wish to replace A with a different layer. Is there a way to do that, without having to tediously and manually reconnect the N outputs? Similarly, what if I wish to insert a layer B, so that A now feeds only to B but it is B that fans out to C1,C2,C3,…CN. Again, can this be done without breaking all N connections and having to manually redraw them?
(3) Is there a way to undo a change to a network (something like Ctrl-Z)?
(4) Is there a way to reorient the pane so the network flows left to right, instead of top to bottom (something like a landscape mode, rather than portrait). deepnetworkdesigner MATLAB Answers — New Questions

​