24小时热门版块排行榜    

查看: 1910  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lwjcug

金虫 (小有名气)

[求助] matlab中画a=1.96直线线求助

小弟进行一个突变检验(0.05水平),其中要画两条a=1.96和a=-1.96的置信区间线。如下:
a=1.96;
a2=-1.96;
year=1981:2011;
plot(year,a,'k-','LineWidth',1);
plot(year,a2,'k-','LineWidth',1);

但是画出来的直线一直是点,都是单个的点,成不了线,我使用的是实线参数,而调整线宽、颜色参数都不管用
但是如果plot(year,a2,'k*','LineWidth',1);换为符号型参数,就可以调其他的
现在需要画实线或虚线点型的不好看
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

学术蠢驴
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军


dubo: 金币+1, 欢迎讨论 2012-06-22 21:57:51
引用回帖:
7楼: Originally posted by lwjcug at 2012-06-22 11:23:36
哥们,追加一问,如本贴图,怎么把那个X轴移到中间y=0的位置啊,现在的样子是X轴在下面。我还是菜鸟...

行是行,就是好麻烦啊.
我索性给你把x,y全部移动到了中心,变成了四象限那种图,
如果需要移动x,y轴位置,根据我的注释,自己改动下数值吧
CODE:
% original data
year = 1981:2011;
a1 = 1.96 * ones(1,length(year));
a2 = -1.96 * ones(1,length(year));

% plot data
plot(year,a1,'k',year,a2,'k');

% xlimit and y limit
xlim([1981,2011]);
ylim([-3,3]);

% hold
hold on;

% get original ticks
X = get(gca,'Xtick');
Y = get(gca,'Ytick');

% get original labels
XL = get(gca,'XtickLabel');
YL = get(gca,'YtickLabel');

% new x-y limits
X_Lim = get(gca,'XLim');
Y_Lim = get(gca,'YLim');

% new position of y-axis
x_offset = diff(X_Lim)/2+X_Lim(1);

% plot new x axis
plot(X_Lim,[0 0],'k','linewidth',2);

% plot new y axis
plot([x_offset x_offset],Y_Lim,'k','linewidth',2);

% small triangles
ax = [X_Lim(end),X_Lim(end)-diff(X_Lim)/100,X_Lim(end)-diff(X_Lim)/100;0,diff(Y_Lim)/100,-diff(Y_Lim)/100];
ay = [x_offset,x_offset+diff(X_Lim)/120,x_offset-diff(X_Lim)/120;Y_Lim(end),Y_Lim(end)-diff(Y_Lim)/80,Y_Lim(end)-diff(Y_Lim)/80];
fill(ax(1,:),ax(2,:),'k');
fill(ay(1,:),ay(2,:),'k');

% get ticks offsets
Xoff = diff(X_Lim)./80+x_offset;
Yoff = diff(Y_Lim)./60;

% draw new ticks
for i=1:length(X)
    plot([X(i) X(i)],[0 Yoff],'-k');
end;
for i=1:length(Y)
    plot([Xoff, x_offset],[Y(i) Y(i)],'-k');
end;

% draw new tick labels
text(X-1,zeros(size(X))-2.*Yoff,XL);
text(Xoff+ones(size(Y))*0.5,Y,YL);

% hide original axis
axis off;

% set white background
set(gcf,'color','w');

结果:
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
8楼2012-06-22 17:38:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
lwjcug: 金币+5, ★★★★★最佳答案 2012-06-22 11:20:07
dubo: 金币+1, 欢迎讨论 2012-06-22 21:57:17
兄弟用line,plot是按照你的数据点plot的
CODE:
line([1981;2011],[1.96;1.96],'Color','k','LineWidth',1);
hold on;
line([1981;2011],[-1.96;-1.96],'Color','k','LineWidth',1);
ylim([-3,3]);

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2012-06-21 17:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖


dubo: 金币+1, 欢迎讨论 2012-06-22 21:57:21
非用plot也可,但是你需要将a和a2都扩展到跟year一样长度才行
CODE:
year=1981:2011;
a=1.96*ones(1,length(year));
a2=-1.96*ones(1,length(year));
plot(year,a,'k-','LineWidth',1);
hold on;
plot(year,a2,'k-','LineWidth',1);
ylim([-3,3]);

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2012-06-21 17:33:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yulai860

铁虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
dubo: 金币+1, 欢迎讨论 2012-06-22 21:57:26
用plot 画也可以吧,k后边不加    —     应该就可以啦

[ 发自手机版 http://muchong.com/3g ]
好好学习天天向上
4楼2012-06-22 08:32:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见