[求助]分段函数画图问题

gaochi526

大家好，小弟是新手，特来请教。
当0<s<1.752时，y=(Log(1-0.062879s+0.032498s²-0.00754s³))/-s；
当1.752<s<3.626时，y=(Log(1-0.062879s+0.032498s²-0.00754s³+ 0.009304 (s-1.752)³))/-s.
这个分段函数怎么画啊，急盼答复，谢谢！

tianj1234

只看效果么？
分别在(0,1.752)和(1.752,3.626)区间中离散点，例如区间1中s选取为0.002:0.002:1.750，此时计算出y结果进行画图即可，同样的方式对应区间2

happy

1. s1=0.008:0.008:1.752;
   
2. s2=1.752:0.002:3.626;
   
3. y1=(Log10(1-0.062879.*s1+0.032498.*s1.^2-0.00754.*s1.^3))./(-s1);
   
4. y2=(Log10(1-0.062879.*s2+0.032498.*s2.^2-0.00754.*s2.^3+0.009304.*(s2-1.752).^3))./(-s2);
   
5. s=[s1 s2];
   
6. y=[y1 y2];
   
7. plot(s,y)

复制代码

ike944

1.752处的函数值是不是需要单独处理？应该去掉吧，题设上没有这点

ganwu

t1=0.8:0.00001:0.808;
>> y1=0;
>> t2=0.808:0.00001:0.86;
>> y2=sin(19.23*pi*(t2-0.808));
>> t=[t1 t2];
>> y=[y1 y2];
>> plot(t,y)
这个为什么不能显示呀，
提示??? Error using ==> plot
Vectors must be the same lengths.
怎么才能解决这个问题，还有这是个重复性的图形，怎么才能简单的实现以后重复性的图形呢，就是好说在时间周轴上重复的出现这个分段函数图形

hyacinth

1. t1=[0.8,0.808];
   
2. y1=[0,0];
   
3. t2=0.808:0.00001:0.86;
   
4. y2=sin(19.23.*pi.*(t2-0.808));
   
5. t=[t1 t2];
   
6. y=[y1 y2];
   
7. plot(t,y)

复制代码

ganwu

实在太谢谢了，能否对我后面提的问题给一些建议呢？我实在是等着急用，我是新手，实在是没有经验

happy

实在太谢谢了，能否对我后面提的问题给一些建议呢？我实在是等着急用，我是新手，实在是没有经验

不懂你上面说的是什么意思？

happy

1. t1=[0.8,0.808];
   
2. y1=[0,0];
   
3. t2=0.808:0.001:0.86;
   
4. y2=sin(19.23.*pi.*(t2-0.808));
   
5. t=[t1 t2];
   
6. y=[y1 y2];
   
7. 
   
8. n=600; %你要画的周期数
   
9. tt=[];
   
10. yy=[];
    
11. for i=1:n
    
12. tt=[tt t+(i-1)*0.06];
    
13. yy=[yy y];
    
14. end
    
15. plot(tt,yy)

复制代码

画那么多周期根本就没有意义

vancross

t1=0.8:0.00001:0.808;
>> y1=0;
此处y1的长度没有t1长，即y1跟t1没建立关系

vancross

原帖由 happy 于 2006-3-19 11:05 发表
t1=[0.8,0.808];y1=[0,0];t2=0.808:0.001:0.86;y2=sin(19.23.*pi.*(t2-0.808));t=[t1 t2];y=[y1 y2];
n=600; %你要画的周期数tt=[];yy=[];for i=1:n    tt=[tt t+(i-1)*0.06];    yy=[yy y];endplot(tt,yy)画那么多 ...

有了周期是为了在后面能够对周期函数做傅里叶变换，然后求频谱图，相位图

因为既然这个周期函数有一段是0，那么必然涉及到占空比，从而可以求对此谐波叠加图

[ 本帖最后由 ChaChing 于 2009-2-14 20:22 编辑 ]