| 查看: 327 | 回复: 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的等高线图，做不下去了。。。求大侠指教 |
回复此楼» 猜你喜欢
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有9人回复
-
好奇怪的filecode
已经有5人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有10人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12618.2
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
4楼2009-09-18 17:03:41
zjliu1265
至尊木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 12618.2
- 散金: 200
- 帖子: 642
- 在线: 680.6小时
- 虫号: 585320
- 注册: 2008-08-08
- 性别: GG
- 专业: 光学
2楼2009-09-18 07:03:53
3楼2009-09-18 12:52:53
100