查看: 144|回复: 4|关注: 0

[已解决] 关于syms与subs函数

[复制链接]

新手

18 麦片

财富积分


050


15

主题

36

帖子

0

最佳答案
syms sigma1 sigma2 sigma3  sigma4 row12 row13 row14 row23 row24 row34;
Cx=[sigma1^2 row12*sigma1*sigma2 row13*sigma1*sigma3 row14*sigma1*sigma4;
   row12*sigma1*sigma2 sigma2^2 row23*sigma2*sigma3 row24*sigma2*sigma4;
   row13*sigma1*sigma3 row23*sigma2*sigma3 sigma3^2 row34*sigma3*sigma4;
   row14*sigma1*sigma4 row24*sigma2*sigma4 row34*sigma3*sigma4 sigma(4)^2];
C2=inv(Cx);
D=det(Cx);
dC2dsigma=[diff(C2,sigma1) diff(C2,sigma2) diff(C2,sigma3) diff(C2,sigma4)];
dDdsigma=[diff(D,sigma1) diff(D,sigma2) diff(D,sigma3) diff(D,sigma4)];
dC2drow=[diff(C2,row12) diff(C2,row13) diff(C2,row14) diff(C2,row23) diff(C2,row24)
     diff(C2,row34)];
dDdrow=[diff(D,row12) diff(D,row13) diff(D,row14) diff(D,row23) diff(D,row24)
     diff(D,row34)];
S=subs(dC2dsigma,{sigma1,sigma2,sigma3,sigma4,row12,row13,row14,row23,row24,row34},
{40 40 40 25 0.2 0.3 0 0.2 0 0});
T=subs(dDdsigma/D,{sigma1,sigma2,sigma3,sigma4,row12,row13,row14,row23,row24,row34},
{40 40 40 25 0.2 0.3 0 0.2 0 0});
有没有大佬帮忙看一下这段程序的问题

回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1629

帖子

342

最佳答案
  • 关注者: 81
发表于 2020-9-11 15:22:16 | 显示全部楼层 |此回复为最佳答案
可以运行了,仅供参考
  1. syms sigma1 sigma2 sigma3  sigma4 row12 row13 row14 row23 row24 row34;
  2. Cx=[sigma1^2 row12*sigma1*sigma2 row13*sigma1*sigma3 row14*sigma1*sigma4;
  3.    row12*sigma1*sigma2 sigma2^2 row23*sigma2*sigma3 row24*sigma2*sigma4;
  4.    row13*sigma1*sigma3 row23*sigma2*sigma3 sigma3^2 row34*sigma3*sigma4;
  5.    row14*sigma1*sigma4 row24*sigma2*sigma4 row34*sigma3*sigma4 sigma4^2];
  6. C2=inv(Cx);
  7. D=det(Cx);
  8. dC2dsigma=[diff(C2,sigma1) diff(C2,sigma2) diff(C2,sigma3) diff(C2,sigma4)];
  9. dDdsigma=[diff(D,sigma1) diff(D,sigma2) diff(D,sigma3) diff(D,sigma4)];
  10. dC2drow=[diff(C2,row12) diff(C2,row13) diff(C2,row14) diff(C2,row23) diff(C2,row24) diff(C2,row34)];
  11.      
  12. dDdrow=[diff(D,row12) diff(D,row13) diff(D,row14) diff(D,row23) diff(D,row24) diff(D,row34)];
  13.      
  14. S=subs(dC2dsigma,{sigma1,sigma2,sigma3,sigma4,row12,row13,row14,row23,row24,row34},...
  15. {40 40 40 25 0.2 0.3 0 0.2 0 0});
  16. T=subs(dDdsigma/D,{sigma1,sigma2,sigma3,sigma4,row12,row13,row14,row23,row24,row34},...
  17. {40 40 40 25 0.2 0.3 0 0.2 0 0});
复制代码
回复此楼 已获打赏: 0 积分

举报

新手

18 麦片

财富积分


050


15

主题

36

帖子

0

最佳答案
 楼主| 发表于 2020-9-11 17:20:33 | 显示全部楼层
20141303 发表于 2020-9-11 15:22
可以运行了,仅供参考

谢谢啦,想问一下我是哪里出错了
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1629

帖子

342

最佳答案
  • 关注者: 81
发表于 2020-9-11 17:42:18 | 显示全部楼层
主要是脚本代码的换行问题
回复此楼 已获打赏: 0 积分

举报

新手

18 麦片

财富积分


050


15

主题

36

帖子

0

最佳答案
 楼主| 发表于 2020-9-11 17:43:07 | 显示全部楼层
20141303 发表于 2020-9-11 17:42
主要是脚本代码的换行问题

明白了!
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表