24小时热门版块排行榜    

查看: 344  |  回复: 3
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 bib96356 的 15 个金币

bib96356

金虫 (正式写手)

[交流] 【求助】绝望中的求助

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的等高线图,做不下去了。。。求大侠指教
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjliu1265

至尊木虫 (正式写手)

★ ★ ★ ★ ★ ★
wuguocheng(金币+1,VIP+0): 谢谢你的程序 9-18 07:21
bib96356(金币+5,VIP+0): 9-18 12:53
用下面这段程序调用你写的cheshi2函数就可以画出等高线了
jj=1:20;
kk=1:30;
for m=1:length(jj);
    for n=1:length(kk);
        NN(m,n)=cheshi2(jj(m),kk(n));
    end
end
contour(kk,jj,NN)
2楼2009-09-18 07:03:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bib96356

金虫 (正式写手)

做出来的图不是我想要的啊,麻烦在帮忙看看啊
3楼2009-09-18 12:52:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjliu1265

至尊木虫 (正式写手)

现在趋势也不对,首先需要保证的是你写的cheshi2.m程序的正确性啊。
你先核查一下吧
4楼2009-09-18 17:03:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bib96356 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-31 00:32 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-31 00:21 by hZiFeudZyoGR
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:57 by hZiFeudZyoGR
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:53 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:07 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-30 17:56 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见