LOFTER for ipad —— 让兴趣,更有趣

点击下载 关闭

LOFTER-网易轻博

MATLAB

5066浏览    481参与
zqhh

坐标轴

set(gcf,'unit','centimeters','position',[1,2,20,15])

%其中gcf是控制绘图区位置大小,[1,2]表示绘图区距离显示器左下角1,2厘米,[20,15]表示绘图区大小为宽20,高15厘米。

set(gca,'positon',[0.1,0.1,0.9,0.9])

%其中gca是控制坐标,[0.1,0.1,0.9,0.9]是图形相对绘图区的比例,[0.1,0.1]是图形相对绘图区左下角的位置,[0.9,0.9]是图形相对绘图区的大小比例。

grid on

%可以在画的图像中添加栅格,用命令grid on,这样可以方便你对齐某条线或是对...

set(gcf,'unit','centimeters','position',[1,2,20,15])

%其中gcf是控制绘图区位置大小,[1,2]表示绘图区距离显示器左下角1,2厘米,[20,15]表示绘图区大小为宽20,高15厘米。

set(gca,'positon',[0.1,0.1,0.9,0.9])

%其中gca是控制坐标,[0.1,0.1,0.9,0.9]是图形相对绘图区的比例,[0.1,0.1]是图形相对绘图区左下角的位置,[0.9,0.9]是图形相对绘图区的大小比例。

grid on

%可以在画的图像中添加栅格,用命令grid on,这样可以方便你对齐某条线或是对比比较方便,

text(1.5,0.3,'cos(x)')    

%将cosx这个注解加到坐标中的某个位置

CarryOn
CarryOn
CarryOn
研数科技
時の砂漠

【转】--- 解决.m文件默认打开方式不能设置为MATLAB的方法(MATLAB 2019b)

最近更新matlab2019b后发现无论用什么方法

比如 右键-打开方式-添加matlab.exe等方式,都无法将matlab.exe添加到可以打开的程序中。

所以.m文件只能通过先打开matlab,然后再在matlab资源管理器中找到相应的.m文件,

然后才能打开。

这样太麻烦了。没想到网络上已经有大神给出了解决办法。

这里分享一下,以做备份。


原文地址:https://blog.csdn.net/yaningli/article/details/84538798


这里作为该文章的转载

内容:

1. 将下面代码复制保存为 associateFiles...

最近更新matlab2019b后发现无论用什么方法

比如 右键-打开方式-添加matlab.exe等方式,都无法将matlab.exe添加到可以打开的程序中。

所以.m文件只能通过先打开matlab,然后再在matlab资源管理器中找到相应的.m文件,

然后才能打开。

这样太麻烦了。没想到网络上已经有大神给出了解决办法。

这里分享一下,以做备份。


原文地址:https://blog.csdn.net/yaningli/article/details/84538798


这里作为该文章的转载

内容:

1. 将下面代码复制保存为 associateFiles.m 文件。

2. 在MATLAB中运行 associateFiles.m 文件,即在命令行中输入 associateFiles 然后 回车。

3. 执行步骤2后,程序运行会生成文件 MatlabFileAssocFix.reg

4. 关闭MATLAB,双击运行 MatlabFileAssocFix.reg 。

5. 设置.m文件默认打开程序为MATLAB2019b 完工。



代码:

  1. function associateFiles(action, userExtList, fileStr)

  2. % associateFiles(action, extList, fileStr)

  3. %

  4. % Makes a registry files that can be used to set correct file associantions on

  5. % a windows platform. The following MATLAB file extensions are supported:

  6. % .m, .mat, .fig, .mexw32, .mexw64, .p, .mdl, .mdlp, .slx, .mldatx, .req,

  7. % .sldd, .slddc, .slxp, .sltx, .mn, .mu, .muphlp, .xvc, .xvz, .ssc, .mlapp,

  8. % .mlappinstall, .mltbx, .mlpkginstall, .mlprj

  9. %

  10. % INPUT:

  11. % action  - optional string. 

  12. %           * 'add' (default) adds/rewrites the MATLAB file association registry

  13. %              keysfor this version.

  14. %           * 'delete' deletes the MATLAB file association registry entries for

  15. %              ALL versions of MATLAB (including "old style" ones)

  16. %           * 'deleteadd' is the same as 'delete' followed by 'add'

  17. % extList - optional string or cell array of strings containing the file

  18. %           extensions that should be associated with this version. Default is

  19. %           all MATLAB file extension (see above).

  20. % fileStr - optional string with the name of the registry file to be written

  21. %           (possibly including path). Default is the file

  22. %           'MatlabFileAssocFix.reg' in the current directory.

  23. %

  24. % USAGE:

  25. % 1) Run with desired options (see above). A registry file should have been

  26. %    created. 

  27. % 2) Exit all running instances of MATLAB.

  28. % 3) Make a backup copy of the windows registry if you need to restore the

  29. %    changes, see https://support.microsoft.com/en-us/kb/322756

  30. % 4) Double click on the created file (possibly need to enter a password) and

  31. %    confirm.

  32. % 5) Restart Windows (or explorer.exe).

  33. % 6) The MATLAB files should now be associated with the MATLAB version that the

  34. %    registry file was created in and e.g. m-files should be opened in an

  35. %    already running instance of MATLAB.

  36. %

  37. % EXAMPLES:

  38. % * associateFiles('deleteadd') - Makes a registry files that deletes all

  39. %   previous MATLAB file association registry keysandwrite new ones that

  40. %   associates all MATLAB files with the MATLAB version that the registry file

  41. %   was created in.

  42. % * associateFiles('', {'.m', '.mat', '.fig'}, 'myFile') - Makes a registry file

  43. %   "myFile.reg" that associates m-, mat- and fig-files with the MATLAB version

  44. %   that the registry file was created in. 

  45. %

  46. % VERSION 1.0

  47. % Defualt input

  48. if (nargin < 1 || isempty(action))

  49.   action      = 'add';

  50. end

  51. if (nargin < 2)

  52.   userExtList = {};

  53. end

  54. if (nargin < 3)

  55.   fileStr = '';

  56. end

  57. if (~iscell(userExtList))

  58. if (isempty(userExtList))

  59.     userExtList = {};

  60. else

  61.     userExtList = {userExtList};

  62.   end

  63. end

  64. % Sanity check

  65. if (~ischar(action) || (~strcmpi(action, 'add') && ...

  66.     ~strcmpi(action, 'delete') && ~strcmpi(action, 'deleteadd')))

  67.   error('The action to perform must be ''add'', ''delete'' or ''deleteadd''!')

  68. end

  69. if (~isempty(userExtList) && ~min(cellfun(@ischar, userExtList)))

  70.   error('The file extension list must be a string or a cell array of strings!')

  71. end

  72. if (~ischar(fileStr))

  73.   error('The file to write to must be a string!')

  74. end

  75. % Get the currently running MATLAB version

  76. verStr = regexp(version, '(\d*?\.\d*?\.\d*?)\.', 'tokens');

  77. verStr = verStr{1}{1};

  78. verNum = str2double(regexprep(verStr, '(\d*?\.\d*)[\x0000-\xffff]*', '$1'));

  79. verHex = sprintf('x', str2double(regexprep(verStr, ...

  80. '(\d*?)\.[\x0000-\xffff]*', '$1')), str2double(regexprep(verStr, ...

  81. '\d*?\.(\d*?)\.[\x0000-\xffff]*', '$1')));

  82. % Get 32/64-bit

  83. arch = computer;

  84. switch arch

  85.   case 'PCWIN'

  86.     binFolder = 'win32';

  87.   case 'PCWIN64'

  88.     binFolder = 'win64';

  89. end

  90. binPath = fullfile(matlabroot, 'bin', binFolder);

  91. % Known MATLAB files with possible DDE actions

  92. fileExtCell = {...

  93. 'fig' ,   'MATLAB Figure'              , '-62'                       , ...

  94.   {'Open', 'uiopen(''%1'',1)'}           , []                          ; ...

  95. 'm'     , 'MATLAB Code'                , '-58'                       , ...

  96.   {'Open', 'uiopen(''%1'',1)'}           , {'Run', 'run(''%1'')'}      ; ...

  97. 'mat'   , 'MATLAB Data'                , '-59'                       , ...

  98.   {'Load', 'load(''%1'')'    }           , {'Open', 'uiimport(''%1'')'}; ...

  99. 'mdl'   , 'Simulink Model'             , '-61'                       , ...

  100.   {'Load', 'uiopen(''%1'',1)'}           , []                          ; ...

  101. 'mdlp'  , 'Simulink Protected Model'   , '-72'                       , ...

  102.   []                                     , []                          ; ...

  103. 'mexw32', 'MATLAB MEX'                 , '-63'                       , ...

  104.   []                                     , []                          ; ...

  105. 'mexw64', 'MATLAB MEX'                 , '-63'                       , ...

  106.   []                                     , []                          ; ...

  107. 'mn'    , 'MuPAD Notebook'             , '-66'                       , ...

  108.   {'Open', 'mupad(''%1'')'}              , []                          ; ...

  109. 'mu'    , 'MuPAD Code'                 , '-67'                       , ...

  110.   {'Open', 'uiopen(''%1'',1)'}           , []                          ; ...

  111. 'muphlp', 'MuPAD Help'                 , '-68'                       , ...

  112.   {'Open', 'doc(symengine, ''%1'')'}     , []                          ; ...

  113. 'p'     , 'MATLAB P-code'              , '-60'                       , ...

  114.   []                                     , []                          ; ...

  115. 'slx'   , 'Simulink Model (SLX format)', '-73'                       , ...

  116.   {'Open', 'uiopen(''%1'',1)'}           , []                          ; ...

  117. 'ssc'   , 'Simscape Model'             , '-65'                       , ...

  118.   {'Open', 'uiopen(''%1'',1)'}           , []                          ; ...

  119. 'xvc'   , 'MuPAD Graphics'             , '-69'                       , ...

  120.   {'Open', 'mupad(''%1'')'}              , []                          ; ...

  121. 'xvz'   , 'MuPAD Graphics'             , '-70'                       , ...

  122.   {'Open', 'mupad(''%1'')'}              , []                          ; ...

  123. 'mlapp'       , 'MATLAB Application'              , [] , [], []      ; ... 

  124. 'mltbx'       , 'MATLAB Toolbox'                  , [] , [], []      ; ... 

  125. 'mldatx'      , 'Simulink Scenario'               , [] , [], []      ; ...  

  126. 'req'         , 'Simulink Requirements Link'      , [] , [], []      ; ... 

  127. 'sldd'        , 'Simulink Dictionary'             , [] , [], []      ; ... 

  128. 'slddc'       , 'Simulink Dictionary'             , [] , [], []      ; ...      

  129. 'mlappinstall', 'MATLAB Application'              , [] , [], []      ; ...  

  130. 'mlpkginstall', 'MATLAB Support Package'          , [] , [], []      ; ... 

  131. 'slxp'        , 'Simulink Protected Model Package', [] , [], []      ; ... 

  132. 'sltx'        , 'Simulink Template'               , [] , [], []      ; ... 

  133. 'mlprj'       , 'MATLAB Project'                  , [] , [], []};

  134. % Possibly trim list

  135. if (~isempty(userExtList))

  136.   fileExtCell = fileExtCell(ismember(fileExtCell(:, 1), ...

  137.     regexprep(userExtList, '\.', '')), :);

  138. end

  139. % Make registry file

  140. if (~isempty(fileStr))

  141.   % Possibly add file extension

  142.   [~, ~, tmp] = fileparts(fileStr);

  143. if (isempty(tmp))

  144.     fileStr = [fileStr, '.reg'];

  145.   end

  146.   fid = fopen(fileStr, 'w');

  147. else

  148.   fid = fopen('MatlabFileAssocFix.reg', 'w');

  149. end

  150. if (fid == -1)

  151.   error('Failed to create registry file')

  152. end

  153. % Write intial lines

  154. fprintf(fid, '%s\r\n\r\n', 'Windows Registry Editor Version 5.00');

  155. fprintf(fid, '%s\r\n\r\n', ';FIXES MATLAB FILE ASSOCIATIONS');

  156. % REMOVE OLD KEYS

  157. explorerKey = ['HKEY_CURRENT_USER\Software\Microsoft\Windows\', ...

  158.   'CurrentVersion\Explorer\FileExts'];

  159. % Iterate over file extensions

  160. for fileExtNo = 1 : size(fileExtCell, 1)

  161.   rmKeys  = {};

  162.   fileExt = fileExtCell{fileExtNo, 1};

  163.   % File extension keys

  164.   [status, result] = dos(['reg query HKEY_CLASSES_ROOT /f .', fileExt, ...

  165.     ' /k /e']);

  166.   if (~status)

  167.     keys = regexp(result, '(HKEY_CLASSES_ROOT[\x0000-\xffff]*?)\n', 'tokens');

  168.     rmKeys = [rmKeys, keys{:}];

  169.   end

  170.   % Old style keys without version numbers

  171.   if (~strcmpi(fileExt, 'mexw64'))

  172.     % Uses single DDE key for mex files

  173.     if (strcmpi(fileExt, 'mexw32'))

  174.       fileExtTmp = 'mex';

  175.     else

  176.       fileExtTmp = fileExt;

  177.     end

  178.     [status, result] = dos(['reg query HKEY_CLASSES_ROOT /f ', ...

  179.       fileExtTmp, 'file /k /e']);

  180.     if (~status)

  181.       keys = regexp(result, '(HKEY_CLASSES_ROOT[\x0000-\xffff]*?)\n', ...

  182.         'tokens');

  183.       rmKeys = [rmKeys, keys{:}];

  184.     end

  185.   end

  186.   % New style keys with version number

  187.   if (strcmpi(action, 'add'))

  188.     % Only remove keys related to this version

  189.     [status, result] = dos(['reg query HKEY_CLASSES_ROOT /f MATLAB.', ...

  190.       fileExt, '.', verStr ' /k']);

  191.   else

  192.     % Remove keys related to ALL version

  193.     [status, result] = dos(['reg query HKEY_CLASSES_ROOT /f MATLAB.', ...

  194.       fileExt, '. /k']);

  195.   end

  196.   if (~status)

  197.     keys = regexp(result, '(HKEY_CLASSES_ROOT[\x0000-\xffff]*?)\n', 'tokens');

  198.     rmKeys = [rmKeys, keys{:}];

  199.   end

  200.   % Explorer keys

  201.   [status, result] = dos(['reg query ', explorerKey, ' /f .', fileExt, ...

  202.     ' /k /e']);

  203.   if (~status)

  204.     keys = regexp(result, '(HKEY_CURRENT_USER[\x0000-\xffff]*?)\n', 'tokens');

  205.     rmKeys = [rmKeys, keys{:}];

  206.   end

  207.   % Write to file

  208.   if (~isempty(rmKeys))

  209.     fprintf(fid, '%s\r\n\r\n', [';REMOVES ', upper(fileExt), ...

  210.       ' FILE ASSOCIATIONS']);

  211.     for keyNo = 1 : length(rmKeys)

  212.       key = rmKeys{keyNo};

  213.       fprintf(fid, '%s\r\n\r\n', ['[-', key, ']']);

  214.     end

  215.   end

  216. end

  217. % ADD KEYS

  218. if (~strcmpi(action, 'delete'))

  219.   % Get text Persistent Handler

  220.   [status, result] = dos(...

  221.     'reg query HKEY_CLASSES_ROOT\.txt\PersistentHandler /ve');

  222.   if (~status)

  223.     PersistentHandler = regexp(result, '\{[\x0000-\xffff]*?\}', 'match');

  224.     PersistentHandler = PersistentHandler{1};

  225.   else

  226.     PersistentHandler = '';

  227.   end

  228.   % DDE call

  229.   ddeCall = 'ShellVerbs.Matlab';

  230.   if (verNum > 8)

  231.     % Changed from R2013a

  232.     ddeCall = [ddeCall, '.', verStr];

  233.   end

  234.   % Default icon

  235.   defIcon = 'm';

  236.   if (~exist(fullfile(binPath, 'm.ico'), 'file'))

  237.     defIcon = '';

  238.   end

  239.   % Path to MATLAB binary directory with \\

  240.   binPathStr = regexprep(binPath, '\\', '\\\\');

  241.   % Write Shell Open key

  242.   key = ['[HKEY_CLASSES_ROOT\Applications\MATLAB.exe\shell\open', ...

  243.     '\command]%r', '@="\"', binPathStr, '\\MATLAB.exe\" \"%1\""%r%r'];

  244.   fprintf(fid, '%s\r\n\r\n', ';ADD SHELL OPEN');

  245.   lines = regexp(key, '([\x0000-\xffff]*?)%r', 'tokens');

  246.   for lineNo = 1 : length(lines)

  247.     fprintf(fid, '%s\r\n', lines{lineNo}{1});

  248.   end

  249.   % Iterate over file types

  250.   for fileExtNo = 1 : size(fileExtCell, 1)

  251.     fileExt = fileExtCell{fileExtNo, 1};

  252.     % File extension keys

  253.     key  = ['[HKEY_CLASSES_ROOT\.', fileExt, ']%r@="MATLAB.', fileExt, '.', ...

  254.       verStr, '"%r'];

  255.     if (strcmpi(fileExt, 'm') && ~isempty(PersistentHandler))

  256.       % Add some values

  257.       key = [key, '"Content Type"="text/plain"%r', ...

  258.         '"PerceivedType"="Text"%r'];

  259.     end

  260.     key = [key, '%r'];

  261.     key = [key, '[HKEY_CLASSES_ROOT\.', fileExt, ...

  262.       '\OpenWithProgids]%r"MATLAB.', fileExt, '.', verStr, '"=""%r%r'];

  263.     if (strcmpi(fileExt, 'm') && ~isempty(PersistentHandler))

  264.       key = [key, '[HKEY_CLASSES_ROOT\.', fileExt, ...

  265.         '\PersistentHandler]%r@="', PersistentHandler, '"%r%r'];

  266.     end

  267.     key  = [key, '[HKEY_CLASSES_ROOT\.', fileExt, ...

  268.       '\Versions\MATLAB.', fileExt, '.' verStr, ']%r"FileVersionMS"=dword:', ...

  269.       verHex, '%r"FileVersionLS"=dword:00000000%r%r'];

  270.     % DDE keys

  271.     ddeData = fileExtCell(ismember(fileExtCell(:, 1), fileExt), :);

  272.     key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  273.       ']%r@="', ddeData{2}, '"%r'];

  274.     if (~isempty(ddeData{3}))

  275.       key = [key, '"FriendlyTypeName"="@', binPathStr, '\\matlab.exe', ...

  276.         ',', ddeData{3}, '"%r'];

  277.     end

  278.     key = [key, '%r'];

  279.     % Icon

  280.     icon = fileExt;

  281.     if (~exist(fullfile(binPath, [icon, '.ico']), 'file'))

  282.       icon = defIcon;

  283.     end

  284.     if (~isempty(icon))

  285.       key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  286.         '\DefaultIcon]%r@="', binPathStr, '\\', icon, '.ico,0"%r%r'];

  287.     end

  288.     % Shell actions

  289.     for shellActionNo = 4:5

  290.       ddePar = ddeData{shellActionNo};

  291.       if (~isempty(ddePar))

  292.         key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  293.           '\Shell\', ddePar{1}, ']%r@="', ddePar{1}, '"%r%r'];

  294.         key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  295.           '\Shell\', ddePar{1}, '\command]%r@="\"', binPathStr, ...

  296.           '\\matlab.exe\""%r%r'];

  297.         key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  298.           '\Shell\', ddePar{1}, '\ddeexec]%r@="', ddePar{2}, '"%r%r'];

  299.         key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  300.           '\Shell\', ddePar{1},'\ddeexec\application]%r@="', ...

  301.           ddeCall, '"%r%r'];

  302.         key = [key, '[HKEY_CLASSES_ROOT\MATLAB.', fileExt, '.' verStr, ...

  303.           '\Shell\', ddePar{1},'\ddeexec\topic]%r@="system"%r%r'];

  304.       end

  305.     end

  306.     % Explorer keys

  307.     key = [key, '[', explorerKey, '\.', fileExt, '\OpenWithProgids]%r'];

  308.     if (strcmpi(fileExt, 'm'))

  309.       key = [key, '"m_auto_file"=hex(0):%r'];

  310.     end

  311.     key = [key, '"MATLAB.', fileExt, '.',  verStr, '"=hex(0):%r%r'];

  312.     if (~isempty(ddeData{4}))

  313.       % Add key

  314.       key = [key, '[', explorerKey, '\.', fileExt, ...

  315.         '\OpenWithList]%r"a"="MATLAB.exe"%r"MRUList"="a"%r%r'];

  316.     else

  317.       key = [key, '[', explorerKey, '\.', fileExt, '\OpenWithList]%r%r'];

  318.     end

  319.     % Write to file

  320.     fprintf(fid, '%s\r\n\r\n', [';ADD ', upper(fileExt), ...

  321.       ' FILE ASSOCIATIONS']);

  322.     lines = regexp(key, '([\x0000-\xffff]*?)%r', 'tokens');

  323.     for lineNo = 1 : length(lines)

  324.       fprintf(fid, '%s\r\n', lines{lineNo}{1});

  325.     end

  326.   end

  327. end

  328. % Cloese file

  329. fclose(fid);


MATLAB2019b 亲测有用。



鹈鹕

imadd函数好好玩哦
(图三红月来自微博-不远归期,药郎是不知道从哪存的图,侵删)

imadd函数好好玩哦
(图三红月来自微博-不远归期,药郎是不知道从哪存的图,侵删)

护艳一生

Matlab-4.氢原子轨道轮廓图

采用Matlab绘制氢原子3dz2轨道轮廓图。

figure('numbertitle','off','name','Orbital profile of H atom')

x=-20:0.1:20;

z=-20:0.1:20;

[X,Z]=meshgrid(x,z);   %生成采样点矩阵

A=(2.*Z.^2-X.^2)./81./sqrt(6*pi)./exp(sqrt(X.^2+Z.^2)./3);

subplot(1,2,1);

[C,h]=contour(X,Z,A,6);

clabel(C,'manual');

xlabel('x'...

采用Matlab绘制氢原子3dz2轨道轮廓图。

figure('numbertitle','off','name','Orbital profile of H atom')

x=-20:0.1:20;

z=-20:0.1:20;

[X,Z]=meshgrid(x,z);   %生成采样点矩阵

A=(2.*Z.^2-X.^2)./81./sqrt(6*pi)./exp(sqrt(X.^2+Z.^2)./3);

subplot(1,2,1);

[C,h]=contour(X,Z,A,6);

clabel(C,'manual');

xlabel('x')

ylabel('z')

title('2D profile of 3d_{z^{2}} orbit');

grid on;

x=-25:0.5:25;

y=x;

z=x;

[X,Y,Z]=meshgrid(x,y,z);

V=(2.*Z.^2-X.^2-Y.^2)./81./sqrt(6)./exp(sqrt(X.^2+Y.^2+Z.^2)./3);

subplot(1,2,2)

patch(isosurface(X,Y,Z,V,0.005),'FaceColor','r','EdgeColor','none');  %提取V=0.005的等值面,并设置填充面颜色

hold on

patch(isosurface(X,Y,Z,V,-0.005),'FaceColor','c','EdgeColor','none');   %%提取V=-0.005的等值面,并设置填充面颜色

xlabel('x')

ylabel('y')

zlabel('z')

title('3D profile of 3d_{z^{2}} orbit')

view(3)

camlight

lighting gouraud

grid on



护艳一生

Matlab-3.绘制误差棒图。

       在[50:50:300]℃的温度条件下,分别对光催化剂效率进行三组相同实验,得到氮氧化物的转化率结果(%)分别为[11 35 50 71 77 82];[15 31 53 73 83 87];[17 27 44 66 80 87]。绘制转化率随温度变化的误差棒图。

A=[11 35 50 71 77 82;15 31 53 73 83 87;17 27 44 66 80 87];  %将转化率列入矩阵A中

Y=mean(A,1);  %对A中每一列求平均值

X=[50...

       在[50:50:300]℃的温度条件下,分别对光催化剂效率进行三组相同实验,得到氮氧化物的转化率结果(%)分别为[11 35 50 71 77 82];[15 31 53 73 83 87];[17 27 44 66 80 87]。绘制转化率随温度变化的误差棒图。

A=[11 35 50 71 77 82;15 31 53 73 83 87;17 27 44 66 80 87];  %将转化率列入矩阵A中

Y=mean(A,1);  %对A中每一列求平均值

X=[50 100 150 200 250 300];  %温度数据输入X向量

e=std(A,0,1);  %A矩阵列数据无偏估计标准差

plot(X,Y,'-r','linewidth',1.5);  %绘制转化率与温度折线图

hold on;

errorbar(X,Y,e,'k.','linewidth',1,'markeredgecolor','r','markersize',15);  %绘制误差棒图

axis([0,350,0,100]);  

set(gcf,'units','centimeters','position',[0.5,12,15,9.27]);

set(gca,'linewidth',1,'units','centimeters','position',[2,1.5,10,6.18],'Fontsize',12);

title('Variation of NO_{x} conversion with temperature','Fontsize',12);

xlabel('Temperature(\circC)','Fontsize',12);

ylabel('NO_{x} Conversion(%)','Fontsize',12);

set(gca,'ygrid','on');







护艳一生

Matlab-1.二维绘图与图形设置。

Matlab绘制y=sinx与y=e^x函数图像。

syms x;

f1=sin(x);

f2=exp(x);

fplot(f1,'r-');  %绘制y=sinx图形

hold on;  %保留已有图像

fplot(f2,'b-');  %绘制y=e^x图像

set(gcf,'Units','centimeters','position',[0.5,12,15,9.27]);  %设置图片位置与尺寸,gcf返回当前图形句柄

set(gca,'Units','centimeters','position...

Matlab绘制y=sinx与y=e^x函数图像。

syms x;

f1=sin(x);

f2=exp(x);

fplot(f1,'r-');  %绘制y=sinx图形

hold on;  %保留已有图像

fplot(f2,'b-');  %绘制y=e^x图像

set(gcf,'Units','centimeters','position',[0.5,12,15,9.27]);  %设置图片位置与尺寸,gcf返回当前图形句柄

set(gca,'Units','centimeters','position',[1.5,1.5,10,6.18]);  %设置坐标轴位置与图尺寸,gca返回当前图形坐标轴句柄

xlabel('x','Fontsize',12);  

ylabel('y','Fontsize',12);  %设置坐标轴名称与字体大小

title('sinx与e^x函数图像','Fontsize',12);  %设置图形标题与字体大小

set(gca,'linewidth',1);  %设置坐标轴线宽

set(get(gca,'Children'),'linewidth',1);  %设置图内线宽

axis([-4*pi,4*pi,-2,2]);  %设置坐标轴范围

text(2*pi/3,sin(2*pi/3),'\leftarrowy=sinx','HorizontalAlignment','left','Fontsize',10);  %图形注释y=sinx

text(0,exp(0),'y=e^x\rightarrow','HorizontalAlignment','right','Fontsize',10);  %图形注释y=e^x

set(gca,'ygrid','on');  %显示水平网格线

set(gca,'ytick',[-2:0.5:2]);  %设置网格线间距

legend('sinx','e^x');  %显示图例






hdw2000
Marsh

压缩机十系数仿真matlab代码

clear,clc

x1 = xlsread ('data.xls','B2:K2'); 

x2 = xlsread ('data.xls','A3:A7');

y = xlsread ('data.xls','B3:K7');

n1=length(x1);n2=length(x2);

x1=ones(n2,1)*x1;x1=x1(:);

x2=x2*ones(1,n1);x2=x2(:);

y=y(:); X=[x1,x2]; n=length(y);

str=num2str([1:n]');

fx1=@(c,x1,x2)(c(1)+c(2)*...

clear,clc

x1 = xlsread ('data.xls','B2:K2'); 

x2 = xlsread ('data.xls','A3:A7');

y = xlsread ('data.xls','B3:K7');

n1=length(x1);n2=length(x2);

x1=ones(n2,1)*x1;x1=x1(:);

x2=x2*ones(1,n1);x2=x2(:);

y=y(:); X=[x1,x2]; n=length(y);

str=num2str([1:n]');

fx1=@(c,x1,x2)(c(1)+c(2)*x1+c(3)*x2+c(4)*x1.^2+c(5)*x1.*x2+c(6)*x2.^2+c(7).*x1.^3+c(8)*x2.*x1.^2+c(9)*x1.*x2.^2+c(10)*x2.^3);

fx2=@(c,x1,x2)(c(1)+c(2)*X(:,1)+c(3)*X(:,2)+c(4)*X(:,1).^2+c(5)*X(:,1).*X(:,2)+c(6)*X(:,2).^2+c(7)*X(:,1).^3+c(8)*X(:,2).*X(:,1).^2+c(9)*X(:,1).*X(:,2).^2+c(10)*X(:,2).^3);

c=[8.29927747717989E+03 2.82383157394932E+02 -5.89576900114284E+01 4.10191749116755E+00 -1.50873655516070E+00 -4.11074321839960E-02 2.14355866374489E-02 -2.76297884478486E-02 -6.66250572642417E-03 -8.28569814231031E-04];

for l=1:5

    c=lsqcurvefit(fx2,c,X,y);

    c=nlinfit(X,y,fx2,c);

end

c

figure(1),clf

plot3(x1,x2,y,'o')

stem3(x1,x2,y,'filled')

text(x1,x2,y+.01,str)

hold on

[x11,x22]=meshgrid(min(x1):range(x1)/80:max(x1),min(x2):range(x2)/80:max(x2));

yhat=fx1(c,x11,x22);

surf(x11,x22,yhat)

shading interp

alpha(.8)

axis tight

xlabel('X1');ylabel('X2'),zlabel('Y')

SSy=var(y)*(n-1)

y1=fx1(c,x1,x2);

RSS=(y-y1)'*(y-y1)

rsquare=(SSy-RSS)/SSy

MSe=RSS/(n-length(c))

c=vpa(c,10)


连享会

连享会-空间计量专题现场班 (2019.6.27-30)


[toc]

1. 课程概览
  • 时间: 2019 年 6 月 27-30 (周四-周日)

  • 地点: 西安,西北工业大学国际会议中心 (百度地图 | 搜狗地图)

  • 主讲嘉宾:杨海生 (中山大学),第 1-4 讲;杨超 (上海财经大学),第 5-8 讲

  • 授课方式:

    • PPT + 电子板书

    • 每天 6 小时 (9:00-12:00;14:00-17:00)+半小时答疑

    • 讲义电子版于开课前一周发送;精美讲义由主办方统一印制

    • 软件:Matlab +...


[toc]

1. 课程概览
  • 时间: 2019 年 6 月 27-30 (周四-周日)

  • 地点: 西安,西北工业大学国际会议中心 (百度地图 | 搜狗地图)

  • 主讲嘉宾:杨海生 (中山大学),第 1-4 讲;杨超 (上海财经大学),第 5-8 讲

  • 授课方式:

    • PPT + 电子板书

    • 每天 6 小时 (9:00-12:00;14:00-17:00)+半小时答疑

    • 讲义电子版于开课前一周发送;精美讲义由主办方统一印制

    • 软件:Matlab + Stata

2. 主讲嘉宾简介

杨海生,中山大学岭南学院经济学系副教授,主要研究领域为空间计量经济学理论与应用、实证金融。在 Economic Geography、Ecological Economics、《经济研究》、《管理世界》、《管理科学学报》、《金融研究》等学术刊物上发表多篇论文,主持和参与多项国家自然科学基金、广东省自然科学基金等课题研究。

杨超,经济学博士,副教授,2015 年毕业于美国俄亥俄州立大学,现任教于上海财经大学经济学院。主讲“计量经济学”、“微观计量经济学”、“高等数理统计”、“高级金融经济学”等课程。主要研究方向为空间计量经济学、网络计量经济学、应用微观经济学。已在计量经济学顶级期刊 Journal of Econometrics、空间计量经济学领域权威期刊 Regional Science and Urban Economics 等期刊发表数篇学术论文,担任 Journal of Econometrics, Journal of Business & Economic Statistics 等国际期刊审稿人,获国家自然科学青年基金项目资助。

3. 课程详情3.1 课程介绍
  • 为什么要学空间计量?

    • 唇亡齿寒、城门失火殃及池鱼、近朱者赤近墨者黑

    • 空间相关、空间溢出效应如此普遍

  • 我们将学到什么?

    • 空间自回归模型 (SAR), 空间误差模型 (SEM)

    • 空间静态\动态面板模型

    • 空间交互影响

    • 引力模型

    • 同行效应

    • Probit\Tobit 空间计量模型

    • 20+ 篇 Top 期刊论文展示解读

可以想象,若德克萨斯州的失业率上升 10%,该州的谋杀率也必然会上升。然而,城门失火殃及池鱼,我们更关心的可能是:德州周边各州的谋杀率会如何变化呢?下面这幅 GIF 图形 便展示了基于空间自回归(SAR)模型估计出的空间溢出效应。

显然,在国际贸易、财政支出、区域创新、产业集聚、能源消费、环境污染、公司金融等诸多领域都存在这种空间相关或空间溢出效应。得益于空间计量经济学在过去几十年的快速发展,上述话题得以深入挖掘。

本课程将通过为期四天的讲解,帮助学员们了解并掌握空间相关性的产生缘由,各类空间计量模型的设定、估计和检验方法,并通过讲解发表于 Top 期刊的经典论文来展示这些模型的实际应用和实现过程。

具体来说,我们将回顾三代空间计量模型的由来、发展及演变历程。不仅会介绍各类空间模型的设定和估计方法,也会深入讨论模型参数和结果的经济含义。在此基础上,会进一步介绍面板空间计量模型、Probit/Tobit 空间计量模型,以及能够刻画空间交互影响、空间和社交网络关系的模型。

课程的关注重点在空间计量模型的应用,我们将利用经典文献的数据集对不同的空间计量模型进行演示,并提供分析结果所使用的 Matlab 和 Stata 程序,以便学员们能够使用这些程序来进行自己的实证分析。此外,我们还将介绍基本模型的拓展及其在金融、国际贸易、环境经济学、财政以及区域经济学等领域中的应用。

各个专题的具体介绍如下:

第一讲 首先回顾三代空间计量模型的设定思想和发展演变历程,以及不同空间模型之间的差异,以便各位能够合理设定和筛选模型;进而讲解空间计量模型中的三种主要估计方法:MLE、IV 和 GMM,配以 Matlab+Stata 编程实例。

第二讲 将通过一个热门话题 (同行效应) 的深入剖析来理解经典空间计量模型在金融领域特别是实证金融中的应用,以便各位了解从 “讲故事”到“建模型”这一过程中可能遇到的各种障碍和解决方法。首先,我们会通过对 Peer Effect 的刻画,把空间计量模型的应用扩展到实证金融的诸多领域,然后讲解 Peer Effect 中公司之间交互影响的微观理论基础。进而介绍 Peer Effect 的识别问题和参数估计的 Stata 实现。

第三讲 介绍高阶空间滞后模型。详细讲解有高阶空间滞后项的空间计量模型的空间权重设定、以及相应的 IV 及 GMM 估计原理和 Matlab + Stata 实现,最后简要介绍高阶空间滞后模型的拓展及发展方向。

第四讲 对前三讲介绍的空间计量的基本模型进行拓展,并结合金融、国际贸易、环境经济学、财政及区域经济学的经典文献帮助学员更好的理解和掌握空间计量模型的应用和实操。

第五讲 首先从国际贸易的具体例子 (Anderson, 1979; Anderson and Wincoop, 2004) 引出 「空间交互影响」 和 「引力模型」的应用背景,进而详述如何用矩阵刻画空间关系,最后讲解 「引力模型」 的 MLE 估计方法,并介绍 Matlab 编程方法和结果的解读。

第六讲 介绍空间和社交网络的计量模型与应用,包括:基本的网络结构,网络交互影响模型的微观基础、识别问题和 MLE 估计方法。我们会通过对空间关系的网络刻画,展现实际生活中各种不同的网络实例;进而讲解网络中个体之间产生交互影响的微观理论基础,以及网络中基本的(关于“同辈效应”)识别问题;最后讲解参数模型的估计方法。

第七讲 介绍空间面板数据模型的基本估计方法。对于静态空间面板模型,主要介绍包含固定效应的静态面板模型的直接估计法和间接估计法,以及如何通过 Matlab 编程实现;最后简要介绍动态空间面板模型的估计方法。

第八讲 介绍当被解释变量为二元选择或者截断数据时,空间计量模型的建立和估计,主要包括 Probit 空间计量模型和 Tobit 空间计量模型。这是近年来研究的热点和难点,在课程中首先介绍如何检验空间关联的存在性,再介绍模型的估计,并做实际应用演示。

3.2 课程大纲第 1 讲 何谓空间计量?(3小时)
  • 空间计量模型的由来、发展及演变历程

  • 几类主要的空间计量模型:SEM,SAR,SEM+SAR

  • 空间计量模型的估计方法:MLE、IV 和 GMM

  • 空间计量模型的 Matlab+Stata 实现

第 2 讲 实例剖析:Peer Effect (同群效应)(3小时)
  • Peer Effect 与空间计量:由来与设定

  • 资产定价中的 Peer Effect

  • 公司决策中的 Peer Effect

  • 范例论文:

    • Pirinsky C , Wang Q . Does Corporate Headquarters Location Matter for Stock Returns?[J]. Journal of Finance, 2006, 61(4):1991-2015. [PDF]

    • Leary M T, Roberts M R. Do Peer Firms Affect Corporate Financial Policy?. Journal of Finance, 2014, 69(1):139–178. [PDF]

    • Foucault, T., L. Fresard, Learning from peers' stock prices and corporate investment, Journal of Financial Economics, 2014, 111 (3): 554-577. [PDF]

第 3 讲 高阶空间滞后模型(3小时)
  • 高阶空间滞后模型的由来(Case et al., 1993)

  • 高阶空间滞后模型的估计(IV 和 GMM)

  • 高阶空间滞后模型的 Matlab+Stata 实现及应用

第 4 讲 空间计量的应用实例(3小时)
  • 空间计量在金融中的应用

  • 空间计量在财政、环境经济学的应用

  • 空间计量在国际贸易和区域经济学的应用

  • 范例论文:

    • Jillian Grennan . Dividend payments as a response to peer influence. Journal of Financial Economics, 2019, 131(3): 549-570. [PDF]

    • 白俊红, 蒋伏心. 协同创新、空间关联与区域创新绩效. 经济研究, 2015, 50(07): 174-187. [PDF]

    • 龙小宁, 朱艳丽, 蔡伟贤, 李少民. 基于空间计量模型的中国县级政府间税收竞争的实证分析. 经济研究, 2014(8): 41-53. [CNKI][PDF]

    • 尹恒, 徐琰超. 地市级地区问基本建设公共支出的相互影响. 经济研究, 2011(7): 55-64. [PDF]

    • Aklin M. Re-exploring the Trade and Environment Nexus Through the Diffusion of Pollution. Environmental & Resource Economics, 2016, 64(4): 663-682. [PDF]

第 5 讲 空间交互影响与引力模型(3小时)
  • 贸易成本与引力模型

  • 引力模型的矩阵表示

  • 基本引力模型的拓展

  • 引力模型的估计 – 范例论文

    • Anderson, James. A Theoretical Foundation for the Gravity Model of Factor Flows. No. 85. Boston College Department of Economics, 1978. [PDF]

    • Anderson, James E., and Eric Van Wincoop. "Trade costs." Journal of Economic literature 42.3 (2004): 691-751. [PDF]

    • LeSage, James, and Robert Kelley Pace. Introduction to spatial econometrics. Chapman and Hall/CRC, 2009. [Link]

第 6 讲 网络交互影响模型(3小时)
  • 空间关联与网络链接

  • 网络交互影响的微观基础

  • 同辈效应与识别 (the "Reflection Problem")

  • 基本交互影响模型的估计

  • 范例论文

    • Brock, William A., and Steven N. Durlauf. "Discrete choice with social interactions." Review of Economic Studies 68.2 (2001): 235-260. [PDF]

    • Manski, Charles F. "Identification of endogenous social effects: The reflection problem." Review of Economic Studies 60.3 (1993): 531-542. [PDF][Link]

    • Lee, Lung-fei, Ji Li, and Xu Lin. "Binary choice models with social network under heterogeneous rational expectations." Review of Economics and Statistics 96.3 (2014): 402-417. Link

    • Yang, Chao, and Lung-fei Lee. "Social interactions under incomplete information with heterogeneous expectations." Journal of Econometrics 198.1 (2017): 65-83. [PDF]

    • Yang, Chao, Lung-fei Lee, and Xi Qu. "Tobit models with social interactions: Complete vs incomplete information." Regional Science and Urban Economics 73 (2018): 30-50. ;[PDF]

    • Yang, Chao, and Lung-fei Lee. "Strategical interactions on municipal public safety spending with correlated private information." Regional Science and Urban Economics 72 (2018): 86–102. [PDF]

第 7 讲 空间面板数据模型(3小时)
  • 空间面板数据模型简介

  • 有固定效应的静态面板模型的估计

  • 动态面板模型

  • 范例论文

    • Lee, Lung-fei, and Jihai Yu. "Estimation of spatial autoregressive panel data models with fixed effects." Journal of Econometrics 154.2 (2010): 165-185. [PDF]

    • Lee, Lung-fei, and Jihai Yu. "Some recent developments in spatial panel data models." Regional Science and Urban Economics 40.5 (2010): 255-271. [PDF] [WP]

第 8 讲 离散和受限数据的空间计量模型(3小时)
  • Spatial Probit 与 Spatial Tobit 模型

  • Spatial Tobit 模型的检验和估计

  • 范例论文

    • Qu, Xi, and Lung-fei Lee. "LM tests for spatial correlation in spatial models with limited dependent variables." Regional Science and Urban Economics 42.3 (2012): 430-445. [PDF]

    • Qu, Xi, and Lung-fei Lee. "Locally most powerful tests for spatial interactions in the simultaneous SAR Tobit model." Regional Science and Urban Economics 43.2 (2013): 307-321. [PDF]

    • Xu, Xingbai, and Lung-fei Lee. "Maximum likelihood estimation of a spatial autoregressive Tobit model." Journal of Econometrics 188.1 (2015): 264-280. [PDF]

4. 报名信息
  • 主办方: 太原君泉教育咨询有限公司

  • 标准费用(含报名费、材料费),差旅及食宿费自理:

    • 全价:4200 元/人

    • 预报名价 (6月15日前报名) 4000元/人

    • 团报价:3600元/人(三人及以上)

    • 学生价:3700元/人(报到时需提供学生证原件)

  • 老学员优惠:

    • 参加过一次主办方培训的老学员:3500元/人

    • 参加过两次主办方培训的老学员:3200元/人

  • 友情高校感恩优惠:
    我们诚挚地感谢曾经团报的高校,来自以下高校的老师和学生可享受感恩价 3400 元/人:(按首字母拼音排序) 安徽财经大学、东北林业大学、河北大学、河北经贸大学、暨南大学、金陵学院、南开大学、山西财经大学、山西大学、太原科技大学、 运城学院、中南财经政法大学、中南大学、中山大学、中央财经大学。

  • Note: 以上各项优惠不能叠加使用。

  • 联系方式:

    • 邮箱:wjx004@sina.com

    • 电话 (微信同号):王老师 18903405450 ;李老师 ‭18636102467

    • 对公账户: 35117530000023891 (山西省太原市晋商银行南中环支行)

  • 温馨提示: 按报名顺序挑选\安排座位

长按/扫描二维码报名:

5. 交通和住宿
  • 校内宾馆: 西安西工大正禾宾馆

  • 酒店地址: 西安 莲湖区 友谊西路127号 [百度地图]

  • 住宿小知士: 正禾宾馆位于西北工业大学校内,上课、午休方便快捷,住宿环境可以参考以下图片:


七根毛线
一个负实特征根对应一个衰减,一...

一个负实特征根对应一个衰减,
一对负特征根对应一个衰减振荡,
因为全部特征根都在负实平面,所以正个系统是稳定的

一个负实特征根对应一个衰减,
一对负特征根对应一个衰减振荡,
因为全部特征根都在负实平面,所以正个系统是稳定的

唐诗扬-北洋岩土工学生
iLib

Kronex&#039;s Outline of Digital Signals and Syate

iLib

Calculus and Differential Equations with MATLAB

iLib

Applied Speech and Audio Processing - With MATLAB

LOFTER

让兴趣,更有趣

简单随性的记录
丰富多彩的内容
让生活更加充实

下载移动端
关注最新消息