| 查看: 286 | 回复: 3 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 bib96356 的 15 个金币 | |||
[交流]
【求助】绝望中的求助
|
|||
|
function minpsi=cheshi2(j,k) phi1=0:1:90; phi2=0; bigphi0=0:1:90; for i=1:length(bigphi0) bigphi=bigphi0(i); B(1,1,:,i)=cosd(phi1)*cosd(phi2)-sind(phi1)*sind(phi2)*cosd(bigphi); B(1,2,:,i)=sind(phi1)*cosd(phi2)+cosd(phi1)*sind(phi2)*cosd(bigphi); B(1,3,:,i)=sind(phi2)*sind(bigphi); B(2,1,:,i)=-cosd(phi1)*sind(phi2)-sind(phi1)*cosd(phi2)*cosd(bigphi); B(2,2,:,i)=-sind(phi1)*sind(phi2)+cosd(phi1)*cosd(phi2)*cosd(bigphi); B(2,3,:,i)=cosd(phi2)*sind(bigphi); B(3,1,:,i)=sind(phi1)*sind(bigphi); B(3,2,:,i)=-cosd(phi1)*sind(bigphi); B(3,3,:,i)=cosd(bigphi); end grain1=[2 -1 1;1 1 -1;0 1 1]; %grain1的thita角和psi角 n0=[1,1,1]; n0=n0/norm(n0); N0=n0*grain1; N0=N0/norm(N0); thita0=atan(-N0(1)/N0(3)); psi0=atan(-N0(1)/N0(2)); n1=[1,1,1]; %grain2的thita角和psi角 n1=n1/norm(n1); N1=n1*B(:,:,j,k); N1=N1/norm(N1); psi(1)=atan(-N1(1)/N1(2)); n2=[1,-1,1]; n2=n2/norm(n2); N2=n2*B(:,:,j,k); N2=N2/norm(N2); psi(2)=atan(-N2(1)/N2(2)); n3=[1,1,-1]; n3=n3/norm(n3); N3=n3*B(:,:,j,k); N3=N3/norm(N3); psi(3)=atan(-N3(1)/N3(2)); n4=[-1,1,1]; n4=n4/norm(n4); N4=n4*B(:,:,j,k); N4=N4/norm(N4); psi(4)=atan(-N4(1)/N4(2)); for h=1:4 bigpsi(h)=psi(h)-psi0; end minpsi=min(bigpsi); minpsi 目标是以j和k分别为x和y轴，作关于minpsi的等高线图，做不下去了。。。求大侠指教 |
回复此楼» 猜你喜欢
今年审到国自然15份，谈谈感受
已经有20人回复
-
Sci. Bull. 悲剧经验
已经有5人回复
-
国自然评分
已经有4人回复
-
提交了我也来说说感想
已经有4人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有11人回复
-
评审有感
已经有22人回复
-
反应很差，大量原料没有反应
已经有6人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12147.2
- 散金: 200
- 帖子: 642
- 在线: 679.9小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
2楼2009-09-18 07:03:53
3楼2009-09-18 12:52:53
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12147.2
- 散金: 200
- 帖子: 642
- 在线: 679.9小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
4楼2009-09-18 17:03:41
5