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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab
91/1返回列表
查看: 892  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Fnai

银虫 (初入文坛)

[求助] matlab 已有2人参与

万能的虫友们，理想中要仿真出delta1的曲线，现实是delta1的曲线出不来，拜托大家看看问题出在哪里，谢谢大家
theta=0:0.17:pi;
delta1=0.01:0.01:0.99;
K=length(theta);
N=length(delta1);
y=zeros(1,N);
z=zeros(1,N);
zeta=zeros(1,K);
f=zeros(1,K);
delta=zeros(1,K);
for i=1:K
    zeta(i)=(1-(4*r*cos(theta(i))./(1+h^2+r^2+2*r*cos(theta(i))))).^(alpha/2);
    f(i)=-log(zeta(i));
    for j=1:N
        y(j)=(log(2)/delta1(j)-log(2)/(1-delta1(j))-2*log(delta1(j))+2*log(1-delta1(j)));
       if y(j)-f(i)==0.0001
         delta=delta1(j);
       end
    end
end
plot(delta);
hold on
xlabel('\theta')
ylabel('\deltaA')
回复此楼

» 猜你喜欢
秀而不媚。寒而不清
1楼 2017-03-12 20:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我怀念的TQ

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
delta=delta1(j); 你这个语句最后得到的delta不就只是一个数么，连不成线
赞一下 回复此楼

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

Fnai
MindandHand
2楼2017-03-13 08:48:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Fnai

银虫 (初入文坛)
送红花一朵
引用回帖:
2楼: Originally posted by 我怀念的TQ at 2017-03-13 08:48:35
delta=delta1(j); 你这个语句最后得到的delta不就只是一个数么，连不成线

那怎么将得到的delta1输出来呢

发自小木虫Android客户端
赞一下 回复此楼
秀而不媚。寒而不清
3楼2017-03-13 09:36:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我怀念的TQ

木虫 (著名写手)
引用回帖:
3楼: Originally posted by Fnai at 2017-03-13 09:36:03
那怎么将得到的delta1输出来呢
...

我没有细看你的程序到底是做什么，只是单纯看了语句问题，delta是你最后要得到的曲线，在循环里你应该是想得到它每一次循环对应的值，那么应该是delta（j）=delta1(j);类似这种才可以得到曲线
赞一下 回复此楼

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

Fnai
MindandHand
4楼2017-03-13 09:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Fnai

银虫 (初入文坛)
送红花一朵
引用回帖:
4楼: Originally posted by 我怀念的TQ at 2017-03-13 09:57:06
我没有细看你的程序到底是做什么，只是单纯看了语句问题，delta是你最后要得到的曲线，在循环里你应该是想得到它每一次循环对应的值，那么应该是delta（j）=delta1(j);类似这种才可以得到曲线...

大神，我用的是穷举搜索法，但是好像没用对。我想根据f(i)，在y(j)中找到与f(i)相似的数，然后并输出那个数，可是图形出来的不对，大神，你再帮我看看吧，谢谢啦

发自小木虫Android客户端
赞一下 回复此楼
秀而不媚。寒而不清
5楼2017-03-13 10:08:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqinly

专家顾问 (著名写手)

至人无己，神人无功，圣人无名

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
5楼: Originally posted by Fnai at 2017-03-13 12:08:44
大神，我用的是穷举搜索法，但是好像没用对。我想根据f(i)，在y(j)中找到与f(i)相似的数，然后并输出那个数，可是图形出来的不对，大神，你再帮我看看吧，谢谢啦
...

这位哥不是说了嘛，改成delta（j）=delta1(j);

发自小木虫Android客户端
赞一下 回复此楼
学而时习之，不亦说乎？有朋自远方来，不亦乐乎？人不知而不愠，不亦君子乎？
6楼2017-03-13 10:25:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Fnai

银虫 (初入文坛)
引用回帖:
6楼: Originally posted by qqinly at 2017-03-13 10:25:39
这位哥不是说了嘛，改成delta（j）=delta1(j);
...

不好意思，本人比较愚笨，改完后图形还是不对

发自小木虫Android客户端
赞一下 回复此楼
秀而不媚。寒而不清
7楼2017-03-13 10:27:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我怀念的TQ

木虫 (著名写手)
引用回帖:
5楼: Originally posted by Fnai at 2017-03-13 10:08:44
大神，我用的是穷举搜索法，但是好像没用对。我想根据f(i)，在y(j)中找到与f(i)相似的数，然后并输出那个数，可是图形出来的不对，大神，你再帮我看看吧，谢谢啦
...

不是大神。对于你的程序，首先，h是啥？其次，alpha用法你用对了么？这些问题你得自己找，因为不出结果的原因可能有很多，我只能提及我看到的。如果自己实在找不出来问题，就找个旁边的朋友 同学问问吧。
赞一下 回复此楼
MindandHand
8楼2017-03-13 10:36:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Fnai

银虫 (初入文坛)
谢谢

发自小木虫Android客户端
回复此楼
秀而不媚。寒而不清
9楼2017-03-13 10:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Fnai 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭