| 查看: 348 | 回复: 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的等高线图，做不下去了。。。求大侠指教 |
回复此楼» 猜你喜欢
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
3楼2009-09-18 12:52:53
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12623.7
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
2楼2009-09-18 07:03:53
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12623.7
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
4楼2009-09-18 17:03:41
30