应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » MATLAB初学者，这样的式子怎么就是画不了图，总是出错!求指导!!!!!!!!!!!!!
51/1返回列表
查看: 2062  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhiy5156

木虫 (正式写手)

[求助] MATLAB初学者，这样的式子怎么就是画不了图，总是出错!求指导!!!!!!!!!!!!!

syms   t;
k1=8490;
k2=2680;
w=12560000;
R12=24.054;  
R21=0.0416;
D=3/8-1/8*sin(w*t);
s1=k2*D;s2=w*t-1/2*k1*D;
fenzi=2000*{[2*cos(s1).*cos(s2)+(R12+R21)*sin(s1).*sin(s2)]+j*[2*cos(s1).*sin(s2)-(R12+R21)*sin(s1).*cos(s2)]};
fenmu=4*[cos(s1)]^2+(R12+R21)^2*[sin(s1)]^2;
y=fenzi/fenmu;
plot(t,y)
总是出错，很郁闷，跪求指导！！！！！！！！！！！！！！

[ Last edited by ben_ladeng on 2011-10-24 at 12:53 ]
回复此楼

» 收录本帖的淘帖专辑推荐

matlab学习精华

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-23 17:33:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

纳斯茵

银虫 (小有名气)
★ ★
微尘、梦想(金币+2): 谢谢参与应助，欢迎常来 2011-10-24 17:50:45
clear;
clc;
%%syms   t;
t=[-100:1:100];
k1=8490;
k2=2680;
w=12560000;
R12=24.054;  
R21=0.0416;
D=3/8-1/8*sin(w*t);
s1=k2.*D;
s2=w*t-1/2*k1*D;
fenzi=2.*cos(s1).*cos(s2)+(R12+R21).*sin(s1).*sin(s2)+j*[2.*cos(s1).*sin(s2)-(R12+R21).*sin(s1).*cos(s2)];
fenmu=4.*[cos(s1)].^2+(R12+R21).^2.*[sin(s1)].^2;
y=(2000*fenzi)./fenmu;
%%figure;ppatterns(t,y);
plot(t,y);

新手哦，不知道符不符合你要的要求哦~
赞一下(1人) 回复此楼
Waiting
6楼2011-10-24 15:13:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

_aqing

银虫 (正式写手)

【答案】应助回帖

★
xzhdty(金币+1): 欢迎常来 2011-10-23 21:45:27
zhiy5156(金币+2): 非常感谢，现在懂了，感谢！！！！！ 2011-10-24 11:15:15
function kk
clc
t=0.1:200
k1=8490;
k2=2680;
w=12560000;
R12=24.054;  
R21=0.0416;
D=3/8-1/8*sin(w.*t);
s1=k2*D;
s2=w.*t-1/2*k1*D;
fenzi=2000.*((2.*cos(s1).*cos(s2)+(R12+R21).*sin(s1).*sin(s2))+i*(2*cos(s1).*sin(s2)-(R12+R21).*sin(s1).*cos(s2)));
fenmu=4*(cos(s1)).^2+(R12+R21).^2*(sin(s1)).^2;
y=fenzi./fenmu;
plot(t,real(y),t,imag(y))
用这个做m文件试试
赞一下(1人) 回复此楼
hehe
2楼2011-10-23 19:14:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★
jjdg(金币+1): 感谢参与 2011-10-24 01:17:02
ben_ladeng: 专家考核存档 2011-10-24 12:53:27
问问题最好把报错信息一并贴上,否则还得开个matlab给你跑代码

1.fenzi这一行,大括号{}在matlab是构造cell数组,你是想要用()吧,改了
2.plot不能plot符号变量,去掉倒数第二行最后的分号,可以在控制台输出y的函数式
CODE:
y =
(4000*cos(12560000*t + (4245*sin(12560000*t))/8 - 12735/8)*cos(1005 - 335*sin(12560000*t)) + (240956*sin(12560000*t + (4245*sin(12560000*t))/8 - 12735/8)*sin(1005 - 335*sin(12560000*t)))/5 - (240956*i*cos(12560000*t + (4245*sin(12560000*t))/8 - 12735/8)*sin(1005 - 335*sin(12560000*t)))/5 + 4000*i*sin(12560000*t + (4245*sin(12560000*t))/8 - 12735/8)*cos(1005 - 335*sin(12560000*t)))/(4*cos(1005 - 335*sin(12560000*t))^2 + (2553496741556419*sin(1005 - 335*sin(12560000*t))^2)/4398046511104)
赞一下(2人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2011-10-23 23:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★
jjdg(金币+1): 感谢参与 2011-10-24 01:16:56
ben_ladeng: 专家考核存档 2011-10-24 12:53:37
引用回帖:
3楼: Originally posted by libralibra at 2011-10-23 23:50:18:
问问题最好把报错信息一并贴上,否则还得开个matlab给你跑代码

1.fenzi这一行,大括号{}在matlab是构造cell数组,你是想要用()吧,改了
2.plot不能plot符号变量,去掉倒数第二行最后的分号,可以在控制台输出y的函数 ...

3.如果非要plot,就不用t做符号,像2楼给出的建议,命名t为一个矩阵
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

zhiy5156
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2011-10-23 23:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 15:18 by KM5EcsNQRBPn
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 14:30 by KM5EcsNQRBPn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
信息提示
关闭
请填处理意见
关闭