数模论坛

 找回密码
 注-册-帐-号
搜索
热搜: 活动 交友 discuz
查看: 2929|回复: 7

matlab程序作图中的一个小问题(求教)

[复制链接]
发表于 2004-6-1 03:01:29 | 显示全部楼层 |阅读模式
<>本人在上数学实验时,有如下程序</P>
<>clear all;
r=100;
li=300;
w=4;
%jiao jiasudu;
for i=1:13
    q=pi/12*(i-1);
      b1=(sin(q))^3-sin(2*q)*cos(q);
      b2=r/li*sin(q)+r^3/(2*li^3)*b1;
    b(i)=-w*w*b2;
end
for i=1:13
     q(i)=pi/12*(i-1);
end
plot(q,b)</P>
<>运行得到期望的数据结果和图形,可是如果改成下面的程序就出现问题,求教各位高手</P>
 楼主| 发表于 2004-6-1 03:02:55 | 显示全部楼层
clear all;
r=100;
li=300;
w=4;
%jiao jiasudu;
for i=1:13
    q=pi/12*(i-1);
      b1=vpa((sin(q))^3-sin(2*q)*cos(q),6);
      b2=r/li*sin(q)+r^3/(2*li^3)*b1;
    b(i)=-w*w*b2;
end
for i=1:13
     q(i)=pi/12*(i-1);
end
plot(q,b)
 楼主| 发表于 2004-6-1 03:05:22 | 显示全部楼层
<>其中,b1=vpa((sin(q))^3-sin(2*q)*cos(q),6)与上面的程序不同,我认为所用的函数只是将结果的精度做了限制,可为什么会出错?</P><>有如下错误提示:??? Error using ==&gt; plot
                   Conversion to double from sym is not possible.</P><>求助</P>[em08]
发表于 2004-6-1 04:43:49 | 显示全部楼层
clc;clear all;
r=100;
li=300;
w=4;
%jiao jiasudu;
for i=1:13
    q=pi/12*(i-1);
      b1=vpa((sin(q))^3-sin(2*q)*cos(q),6);
      
      b2=r/li*sin(q)+r^3/(2*li^3)*b1;
      
    b(i)=-w*w*b2;
   
end
for i=1:13
     q(i)=pi/12*(i-1);
end
bn=sym2poly(b)
plot(q,bn)
发表于 2004-6-1 04:45:16 | 显示全部楼层
<>b是sym object</P><>用bn=sym2poly(b)转化就行拉</P><>改后如上
</P>
 楼主| 发表于 2004-6-2 00:38:03 | 显示全部楼层
上面的文字中,sym object是什么意思?
顺便问一下,如何在matlab里查看函数的说明
请赐教!!
不胜感谢!
发表于 2004-6-2 05:04:02 | 显示全部楼层
<>sym object是符号对象</P><>查看说明用help yourfunction或doc yourfunction</P>
发表于 2004-7-13 21:51:22 | 显示全部楼层
5楼的朋友,你确定你没错吧[em01][em02]
您需要登录后才可以回帖 登录 | 注-册-帐-号

本版积分规则

小黑屋|手机版|Archiver|数学建模网 ( 湘ICP备11011602号 )

GMT+8, 2024-11-27 14:39 , Processed in 0.055917 second(s), 18 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表