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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » matlab程序， 被调用函数中使用 if 语句出错
51/1返回列表
查看: 1053  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者hylb3020将赠送您 30 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hylb3020

木虫 (正式写手)

流体力学爱好者

[求助] matlab程序， 被调用函数中使用 if 语句出错 已有1人参与

clear all; close all;
l=45;
alpha1=-0.25;
lambda1=0.25;
n=128;
x=l/n*[-n/2:n/2-1]; y=x;
kx=(2*pi/l)*[0:n/2-1 -n/2:-1];
kx(1)=1e-6; ky=kx;
[x,y]=meshgrid(x,y);
[kx,ky]=meshgrid(kx,ky);
h=2*pi/n;

column=[00.5*(-1).^(1:n-1).*cot((1:n-1)*h/2)]';
d=(2*pi/l)*toeplitz(column,column([1n:-1:2]));
column=[-pi^2/(3*h^2)-1/6-0.5*(-1).^(1:n-1)./sin(h*(1:n-1)/2).^2];
d2=(2*pi/l)^2*toeplitz(column);
column=[0 (-1).^(1:n-1).*cot((1:n-1)*h/2).*(-pi^2/(2*h^2)+3/4*csc((1:n-1)*h/2).^2)]';
w=1+0.1*cos(0.31*x)+0.1*cos(0.31*y);%初始条件

t=[0 17 17.5 18 19 20 21 22 23 25 30 40]

[t,wsol]=ode45('advection_diffusion', t, w( : ),[],n,d,d2,alpha1,lambda1,x);


advection_diffusion.m
function dw=advection_diffusion(t,w,dummy,n,d,d2,alpha1,lambda1, x);
w=reshape(w,n,n); wt=fft2(w);
if (t<=20)
dw=reshape(i*(alpha1*d2*w+lambda1*w*d2')+t,n^2,1);
else
dw=reshape(i*(alpha1*d2*w+lambda1*w*d2')-t,n^2,1);
end


我想求解一个偏微分方程，变量t 取了12个值，在被调用的函数 advection_diffusion 中需要用到一个 if 语句，对 t进行分类，程序一直运行不出来。希望得到大家的指点，不胜感激。
回复此楼

» 猜你喜欢
人生得意须尽欢
1楼 2021-06-29 17:17:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hylb3020

木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by 独孤神宇 at 2021-06-29 19:41:51
if 后面t&lt;20的括号去掉

4楼2021-06-29 22:08:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
if 后面t<20的括号去掉

发自小木虫Android客户端
赞一下(1人) 回复此楼

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

hylb3020 hylb3020
数值计算
2楼2021-06-29 19:41:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hylb3020

木虫 (正式写手)

流体力学爱好者
送红花一朵
引用回帖:
2楼: Originally posted by 独孤神宇 at 2021-06-29 19:41:51
if 后面t&lt;20的括号去掉

呀，确实是这个原因。感谢啊  很久没有上小木虫了，也不知道怎样把金币给你。
赞一下 回复此楼
人生得意须尽欢
3楼2021-06-29 22:06:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +13 xhuama 2026-03-02 14/700 2026-03-05 00:04 by 雪霓裳
[考研] 085600材料调剂 总分330 +4 池池丶 2026-03-03 4/200 2026-03-04 21:15 by wutongshun
[考研] 080500材料科学与工程 +8 202114020319 2026-03-03 8/400 2026-03-04 20:54 by hmn_wj
[考研] 085600材料与化工 298 调剂 +8 小西笑嘻嘻 2026-03-03 8/400 2026-03-04 20:53 by hmn_wj
[考研] 286 +6 ksncj 2026-03-04 6/300 2026-03-04 20:49 by lature00
[考研] 0703化学求调剂，总分320分，一志愿中南民族大学 +3 19286712440 2026-03-04 3/150 2026-03-04 18:30 by houyaoxu
[考研] 0703化学 一志愿东北大学 279分求调剂 +4 聂乐平 2026-03-03 4/200 2026-03-04 16:31 by zhukairuo
[考研] 0703化学求调剂 +3 将未有期 2026-03-04 3/150 2026-03-04 16:19 by sslc1985
[考研] 0703化学306调剂 +4 26要上岸 2026-03-03 4/200 2026-03-03 20:08 by yoyohj0427
[考研] 290求调剂 +9 ErMiao1020 2026-03-02 9/450 2026-03-03 18:03 by linlonghao
[考研] 266材料化工求调剂 +3 哇塞王帅 2026-03-03 3/150 2026-03-03 17:58 by houyaoxu
[考研] 289求调剂 +7 BrightLL 2026-03-02 9/450 2026-03-03 16:02 by tgxtgxtgx9
[考研] 268求调剂 +6 好运连绵不绝 2026-03-02 6/300 2026-03-03 13:03 by 秋收
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 16/800 2026-03-03 00:36 by 好好好1233
[考研] 0856材料调剂 +5 沿岸有贝壳OUC 2026-03-02 5/250 2026-03-02 20:31 by hypershenger
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 哈工大计算机刘劼团队招生 +4 hit_aiot 2026-03-01 6/300 2026-03-02 11:53 by 一声问好
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪！？！
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[基金申请] 刚录用，没有期刊号，但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
信息提示
关闭
请填处理意见
关闭