24小时热门版块排行榜    

CyRhmU.jpeg
查看: 211  |  回复: 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的回帖

bib96356

金虫 (正式写手)

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

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的回帖

zjliu1265

至尊木虫 (正式写手)

现在趋势也不对,首先需要保证的是你写的cheshi2.m程序的正确性啊。
你先核查一下吧
4楼2009-09-18 17:03:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见